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Abstract 

We  examine  the  electromagnetic  scattering  from  an  arbitrarily  shaped  open  cavity  embed¬ 
ded  in  a  perfectly  conducting,  infinite  ground  plane.  The  cavity  is  filled  with  an  arbitrary  linear, 
isotropic,  homogeneous,  lossy  material.  The  field  everywhere  in  the  cavity  interior  and  above  the 
ground  plane  are  expressed  in  terms  of  the  tangential  fields  on  the  cavity  surface  and  aperture.  We 
develop  a  coupled  set  of  three  integral  equations  to  determine  the  tangential  electromagnetic  fields 
on  the  aperture  and  cavity  surface.  The  integral  equations  are  based  on  fundamental  mathematical 
laws  obeyed  by  the  fields  rather  than  on  the  usual  surface  or  volume  equivalence  principles.  The 
unknown  functions  are  the  tangential  electric  and  magnetic  fields  on  the  cavity  aperture,  and  the 
tangential  magnetic  field  on  the  cavity  surface.  The  support  of  the  unknown  tangential  fields  is 
finite. 

We  design  a  moment-method  based  algorithm  to  find  an  approximate  solution  to  the  integral 
equations  for  plane  wave  illumination,  and  implement  it  on  a  digital  computer.  The  algorithm  is 
specialized  for  the  case  of  axisymmetric  cavities,  and  uses  the  same  basis  functions  as  are  used  in 
traditional  body-of-revolution  scattering  codes.  The  unknown  tangential  fields  are  expanded  using 
piecewise-linear  functions  in  the  elevation  plane  and  complex  exponentials  in  the  azimuth  plane. 
Orthogonality  of  the  complex  exponentials  is  exploited  to  reduce  the  size  of  the  moment  method 
matrix. 

The  algorithm  yields  a  well-conditioned  numerical  solution  for  the  tangential  fields.  The 
tangential  fields  obey  the  edge  condition  at  the  aperture  rim.  The  integral  equations  are  uniquely 
solvable  at  frequencies  where  other  integral  equation-based  techniques  admit  spurious  solutions. 
Radar  cross  section  calculations  are  compared  to  experimental  measurements  of  full-scale  physical 
models.  Results  show  that  an  open  cavity  can  serve  as  an  effective  RCS  enhancement  device. 
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Electromagnetic  Scattering  from  a  Cavity  in  a  Ground  Plane:  Theory  and  Experiment 


7.  Introduction 

1.1  Motivation 

Methods  of  radar  cross  section  (RCS)  reduction  have  been  extensively  investigated  over 
the  past  several  decades  and  have  reached  a  state  of  maturity  sufficient  to  find  many  practical 
applications.  On  the  other  hand  methods  of  RCS  enhancement  have  received  only  limited  attention. 
An  example  of  RCS  enhancement  is  the  use  of  a  string  of  comer  reflectors  to  locate  mnways 
by  approaching  aircraft  in  bad  weather.  A  similar  example,  albeit  in  the  optical  regime,  is  the 
placement  of  cat’s  eye  lane  demarkers  on  highways  to  help  motorists  at  night.  Such  aids  represent 
“bmte  force”  approaches  employing  already  existing  technology,  and  are  not  constrained  by  other 
factors  such  as  the  requirement  to  minimize  size  or  weight.  As  such,  they  do  not  constitute  the 
product  of  a  serious,  comprehensive  study  of  RCS  enhancement  techniques.  Hence  it  would  appear 
that  an  investigation  of  RCS  enhancement  techniques  is  warranted.  Before  one  begins  to  undertake 
such  studies,  it  is  worthwhile  to  identify  a  few  practical  applications  to  justify  any  such  efforts. 

•  RPV  Tracking:  Remotely  piloted  vehicles  (RPVs)  are  an  inexpensive  means  of  gathering 
intelligence  data  and  are  routinely  used  in  tactical  scenarios.  It  goes  without  saying  that 
the  user  must  maintain  constant  contact  with  the  RPV.  Since  the  operation  is  supposed  to 
be  clandestine,  the  use  of  tracking  beacons  is  ill-advised.  To  enhance  the  ability  to  track 
the  RPV  by  radar,  RCS  enhancement  without  compromising  flying  qualities  appears  attrac¬ 
tive.  Comer  reflectors  or  Luneberg  lenses  are  devices  that  can  enhance  RPV  RCS,  but  such 
devices  invariably  increase  the  drag  of  an  otherwise  streamlined  air  vehicle.  The  answer 
obviously  lies  in  working  from  the  inside  with  suitable  coupling  to  the  outside  by  means  of 
RF  windows  which  blend  smoothly  with  the  skin  of  the  RPV. 

•  UAV  Decoys:  Unmanned  air  vehicles  (UAVs)  can  serve  as  an  inexpensive  decoy.  Such  de¬ 
coys  might  be  launched  in  large  numbers  to  saturate  enemy  air  defenses,  thus  increasing  the 
survivability  of  friendly  aircraft.  A  small,  inexpensive  UAV  would  typically  have  a  smaller 
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RCS  than  that  of  a  fighter-size  aircraft,  thereby  providing  an  adversary  with  a  way  to  dis¬ 
criminate  between  the  two.  An  effective  decoy  should  mimic  the  strike  aircraft  in  as  many 
ways  as  possible,  so  we  desire  a  method  to  increase  the  decoy’s  RCS  to  the  level  of  the  strike 
aircraft.  In  this  way,  the  effectiveness  of  a  UAV  decoy  could  be  considerably  increased. 

•  Camouflage  of  Stealthy  Aircraft:  As  low-RCS  aircraft  enter  the  inventory,  the  need  to  con- 
trollably  increase  their  RCS  has  become  evident.  In  many  cases,  a  stealthy  aircraft’s  low  RCS 
might  impair  its  capability,  especially  in  a  peacetime  training  environment.  For  instance,  air 
traffic  control  radars  must  be  able  to  detect  low-RCS  aircraft  in  their  areas  of  responsibil¬ 
ity  to  ensure  air  safety;  and  reliance  on  a  transponder  may  introduce  safety-of-flight  issues. 
Additionally,  one  may  wish  to  protect  knowledge  of  a  low-RCS  aircraft’s  true  radar  signa¬ 
ture,  so  the  ability  to  increase  it  in  a  controlled  manner  may  prevent  the  compromise  of  vital 
operational  capabilities. 

•  Identification:  Current  non-cooperative  target  recognition  (NCTR)  research  efforts  exploit 
target  RCS  signatures,  so  it  is  conceivable  that  controllable,  distinctive  RCS  enhancement 
could  be  used  as  a  means  to  identify  friendly  aircraft.  Currently,  friendly  aircraft  are  nor¬ 
mally  identified  via  Identification  Friend  or  Foe  (IFF)  transponders  (which  require  an  air¬ 
craft  to  broadcast  its  presence  using  coded  signals  over  specific  frequencies)  or  through 
non-cooperative  techniques  that  fail  at  some  aspect  angles.  An  RCS  enhancement  device 
may  create  an  RCS  signature  that  is  distinctive  enough  to  make  an  unerring  friendly  identi¬ 
fication  at  all  aspect  angles  possible. 

Several  techniques  have  been  considered  to  accomplish  effective  RCS  enhancement.  Imped¬ 
ance  loading  has  been  investigated  to  control  RCS  [55],  although  most  research  in  impedance 
loading  techniques  has  been  focused  on  RCS  reduction  and,  furthermore,  the  technique  has  been 
found  to  be  effective  only  over  fairly  narrow  bandwidths.  RCS  enhancement  has  traditionally  been 
accomplished  using  comer  reflectors,  bicones,  or  dielectric  lenses  [12,  Chap  8].  These  devices 
are  effective  RCS  enhancement  devices  because  they  scatter  strongly  over  large  bandwidths,  broad 
aspect  angles,  and  varied  polarizations.  However,  when  applied  to  an  air  vehicle,  they  exhibit  the 
following  drawbacks: 

•  Decreased  aerodynamic  performance  through  increased  drag, 

•  Additional  parasitic  weight,  and 
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•  Lack  of  in-flight  switching  capability 

It  would  be  useful  if  an  alternative  method  of  RCS  enhancement  could  be  found  to  mitigate 
these  drawbacks,  while  effectively  increasing  RCS. 


Figure  1  An  open  cavity  embedded  in  the  wing  of  a  UAV 


1.1.1  A  Viable  Candidate  for  RCS  Enhancement.  A  potentially  viable  candidate  design 
for  an  RCS  enhancement  device  consists  of  one  or  more  open  cavities  embedded  behind  the  air 
vehicle’s  exterior  skin.  Figure  1  shows  a  possible  location  of  an  open  cavity  on  an  unmanned 
air  vehicle  (UAV)  -  on  the  upper  surface  of  its  wing.  Here,  the  open  cavity  would  provide  RCS 
enhancement  for  aspects  above  the  waterline  of  the  UAV;  RCS  enhancement  below  the  waterline 
could  be  accomplished  by  a  similar  open  cavity  mounted  on  the  wing’s  underside.  Unlike  an  ex¬ 
ternal  radar  reflector,  an  open  cavity  would  not  extend  into  the  airstream  and  create  additional 
aerodynamic  drag.  The  additional  weight  to  the  air  vehicle  is  negligible.  An  open  cavity  could 
conceivably  be  turned  on  or  off  by  electrically  controlling  the  conductivity  of  a  membrane  cov¬ 
ering  its  aperture  thus  providing  an  in-flight  switching  capability  ( [52, 66]).  *  Additionally,  it  is 
anticipated  that  an  open  cavity  would  be  an  effective  RCS  enhancement  device,  scattering  strongly 

*In  the  “off  state,  the  membrane’s  conductivity  would  be  very  high,  blending  the  aperture  of  the  open  cavity  into 
the  conducting  skin  of  the  aircraft  and  effectively  isolating  the  open  cavity  from  the  incident  field.  In  this  state,  the  open 
cavity  would  not  affect  the  normal  RCS  of  the  air  vehicle.  In  the  “on”  state,  the  membrane’s  conductivity  would  be  very 
low,  and  the  membrane  would  be  nearly  transparent  to  the  incident  field.  In  this  state,  the  open  cavity  would  scatter 
strongly,  overriding  the  natural  RCS  of  the  aircraft. 
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over  a  large  bandwidth,  broad  range  of  aspect  angles,  and  arbitrary  polarization;  it  is  the  goal  of 
this  dissertation  to  characterize  its  scattering  characteristics. 

1.1.2  Open  Cavity  RCS  Enhancement  Device.  For  the  purposes  of  this  dissertation,  we 
define  an  open  cavity  RCS  enhancement  device  (OCRED)  to  consist  of  an  open  hox  or  enclosure 
whose  mouth  or  aperture  is  set  in  a  perfectly  conducting  plane.  A  “shoebox”  shaped  example  of 
this  geometry  is  shown  in  Figure  3(a).  The  interior  of  the  open  cavity  couples  to  the  half  space 
above  via  an  opening  (“aperture”)  in  the  ground  plane.  The  ground  plane  is  perfectly  conducting, 
as  are  the  walls  of  the  open  cavity,  and  the  interior  of  the  open  cavity  is  filled  with  a  homogeneous, 
isotropic,  linear  dielectric  material.  The  half-space  above  the  ground  plane  is  a  free-space  region. 
The  scattering  properties  of  an  OCRED  can  be  studied  using  analytical  methods  complemented  by 
experimental  measurements.  There  is  a  woeful  lack  of  experimental  data  in  the  open  literature,  and 
each  analytical/computational  method  in  common  use  is  limited  in  some  way. 

In  the  context  of  this  document,  the  OCRED  is  a  canonical  problem  of  the  more  complex 
situation  in  the  real  world.  For  example,  in  the  real  world  the  OCRED  would  he  mounted  on  a 
finite-sized  air  vehicle  rather  than  in  an  infinite  ground  plane,  and  the  cavity  walls  and  surrounding 
air  vehicle  skin  may  not  be  perfectly  conducting.  We  will  deal  here  with  the  idealized,  canonical 
problem  as  a  means  to  gain  insight  into  the  real  world  one. 

1. 1.2. 1  Analytical  Predictions.  Predicting  the  scattering  characteristics  of  the  OC¬ 
RED  is  a  challenging  problem  in  electromagnetic  theory,  especially  in  the  resonance  regime.  In  the 
low-frequency  and  high-frequency  limits,  certain  approximations  can  be  made  to  make  the  prob¬ 
lem  more  tractable,  but  when  the  dimensions  of  the  OCRED  are  on  the  order  of  a  wavelength,  these 
approximations  generally  fail.  Most  current  general-purpose  RCS  prediction  codes  are  based  on 
either  high-frequency  asymptotic  methods,  or  the  method  of  moments.  High-frequency  asymptotic 
methods  include  physical  optics  (PO),  the  physical  theory  of  diffraction  (PTD),  geometrical  optics 
(GO),  and  the  uniform  theory  of  diffraction  (UTD),  and  view  scattering  as  a  phenomenon  depen¬ 
dent  only  on  the  incident  field  and  local  geometry.  But  in  a  cavity,  the  scattering  is  a  complex 
function  of  the  interactions  between  the  aperture,  the  cavity  interior,  and  the  ground  plane.  For 
this  reason,  high-frequency  asymptotic  methods  do  not  predict  OCRED  scattering  accurately.  The 
conventional,  general-purpose  integral  equation  approach  has  the  inherent  capability  to  predict  OC- 
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RED  scattering,  but  the  need  to  model  the  equivalent  currents  on  the  infinite  ground  plane  makes 
it  impractical.  (An  important  variation,  the  generalized  network  formulation,  will  be  discussed  in 
Chapter  E.) 


1.1. 2. 2  Experimental  Measurements.  Building  and  measuring  an  OCRED  model 
has  the  advantage  of  providing  high-quality  data  and  is  easily  within  the  capabilities  of  today’s 
fabrication  technology  and  RCS  measurement  ranges.  Unfortunately,  it  is  not  cost-effective  to 
build  and  measure  many  such  models  with  the  aim  of  obtaining  insight  into  the  viability  or  design 
of  an  effective  OCRED,  and  there  is  a  surprising  lack  of  measured  data  in  the  open  literature. 

1.2  Organization 

This  dissertation  represents  an  attempt  to  understand  the  electromagnetic  scattering  charac¬ 
teristics  of  the  OCRED  geometry.  In  Chapter  E,  we  summarize  the  pertinent  related  work  that  has 
appeared  in  the  open  literature.  In  Chapter  IE,  we  develop  a  new  set  of  integral  equations  to  predict 
the  scattering  from  an  arbitrary  OCRED  geometry.  In  Chapter  IV,  we  develop  an  algorithm  to  find 
an  approximate  solution  to  these  integral  equations  for  the  case  of  axially  symmetric  (body  of  rev¬ 
olution)  geometries.  In  Chapter  V  we  present  the  results  of  implementing  the  numerical  algorithm 
and  making  experimental  measurements  of  several  physical  models.  In  Chapter  VI,  we  conclude 
the  dissertation  with  a  brief  summary  of  the  findings  and  recommendations  for  further  research. 
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11.  Related  Work 


The  OCRED  geometry  consists  of  two  canonical  geometries:  an  open  cavity  and  an  aperture  in  a 
ground  plane.  Each  of  these  forms  has  been  investigated  extensively,  but  the  composite  problem 
has  received  far  less  attention. 

2.1  Open  Cavities 

A  large  body  of  work  exists  dealing  with  the  problem  of  scattering  from  an  open,  perfectly 
conducting  cavity  in  free  space  [17, 27, 38, 58].  Most  of  this  work  makes  use  of  Harrington  and 
Mautz’  generalized  network  formulation  [23].  This  method  partitions  the  problem  into  an  exterior 
problem  and  an  interior  problem,  as  shown  in  Figure  2.  Each  (simpler)  problem  is  treated  indi¬ 
vidually,  and  the  composite  problem  is  solved  by  enforcing  tangential  field  continuity  across  the 
aperture.  If  there  are  no  sources  present,  the  electric  field  for  the  interior  problem  will  satisfy  the 
homogeneous  vector  wave  equation  and  the  perfect  electrical  conductor  (PEC)  boundary  condi¬ 
tion.  For  an  arbitrarily  shaped  interior  problem,  it  can  be  shown  [11,68]  that  the  interior  fields  are 
not  uniquely  determinable  at  a  countably  infinite  set  of  discrete,  resonant  frequencies.  At  these 
frequencies,  the  generalized  network  formulation  fails  to  provide  a  unique  solution  for  the  fields  in 
the  cavity  interior  and  surface  [35].  This  is  a  very  serious  limitation,  especially  for  cavities  of  arbi¬ 
trary  shape,  since  it  is  very  difficult  to  find  the  resonant  frequencies.  Only  for  a  limited  number  of 
simple  geometries  (such  as  rectangular  boxes,  spheres,  and  circular  cylinders)  one  can  analytically 
determine  the  resonant  frequencies.  It  is  therefore  not  apparent  whether  the  generalized  network 


Figure  2  Partitioning  of  an  open  cavity  problem  into  an  exterior  problem  and  an  interior  prob¬ 
lem,  showing  equivalent  magnetic  current  M. 
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formulation  will  break  down  for  a  particular  frequency  of  interest.  The  resonance  phenomenon  ex¬ 
hibited  by  the  interior  problem  does  not  mean  that  the  OCRED  fields  are  not  uniquely  solvable  at 
all  frequencies.  The  source-free  interior  problem  artificially  enforces  a  vanishing  tangential  electric 
field  across  the  aperture,  whereas  the  complete  OCRED  problem  implies  no  such  requirement.  It  is 
this  “spurious  resonance”  that  motivates  us  to  find  an  alternate  method  to  the  generalized  network 
formulation  technique. 

A  related  problem  is  that  of  scattering  by  perfectly  conducting  finite  obstacles  in  free  space. 
This  problem  is  typically  solved  using  the  equivalence  principle  [21]  which  involves  determining 
equivalent  electric  currents  on  the  boundary,  satisfying  an  integral  equation.  Typically,  the  integral 
equation  is  either  the  electric  field  integral  equation  (EEIE)  or  the  magnetic  field  integral  equation 
(MFIE)  [40].  Neither  formulation  is  uniquely  solvable  at  interior  resonant  frequencies  [57].  The 
root  of  the  problem  is  the  failure  of  the  EFIE  and  MFIE  to  enforce  proper  boundary  conditions 
on  the  scatterer  at  resonance  [51, 68].  This  results  in  a  nontrivial  solution  to  the  homogeneous 
(source-free)  problem,  which  contaminates  the  true  solution.  More  complicated  formulations  have 
been  proposed  to  “fix”  the  EFIE  and  MFIE  [42,43,67,68].  Unfortunately,  none  of  these  approaches 
apply  to  the  interior  problem  because  the  resonances  of  the  interior  cavity  are  real. 

Despite  its  problems  with  “spurious  resonances,”  the  generalized  network  formulation  has 
been  successfully  used  by  many  investigators  (see,  for  example,  [63]).  In  all  these  cases,  particular 
care  is  taken  to  stay  away  from  frequencies  corresponding  to  known  internal  resonances. 

2.1.1  Open  Spherical  Shells.  A  particular  canonical  open  cavity  which  can  be  ana¬ 
lyzed  without  resorting  to  the  genereilized  network  formulation  is  the  open,  infinitesimally  thin, 
conducting  spherical  shell  [5,56,58,69].  Because  this  geometry  is  separable,  it  can  be  analyzed 
by  expanding  the  exterior  and  interior  fields  using  spherical  vector  wave  functions  [47].  The  re¬ 
sulting  infinite  series  is  typically  truncated,  and  a  linear  system  of  equations  solved  to  determine 
the  unknown  coefficients.  Unfortunately,  this  technique  applies  only  to  geometries  natural  to  the 
spherical  coordinate  system. 
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2.2  Apertures  in  a  Ground  Plane 

Scattering  by  and  transmission  through  an  arbitrary  aperture  in  a  ground  plane  has  been  ex¬ 
tensively  investigated.  An  excellent  review  of  techniques  developed  through  the  late  1970s  can 
be  found  in  Butler  et  al  [8].  More  recently,  investigators  have  explored  other  techniques.  Lueb- 
bers  [39]  uses  the  finite-difference  time  domain  (FDTD)  method  to  find  the  scattering  from  an  infi¬ 
nite  slot  in  a  thick  ground  plane.  This  method  can  be  extended  to  three-dimensional  apertures,  but 
the  geometry  must  be  approximated  by  a  rectangular  grid.  The  finite  element  method  (FEM)  has 
become  popular  [16, 26, 28, 29],  but  requires  modeling  the  interior  of  the  computational  domain 
and  matching  the  solution  to  a  boundary  integral  equation  that  implicitly  contains  the  radiation 
condition.  Several  authors  [10, 13]  have  considered  the  problem  of  electrically  small  apertures  and 
expressed  the  solution  as  a  power  series  in  frequency.  When  the  frequency  is  sufficiently  small  the 
first  term  of  the  series  yields  reasonably  accurate  results. 

2.3  Open  Cavity  in  a  Ground  Plane 

The  composite  OCRED  problem  represents  a  combination  of  an  open  cavity  and  an  aperture 
in  a  ground  plane.  The  analogous  two-dimensional  problem  of  scattering  by  an  infinite  trough  in  a 
ground  plane  has  been  treated  by  several  investigators  [3,14,20,24,25].  The  approaches  which  have 
been  taken  can  be  broadly  categorized  as  low-frequency,  high-frequency,  and  resonant-frequency 
techniques,  based  on  the  assumed  value  of  kd  where  k  is  the  wavenumber  and  d  is  a  characteristic 
linear  dimension  of  the  cavity  or  aperture. 

2.3.1  Low-Frequency  Techniques,  kd  1.  Hansen  and  Yaghjian  [20]  examine  the 
two-dimensional  trough  in  a  ground  plane,  using  a  truncated  power  series  in  kd.  This  represents  a 
perturbation  of  the  corresponding  electrostatic  problem.  Kempel  and  Senior  [32]  model  an  OCRED 
as  an  impedance  disk  on  a  ground  plane.  This  technique  can  only  be  applied,  however,  to  cavity 
geometries  that  can  be  modeled  as  sections  of  waveguides.  Since  the  size  of  the  cavity  is  assumed 
to  be  small  compared  to  the  wavelength,  low-frequency  techniques  are  primarily  useful  for  finding 
the  scattering  from  imperfections  such  as  cracks,  gaps,  and  rivet-holes  in  the  skins  of  aircraft.  For 
larger  cavities,  low-frequency  techniques  are  not  useful. 
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2.3.2  High-Frequency  Techniques,  kd  »  1.  Nearly  all  high  frequency  techniques  may 
be  grouped  under  the  term  “ray”  techniques,  because  they  take  the  viewpoint  that  electromag¬ 
netic  waves  propagate  like  rays  of  optical  light.  Ray  techniques  are  valid  if  scattering  is  a  local 
phenomenon  and  physically  distant  portions  of  the  geometry  do  not  significantly  interact.  Many 
investigators  have  used  ray  techniques  and  in  several  instances  merged  them  with  modal  methods 
when  the  cavity  has  properties  of  a  waveguide  [7, 9, 15, 36, 37, 49].  High-frequency  techniques 
have  application  in  the  prediction  of  scattering  from  jet  engine  inlet  ducts.  When  the  cavity  is  not 
electrically  large,  interactions  between  the  aperture  and  cavity  termination  become  important  and 
the  accuracy  of  high-frequency  techniques  is  questionable. 

2.3.3  Resonant-Frequency  Techniques,  kd  =  0{1).  Resonant-frequency  techniques 
employ  the  finite  element,  finite-difference  time  domain,  modal,  or  integral  equation  method,  or 
some  hybridization  of  them.  These  are  discussed  below. 

Hybrid  Finite  Element  Method.  The  finite  element  method  (FEM)  has  been  combined  with 
the  boundary  integral  equation  method  to  solve  OCRED-like  problems  [3].  It  has  been  used  to 
model  arbitrary  geometries  [27].  However,  “penalty  functions”  [53]  must  be  introduced  into  the 
three-dimension£il  FEM  variational  equations  to  suppress  spurious  solutions.  Sun  et  al  show  that 
these  spurious  solutions  are  caused  by  an  incorrect  component  of  the  solution  associated  with  the 
static  electric  potential  (i.e.,  the  null  space  of  the  curl  operator),  rather  than  “the  nonsolenoidal 
nature  of  finite-element-approximation  procedures”  [60].  Webb  [64]  uses  special  basis  functions 
to  model  field  behavior  at  edges,  increasing  the  complexity  of  the  method. 

Finite-Difference  Time  Domain  (FDTD)  Method.  The  FDTD  method  has  been  employed  to 
predict  the  scattering  from  open  cavities.  Penney  and  Luebbers  [50]  use  FDTD  to  analyze  cavity- 
backed  antennas,  while  Lee  et  al  [34]  apply  it  to  find  the  scattering  from  a  jet  engine  inlet  duct. 
A  single  FDTD  simulation  can  generate  broadband  scattering  data,  but  only  for  a  single  angle 
of  incidence.  Additionally,  FDTD  meshes  are  restricted  to  be  rectangular,  and  special  absorbing 
boundary  conditions  must  be  used  to  approximate  the  radiation  condition  for  the  scattered  fields. 

Modal  Methods.  For  rectangular  geometries,  the  fields  in  the  OCRED  cavity  can  be  ex¬ 
panded  in  waveguide  modes  [4].  Large  field  magnitudes  in  the  aperture  have  been  observed  when 
the  aperture  is  electrically  small  and  the  frequency  is  near  a  resonant  frequency  for  the  interior 
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problem  [35].  Modal  methods  are  restricted  to  the  narrow  class  of  geometries  for  which  analytical 
expressions  for  the  modes  are  available. 

Integral  Equation  Methods.  Most  integral  equation  approaches  to  the  OCRED  scattering 
problem  employ  the  generalized  network  formulation,  as  described  in  Section  2.1,  and  are  not 
uniquely  solvable  at  some  discrete  frequencies.  An  interesting  and  very  promising  boundary  inte¬ 
gral  equation  approach  was  recently  proposed  by  Asvestas  &  Kleinman  [2].  Their  method  requires 
solving  for  surface  currents  on  the  open  cavity  surface,  and  the  authors  claim  that  the  integral 
equations  are  uniquely  solvable  at  all  frequencies,  but  provide  no  substantiation.  Additionally,  the 
formulation  only  applies  to  empty  cavities  and  the  theory  must  be  extended  to  account  for  material- 
filled  cavities.  Recently,  the  method  was  extended  to  handle  a  dyadic  impedance  boundary  condi¬ 
tion  imposed  on  the  cavity  surface,  though  the  ground  plane  is  still  perfectly  conducting  [1]. 

To  select  an  appropriate  technique  to  model  OCRED  scattering,  we  look  at  three  desirable 
attributes.  First,  the  technique  must  be  valid  for  cavities  that  are  neither  electrically  large  nor  elec¬ 
trically  small.  Second,  the  technique  should  not  suffer  from  the  problem  of  spurious  resonances. 
Third,  the  technique  should  minimize  the  computational  domain  so  that  computational  resources 
are  conserved.  The  technique  that  best  satisfies  these  three  attributes  is  the  Asvestas  boundary 
integral  equation  technique,  which  will  be  further  explored  in  the  remainder  of  this  dissertation. 
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Ill  Theory 

In  this  chapter,  we  will  cast  the  problem  of  scattering  from  an  open  cavity  embedded  in  a  perfectly 
conducting  ground  plane  in  terms  of  integral  equations.  We  shall  first  define  the  specific  geometry 
for  the  field  components.  Following  the  approach  of  Asvestas  and  Kleinman  [2],  we  shall  then 
obtain  the  final  integral  equations. 

In  Section  3.1,  we  introduce  the  OCRED  geometry.  In  Section  3.2,  we  define  the  electro¬ 
magnetic  fields  for  the  OCRED  problem,  along  with  the  partial  differential  equations  and  boundary 
conditions  they  satisfy.  In  Section  3.3,  we  introduce  several  scalar  and  dyadic  functions  and  de¬ 
fine  their  properties.  In  Section  3.4,  we  develop  radiation  integrals  that  define  the  electromagnetic 
fields  everywhere  in  the  OCRED  domain  in  terms  of  the  tangential  field  components  on  the  cavity 
surface  and  aperture.  Finally,  in  Section  3.5,  we  transform  the  radiation  integrals  into  a  set  of  three 
coupled  surface  integral  equations  involving  these  tangential  field  components  alone. 

3.1  Geometry 

The  general  OCRED  geometry  is  shown  in  Figure  3.  The  entire  a:y-plane  is  partitioned  into 
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(a)  A  shoebox-shaped  OCRED  in  a  finite  ground  (b)  Cross-sectional  view  of  an  OCRED, 

plane.  For  an  infinite  ground  plane,  let  w  -¥  oo. 

Figure  3  The  geometry  of  an  OCRED.  Figure  3(a)  shows  a  three-dimensional  representation. 
Figure  3(b)  shows  a  two-dimensional  view,  with  regions  and  surfaces  defined. 

two  regions:  a  bounded,  simply-connected  region,  cr,  and  its  complement  in  the  xy-plane,  The 
cavity  volume,  D,  is  bounded  by  the  closed  regular  surface  dD,  which  is  composed  of  a  and  the 
open  regular  surface,  S.  In  the  following  development,  we  shall  assume  the  following: 

•  The  cavity  surface  S  and  the  aperture  complement  are  perfect  electric  conductors. 
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•  The  boundary  of  a  (i.e.  the  rim  of  the  aperture),  da,  is  a  closed,  bounded,  piecewise  smooth 
curve  in  the  a;y-plane. 

•  The  cavity  surface  S  does  not  extend  into  the  upper  half  space,  and  intersects  the  xy-plane 
only  at  da. 

•  The  unit  normal  vector,  n,  is  defined  almost  everywhere  (except  for  a  finite  number  of  edges, 
comers,  or  tips)  on  S  and  is  oriented  away  from  D. 

•  The  cavity  volume  D  is  filled  with  a  homogeneous,  isotropic,  linear  material  characterized 
by  its  permittivity  6i  and  permeability  /xi,  one  or  both  of  which  may  be  complex. 

•  The  upper  half  space  is  filled  with  free  space. 

Before  proceeding  further,  we  should  clarify  that  the  symbol  r  will  be  used  to  denote  either 
the  position  vector  xx  +  yy  +  zi  or  the  corresponding  point  (x,  y,  z).  The  distinction  will  be  clear 
from  the  context  in  which  r  is  used. 


Figure  4  Position  vector,  r,  and  its  image,  Fj 


In  what  follows  we  shall  make  extensive  use  of  image  theory.  We  define  the  image  of  r 
with  respect  to  the  xy-plane,  as  fj  =  xx  +  yy  —  zz  =  r  •  Ij,  where  Ij  =  xx  +  yy  —  zz  is  the 
image  dyadic.  Similarly  for  a  field  vector  A(f)  ==  Aa;(f)x  +  Ay(r)y  +  Az(r)z,  we  have  the  image 
Ai(r)  =  Aj;(r)x  +  Ay(r)y  -  Az(f)z  =  A(r)  •  Ij.  To  summarize,  the  image  of  any  vector,  be  it  a 
position  vector  or  a  field  vector  function,  is  simply  that  vector  dotted  into  the  image  dyadic. 

Image  surfaces  and  image  volumes  are  defined  analogously  to  image  vectors:  the  image  of 
the  cavity  surface  S  is  denoted  Si  and  defined  as  Si  =  (r  :  fj  e  S},  while  the  image  of  the  cavity 
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volume  D  is  denoted  Di  and  is  defined  as  Di  —  {r  :  Vi  E  D},  The  portion  of  the  upper  half  space 
less  those  points  that  are  not  contained  in  either  volume  Di  or  surface  Si  is  denoted  Df. 

3.2  Fields 

We  restrict  all  impressed  sources  to  exist  only  in  Df  so  that  the  total  fields  in  D  and  Di 
satisfy  the  following  Maxwell’s  equations  for  homogeneous,  linear,  isotropic,  source-free  media 
time  convention)  [21]: 


V  X  E  ==  —jiVfiM 

(la) 

V  X  H  =  j(jt}eE 

(lb) 

VE  =  0 

(Ic) 

0 

II 

t> 

(Id) 

where  E  =  E(r)  is  the  electric  field  intensity  in  volts/meter,  and  H 

=  H(r)  is  the  magnetic 

field  intensity  in  Amperes/meter,  a;  is  the  frequency  in  radians/second,  /u  is  the  permeability  of 
the  medium  in  Henrys/meter,  e  is  the  permittivity  of  the  medium  in  Farads/meter,  and  j  = 

Taking  the  curl  of  equations  (la)  and  (lb)  and  combining  them,  we  see  that  both  E  and  H  satisfy 
the  homogeneous  wave  equation  [59] 

V  X  V  X  E(f)  -  k^B{r)  =  0  (2a) 

V  X  V  X  H(r)  -  fc2H(r)  =  0  (2b) 

where  k  =  Uy/fle  is  the  propagation  constant  of  the  medium  in  inverse  meters.  We  assume  that 
the  known  incident  fields,  and  impinge  on  the  open  cavity,  giving  rise  to  the  reflected 
fields,  E*®^  and  and  the  diffracted  or  scattered  fields,  E®®®*  and  H®®®*.  The  given  incident  fields 
are  defined  to  be  those  that  exist  in  unbounded  free  space.  The  reflected  fields  are  defined  to  be 
those  scattered  by  an  unbroken,  perfectly-conducting  ground  plane  located  at  2:  =  0.  The  scattered 
fields  represent  a  perturbation  due  to  the  presence  of  the  cavity,  and  are  defined  as  the  remainders: 
gscat  =  E  -  E‘"®  -  E'®f  and  H®®®‘  =  H  -  H‘"®  -  H'®^ 
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In  addition  to  equations  (1)  and  (2),  the  fields  must  also  satisfy  certain  boundary  conditions. 
The  tangential  component  of  the  total  electric  field  must  vanish  on  and  S: 


fi  X  E(r)  =  0 

for  r  G  S' 

(3a) 

z  X  E(r)  =  0 

for  r  G 

(3b) 

The  tangential  components  of  the  total  electric  and  magnetic  fields  must  be  continuous  across  the 
aperture  a: 


lim  [z  X  E(r -|-  J  z)  —  z  x  E(r  —  ^  z)l  =0 

for  r  G  cr 

(4a) 

Im  [z  X  H(f  +  z)  -  z  X  H(r  -  5  z)]  =  0 

for  r  G  a 

(4b) 

Finally,  the  scattered  fields  must  satisfy  the  Sommerfeld’s  radiation  conditions  at  infinity  ( [59,  pp. 
485-6]  and  [31,  p.  56]) 


lim 

|r|— >c»  L 

lim 

|f|->00 


F.  +  jifc|r|E*“‘(r)]  =  0 

r  •  VH*=®(r)  -f  jA;ir|H'“*(r)]  =  0 


as  well  as  the  regularity  conditions  [59,  pp.  485-6] 


(5a) 

(5b) 


lim  ir|E*‘^‘“(r)  =  0 

(6a) 

|r|— >oo 

Urn  |r|H*«“(r)  =  0 

r|— >00 

(6b) 

3. 3  Green ’s  Functions 

We  introduce  the  scalar  Green’s  function 


G(fc;r,r')  = 


-jkR 


R=\r-r'\ 


(7) 


where  r  =  xx  +  yy  +  zz  and  r'  =  x'x  -I-  y'y  -|-  z'z  are  position  vectors  in  unprimed  and 
primed  coordinates,  respectively.  The  scalar  Green’s  function  satisfies  the  distributional  differential 
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equation 


V^G  +  k'^G  =  -S{r-r') 

=  -5{x  -  x')  5{y  -  y')  6{z  -  z') 


(8) 


as  well  as  the  scalar  analogues  of  the  radiation  and  regularity  conditions  expressed  by  equations  (5) 
and  (6).  For  readability,  we  will  suppress  the  dependence  of  G{k-,  r,  r ')  on  k,  r,  and  r '  where  their 
presence  is  clearly  implied. 

For  convenience  we  define  the  auxiliary  scalar  functions 


GDik-,r,r')  =  G(A:;r,r')  -  G{k;r,ri)  (9a) 

GN{k\r,r')  =  G{k;r,r')  +  G(A:;r,r/)  (9b) 

It  is  easy  to  verify  that  Gd  and  Gn  satisfy  the  homogeneous  boundary  conditions 

(?£)(/!;?,?')  =  0  forz  =  0  (10a) 

Q 

— Giv(A;;r,rO  =  0  forz  =  0  (10b) 

oz 


We  shall  refer  to  as  the  Dirichlet  scdar  half-space  Green’s  function,  and  Gat  as  the  Neumann 
scalar  half-space  Green’s  function.  We  recognize  G^  as  the  Green’s  function  for  a  horizontal 
electric  dipole  over  a  ground  plane,  and  Gat  as  the  Green’s  function  for  a  vertical  electric  dipole 
over  a  ground  plane  [62,  page  404].  Next,  we  introduce  the  dyadic  function 

r(A:;r,  r')  =  — jA:VG(/!:;r,  ?')  X  I  (11) 

where  I  =  xx  -t-  yy  +  zz  is  the  idemfactor.  One  can  show  that  F  satisfies  the  distributional 
differential  equation  (Appendix  A) 


V  X  V  X  r  -  =  -jkVd{r-r')  x  I 


(12) 
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The  curl  of  F  can  be  expressed  in  the  alternate  forms 

V  X  r(A:;  r,  r ')  -  -jk  [k‘^G{k-,  r,  r ')  +  VVG(fc;  r,  r ')] 

=  -jk  I  [3  +  3jkR  +  (jkRf^  RR 

-  [1  +  jkR  +  {jkRf]  !}  (13) 

where  R  =  (r  —  r ')  /R.  Similarly,  we  define  the  auxiliary  dyadic  functions 

Fi  =  —jk  (vGn  X  It  +  VGd  X  (14a) 

Ta  =  -jk  (vGd  X  It  +  VGn  x  zz^  (14b) 

where  It  =  xx  +  yy  is  the  transverse  idemfactor.  It  is  easy  to  show  that  the  auxiliary  dyadic 
functions  satisfy  the  homogeneous  boundary  conditions 

z  X  Fi  =  0  for  «  =  0  (15a) 

z  X  V  X  F2  =  0  for  z  =  0  (15b) 

as  well  as  the  equations 

z  •  [^A(f)  X  V  X  Fi(A:;  r,  f ')j  =  2z  •  [^A(r)  x  V  x  F(A:;  r,  r ')j  for  z  =  0  (16a) 

z  •  j^A(r)  X  F2(A:;  r,  r ')j  =  2z  •  j^A(r)  x  F(A;;  f,  r ')j  for  z  =  0  (16b) 

where  A(f)  is  any  vector  function  of  position. 


3.4  Radiation  Integrals 

The  scalar  and  dyadic  functions  defined  in  the  previous  section  provide  the  building  blocks 
for  three  theorems  which  are  fundamental  to  the  development  of  the  integral  equations  of  our 
problem.  All  of  these  theorems  are  presented  in  [2]  and  proved  in  Appendix  A.  We  state  them 
here. 
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Theorem  1  Let  V  be  a  homogeneous  region  with  regular  boundary  dV,  n  be  the  outward  unit 
normal  vector  on  dV,  and  F  =  r(A:;r,r')  be  defined  by  equation  (11).  If  A  =  A(r)  satisfies 
V  X  V  X  A  —  k'^A  =  0  V  r  G  V,  then 


n 


|a  X  V  X  r]  +  [v  X  a]  X  r} 


ds  = 


I jkV  X  A(r') 


10 


for  r'  eV 
forr' 


(17) 


where  V  is  the  closure  ofV;  that  is,  the  union  ofV  and  its  boundary. 


Theorem  2  Let  V,  V,  dV,  n,  and  A  be  defined  as  in  Theorem  1.  Let  T\  be  defined  by  equa¬ 
tion  (14a),  Vi  be  the  image  ofV  with  respect  to  the  xy-plane,  and  Vi  be  the  closure  of  Vi.  Then 


j  n  •  I A  X  I^V  X  Fij  +  [^V  X  Aj  x  Fij  ds=  < 


jkV' X  A{r')  forr' eV 

:;fc[v^xA(ri')].  forr' eV 
0  forT'^V\JVi 


(18) 


Theorem  3  Let  V,  V,  dV,  Vi,  V i,  h,  and  A  be  defined  as  in  Theorem  2.  Let  r2  be  defined  by 
equation  (14b).  Then 


J  n  •  |a  X  [v  X  Fa]  +  [v  X  a]  X  Fa}  rfs  =  | 


jkV  X  A(r') 

-jk  [V'  X  A(r/)]^ 


0 


for  r’  eV 
forr'EVi 
forr'  ^VKJVi 
(19) 


The  operator  V  differentiates  with  respect  to  unprimed  coordinates,  while  V'  differentiates 
with  respect  to  primed  coordinates.  The  differential  operator  Vj  is  the  same  as  V'  except  that  ^ 
is  replaced  by  Theorems  1-3  do  not  apply  when  f '  lies  exactly  on  dV,  but  are  valid  as  r ' 
approaches  dV  from  either  the  interior  or  exterior  of  V. 

Much  like  the  source-free  form  of  the  Stratton-Chu  equations  [59],  Theorems  1-3  relate  a 
vector  function  evaluated  at  a  point  to  an  integral  over  a  closed  surface.  We  may  think  of  A  and 
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V  X  A  on  the  surface  as  effective  sources  giving  rise  to  V  x  A  in  the  volume  interior,  and,  in  the 
case  of  Theorems  2  and  3,  V  x  A  in  the  image  of  the  volume  interior.  However,  Theorems  1-3 
are  stated  in  general  terms;  the  only  requirements  are  that  the  vector  function  A(r)  satisfy  the 
homogeneous  wave  equation  (2)  and  that  the  outward  unit  normal  n  be  defined  almost  everywhere 
on  dV.  In  this  section,  we  will  specialize  Theorems  1-3  to  the  OCRED  scattering  problem  shown 
in  Figure  3  by  associating  A(r)  with  various  OCRED  field  components  and  associating  V  with 
either  the  cavity  volume  D  or  the  upper  half  space,  z  >  0. 

A  notational  distinction  arises  when  the  cavity  volume  D  is  filled  with  a  material  whose  con¬ 
stitutive  parameters  differ  from  those  of  free  space.  When  we  associate  V  with  D,  the  wavenumber 
for  the  cavity  medium  appears  in  the  wave  equation  and  the  scalar  and  dyadic  Green’s  functions; 
that  is,  k  =  ki  —  oj^Z/xieT.  When  we  associate  V  with  the  upper  half  space,  the  wavenumber  is 
the  free  space  wavenumber,  k  =  ko  = 

An  exhaustive  list  of  the  possible  equations  arising  from  the  application  of  Theorems  1-3 
to  the  OCRED  problem  appears  in  Appendix  B.  Here,  we  will  only  present  those  results  which 
directly  contribute  to  the  final  integral  equations  of  Section  3.5. 

Application  1.  Let  V  =  D,  A  =  H,  and  r '  €  £>  in  Theorem  1.  Making  use  of  equations  (lb) 
and  (3a),  we  obtain 

kjYiEir')  =  -  j  [n  X  H(f)]  •  V  x  f{ki-,  f,  f ')  ds 

-  jkiYi  j  [zx  E(r)]  •  r(A:i;  r,  f ')  da  for  f'eD  (20) 

where  T i  =  \/ei/pi  is  the  intrinsic  admittance  of  the  cavity  material.  Equation  (20)  shows  that 
the  total  electric  field  at  any  point  in  the  cavity  volume  is  a  function  of  the  tangential  electric  and 
magnetic  fields  on  the  cavity  surface  and  aperture. 
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Application  2.  Let  V  =  L>,  A  =  E,  and  r '  G  £>  in  Theorem  1.  Making  use  of  equations  (la) 
and  (3a),  we  obtain 

klZ{R{r')  =  y  [z  X  E(r)]  •  V  x  r(A:i;r,r')  ds 

-  jkiZi  j  [n  X  H(r)]  •  f{ki;  r,  r ')  da  for  r'  e  D  (21) 

where  Zi  =  y/p^Je[  is  the  intrinsic  impedance  of  the  cavity  material.  Equation  (21)  shows  that 
the  total  magnetic  field  at  any  point  in  the  cavity  volume  is  a  function  of  the  tangential  electric  and 
magnetic  fields  on  the  cavity  surface  and  aperture. 

Application  3.  Let  V  be  the  upper  half  space,  A  =  and  >  0  in  Theorem  2.  Making 
use  of  equations  (la),  (3b),  (15a),  and  (16a),  we  obtain 


^0  ^0  xTScat 


(r  0  =  f  zx  E(f)  •  V  X  r{jto;  r,  r ')  da  for  z' >  0 
J  a 


Equation  (22)  shows  that  the  scattered  magnetic  field  at  any  point  in  the  upper  half  space  is  a 
function  of  the  tangential  electric  field  on  the  cavity  aperture. 

Application  4.  Let  V  =  L>,  A  =  E,  and  r '  G  in  Theorem  3.  Making  use  of  equations  (la), 
(3a),  (15b),  £md  (16b),  we  obtain 

jkiH(f')  =  2^  [2  X  H(r)]  •r(ki;r,r')  da 

+  J  ^hx  H(r)j  ■  r2(fci;r, ?')  ds  for?'  G  D  (23) 

Equation  (23)  shows  that  the  total  magnetic  field  at  any  point  in  the  cavity  volume  is  a  function  of 
the  tangential  magnetic  field  on  the  cavity  surface  and  aperture. 

Application  5.  Let  V  be  the  upper  half  space,  A  =  and  z'  >  0  in  Theorem  3.  Making 
use  of  equations  (lb),  (3b),  (15b),  and  (16b),  we  obtain 


•ft  scat 


(?')  =  f  z  X  E(r)  •  r(A:o;r)rO  da  for  z'  >  0 
J  a 
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Equation  (24)  shows  that  the  scattered  electric  field  at  any  point  in  the  upper  half  space  is  a  function 
of  the  tangential  electric  field  on  the  cavity  aperture.  This  equation  serves  as  the  overall  radiation 
integral  for  the  OCRED  scattering  problem.  Once  z  x  E  on  ct  is  known,  equation  (24)  allows  us 
to  calculate  the  OCRED  echo  area  Ag,  (also  known  as  radar  cross  section)  using  the  well-known 
definition  [21,  page  116] 


Ag  =  lim  47r  |r  T 
Ir'Kcx)  '  ' 


j,scat(-/) 


(25) 


where  the  incident  field  is  a  uniform  plane  wave  having  unit  amplitude. 


3.5  Surface  Integral  Equations 

Equations  (20)-(24)  define  the  total  fields  in  the  cavity  volume  and  the  scattered  fields  in 
the  upper  half  space  in  terms  of  the  tangential  field  components  on  the  cavity  surface  and  aperture. 
Thus  it  is  sufficient  to  find  these  tangential  field  components  to  solve  the  entire  OCRED  scattering 
problem.  In  this  section,  we  develop  surface  integral  equations  which  may  be  used  to  find  these 
tangential  field  components. 

We  notice  that  the  left  and  right  sides  of  equations  (20)-(23)  involve  different  types  of  field 
quantities.  Specifically,  we  see  that  the  right  sides  of  equations  (20)-(23)  involve  the  tangential 
components  of  E  and  H  evaluated  on  dD,  but  that  the  left  sides  of  these  equations  involve  E  and 
H  evaluated  at  points  away  from  dD. 

To  make  equations  (20)-(23)  more  useful,  we  will  let  r '  approach  dD,  and  then  “sift  out” 
the  tangential  component  of  the  result.  First,  we  fix  a  point  r  on  dD,  so  that  n  is  the  outward  unit 
normal  vector  at  r,  as  shown  in  Figure  5.  Second,  we  cross  multiply  the  equations  by  n.  Finally, 
we  let  r '  =  r  -f  (5  n  and  evaluate  the  limit  as  5  — >  0.  In  so  doing,  we  make  use  of  the  following 
theorem  ( [2],  [30,  p.  334],  [62,  p.  354],  [48,  p.  205]): 

Theorem  4  Let  V  be  a  volume  with  regular  bourulary  dV,  t  be  a  fixed  point  on  dV,  n  be  the 
outward  unit  normal  vector  at  r,  and  r(A:;  r,  r')  be  defined  by  equation  (11).  ^A(r)  is  continuous 
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Figure  5  A  portion  of  dD  in  the  neighborhood  of  r.  The  point  r '  approaches  r  along  n  such 
that  r'  —  f  =  6  h. 


on  dV,  then 

lim  n  X  /  A(r")  •  T{k-,r",r')  ds" 

f'->r  Jqy 

=  q:^A(f)  +  nx  [  A{r") -  fit, r ", r')  ds"  (26) 

2  Jav 

where  the  upper  sign  is  taken  ifr'  r  from  the  exterior  ofV,  and  the  lower  sign  is  taken  ifr'  ->  r 
from  the  interior  ofV. 

Theorem  4  provides  a  means  to  treat  the  limit  of  the  integrals  involving  the  dyadic  kernel  F  as 
r'  approaches  dD;  the  limit  as  r '  dD  is  different  from  the  evaluation  of  such  integrals  with 
r '  G  dD.  In  contrast,  integrals  involving  the  dyadic  kernel  V  x  F  exhibit  no  such  discontinuity 
and  behave  as 

lim  fix  f  A{f")  -V  xf{k;r",r')  ds"  =^hx  f  A{r") -V  xf{k;r",r')  ds"  (27) 
?'->?  Jqy  Jqy 

where  V,  dV,  r,  n,  F,  and  A  are  defined  as  in  Theorem  4.  Furthermore,  we  note  that  the  total 
tangential  fields  are  continuous  across  a,  as  shown  in  equations  (4a)  and  (4b). 

We  now  apply  the  procedure  outlined  earlier  in  this  section  to  each  of  equations  (20)-<23). 
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Application  to  equation  (20).  Fix  r  G  cr  so  that  n  =  z,  and  set  r '  =  r  —  5  z,  with  5  >  0.  We 
cross  multiply  equation  (20)  by  z  and  take  the  limit  as  5  — 0.  Applying  Theorem  4,  we  obtain 

— ^i^z  X  E(r')  =  z  X  /  [n  X  H(r)l  •  V  x  r(fei;r,  ?')  ds 
2  JdD  *■ 

+  jkiYiz  X  j  j^z  X  E(r)j  •  r(A;i;f, r')  dcT  forr'Gcr  (28) 

It  is  easy  to  show  that  z  x  |j^z  x  E(r)|  •  r(fci;  r,  r ')|  =  0  when  r  6  a  and  r'  G  a.  Thus  the 
second  integral  in  equation  (28)  vanishes  and  we  are  left  with 

— X  E(r ')  =  z  X  f  [n  X  H(r)l  -  V  x  T{ki; r, r ')  ds  for  r'  E  a  (29) 
2  JdD  ^  -I 

Equation  (29)  shows  that  the  tangential  electric  field  on  a  can  be  expressed  as  a  function  of  the 
tangential  magnetic  field  on  dD. 

Application  to  equation  (21).  Fix  r  G  a  so  that  n  =  z,  and  set  r '  =  r  —  5  z,  with  ^  >  0.  We 
cross  multiply  equation  (21)  by  z  and  take  the  limit  as  (5  0.  Applying  Theorem  4,  we  obtain 


—  jkiZizx  J  n  X  H(r)j  •  r(A:i;r,  r')  ds  forr'Gcr  (30) 


Equation  (30)  shows  that  the  tangential  magnetic  field  on  a  can  be  expressed  as  a  function  of  the 
tangential  magnetic  field  on  S  £md  the  tangential  electric  field  on  a. 

Application  to  equation  (22).  Fix  r  G  cr  so  that  n  =  z,  and  set  f '  =  r  +  z,  with  (5  >  0.  We 
cross  multiply  equation  (22)  by  z  and  take  the  limit  as  5  ^  0.  Applying  Theorem  4,  we  obtain 

— X  H®“‘(f ')  =  z  X  f  z  X  E(r)  •  V  x  r(fco;  f,  r ')  da  for  r '  G  cr  (31) 
^  J  rr 
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Noting  that  H®“*  =  H  —  and  that  z  x  =  z  x  H‘“®  on  a,  this  becomes 

^^^z  X  [H(r  0  -  2H‘"®(r ')]  = 

z  X  j  z  X  E(r)  •  V  X  r(A:o;  r,  r ')  da  for  r '  G  cr  (32) 

J  a 

Equation  (32)  shows  that  the  tangential  magnetic  field  on  a  can  be  expressed  as  a  function  of  the 
tangential  electric  field  on  a  and  the  tangential  incident  magnetic  field  on  a. 

Application  to  equation  (23).  Fix  f  e  S  and  set  r '  =  r  —  5  n,  with  J  >  0.  We  cross  multiply 
equation  (23)  by  n  and  take  the  limit  as  <5  0.  Applying  Theorem  4,  we  obtain 


^n'  X  H(r  0  =  2n'  X  J  [z  x  H(r)]  •  r(fci;  r,  r ')  da 

+  n'  X  J  [n  X  H(f)]  - r2(fci; r, r ')  ds  for  r'  €  S  (33) 

Equation  (33)  shows  that  the  tangential  magnetic  field  on  S  satisfies  a  second-kind  integral  equation 
because  the  unknown  function  n  x  H  appears  both  within  the  integral  and  outside  it. 

Equations  (29),  (30),  (32),  and  (33)  are  coupled  integral  equations  involving  four  functions: 
the  tangential  magnetic  field  on  5,  the  tangential  electric  field  on  a,  the  tangential  magnetic  field 
on  cr,  and  the  incident  tangential  magnetic  field  on  a.  The  last  of  these  functions  is  known  while  the 
first  three  are  not  known  a  priori.  We  associate  these  tangential  field  components  with  equivalent 
current  densities  as  follows: 


3s{r)  =  nx  H(r) 

for  r  G  S 

(34a) 

J<,(r)  =  z  X  H(r) 

for  r  G  a 

(34b) 

M^(r)  =  -z  X  E(r) 

for  r  G  cr 

(34c) 

J°(f)  =  2z  X  H'"®(f) 

for  r  G  cr 

(34d) 

We  call  3s  and  “equivalent  electric  current  densities”  and  the  “equivalent  magnetic  current 
density.”  J*’  is  the  familiar  “physical  optics”  electric  current  density  that  exists  on  an  unbroken 
ground  plane.  With  these  definitions,  equations  (29),  (30),  (32),  and  (33)  become  integral  equations 
involving  three  unknown  current  densities  (3s,  J<t,  and  M^),  along  with  one  known  one  (J°).  We 
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will  find  it  convenient  to  combine  equations  (30)  and  (32)  by  dividing  equation  (30)  by  —jki, 
dividing  equation  (32)  by  —jko,  and  forming  the  difference  of  the  result.  Then,  equations  (29), 
(30),  (32),  and  (33)  become 


z  X 


f  3^{r)-V  xr{ki-,f,r')da  +  zx  f  ^(r)  •  V  x  r(A:i;?,r') 
J(T  JS 


ds 


k^Yi 


Mo^(r')  forr'Gcr 


(35) 


z  X 


Vxr(A:i;r,rO  V  x  r(A:o;r,r') 


-jki 


-jko 


do 


+ 


jkiZi  +  jko  Zo 


3a{r')  -  Z]%  X  ■  r(A:i;r,r')  ds 

j‘^{r')  forr'ea 


jkoZo  ^Of^is 


(36) 


2no  X  f  3tr{r)  ■Tiki]?,?')  da +  no  X  f  3  sir)  ■T2{ki-,T,r')  ds 
Ja  Js 

forr'  e  S 


fun 


(37) 


Equations  (35),  (36),  and  (37)  constitute  a  system  of  three  coupled  integral  equations  involv¬ 
ing  three  unknown  functions:  J^,  J^,  and  M^.  When  this  system  is  solved  and  the  equivalent 
current  densities  J5,  Jcr,  and  are  found,  then  equations  (20)-(24)  can  be  used  to  find  the  fields 
everywhere  in  the  cavity  volume  and  in  the  upper  half  space. 

We  note  that  equation  (35)  may  be  considered  to  be  an  explicit  definition  of  in  terms 
of  3a  and  J5,  which  may  be  used  in  equations  (36)  and  (37)  to  eliminate  as  an  unknown. 
Similarly,  we  note  that  equation  (36)  may  be  considered  to  be  an  explicit  definition  of  3^  in  terms 
of  M(j  and  J5,  which  may  be  used  in  equations  (35)  and  (37)  to  eliminate  as  an  unknown. 
While  these  may  seem  appealing  because  they  reduce  the  number  of  unknown  functions  from 
three  to  two,  this  comes  at  the  price  of  much  more  complicated  integral  equations  involving  iterated 
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surface  integrals.  ^  We  choose  to  keep  the  equations  as  simple  as  possible,  and  thus  we  will  work 
with  all  three  integral  equations  for  the  remainder  of  this  document. 

J.6  The  Degenerate  Case 

We  define  the  degenerate  case  to  be  that  in  which  the  cavity  is  filled  with  free  space.  In  this 
case,  k  =  ki  =  kQ,Y  =  Yi  =10.  and  Z  =  Zi  =  Zq,  and  examination  of  equation  (36)  reveals 
that  the  dyadic  kernel  of  the  integral  involving  vanishes,  so  that  equation  (36)  involves  the 
electric  current  densities  and  J5  alone.  In  other  words,  the  electric  current  densities  and  Js 
effectively  decouple  from  the  magnetic  current  density  M(j.  With  eliminated,  equations  (36) 
and  (37)  form  a  system  of  equations  independent  of  equation  (35).  We  note  that  this  reduced  set  of 
equations  is  identical  to  those  developed  by  Asvestas  and  Kleinman  [2]. 

In  the  degenerate  case,  it  is  possible  to  express  the  scattered  electric  field  ‘mDj  in  terms  of 
J5.  To  see  this,  we  combine  equations  (178)  and  (180)  to  get 

gscat(jf/)  f  -  p/)  for  ?'  G  Df  (38) 

3^  Js 

Then  by  equation  (lb), 

E*“‘(?0  =  ^Vx^Js(f)-ri(A:;f,f')ds  forr'eDf  (39) 

Equations  (36),  (37),  and  (39)  are  all  that  are  needed  to  solve  the  degenerate  scattering  OCRED 
problem. 

3. 7  Summary 

In  this  chapter,  we  have  defined  the  OCRED  geometry  along  with  the  characteristics  of  the 
electromagnetic  fields  for  the  OCRED  scattering  problem.  We  have  introduced  dyadic  functions 
which  form  the  kernels  of  several  integral  equations  defining  the  electromagnetic  fields  everywhere 
in  terms  of  the  equivalent  electric  and  magnetic  current  densities  on  the  cavity  surface  and  aperture. 

'Such  a  reduction  may  be  attractive  when  dealing  with  very  large  geometries  where  matrix  solve  time  dominates  the 
overall  solution  time,  but  in  this  effort  we  will  limit  ourselves  to  smaller  OCREDs. 
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We  have  transformed  these  integral  equations  into  a  set  of  three  coupled  surface  integral  equations 
that  are  satisfied  by  the  equivalent  current  densities.  We  have  shown  that  this  set  is  a  generalization 
of  the  set  developed  by  Asvestas  and  Kleinman. 
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IV.  Numerical  Solution  of  the  Integral  Equations 

The  set  of  three  coupled  integral  equations  (35)-(37)  cannot  be  solved  analytically  for  arbitrary 
cavity  shapes.  However,  we  can  approximate  the  solution  to  an  acceptable  level  of  accuracy  by 
using  the  method  of  moments  [22].  The  method  of  moments  involves  approximating  the  integral 
equations  by  a  system  of  linear  equations,  and  then  solving  the  system  to  find  the  unknown  func¬ 
tions.  Provided  it  is  implemented  judiciously,  the  approximate  solution  produced  by  the  method  of 
moments  approaches  the  exact  solution  if  the  latter  is  available. 

To  start  with  the  method  of  moments  involves  approximating  the  unknown  function  as  a 
linear  combination  of  a  known  set  of  functions,  called  basis  functions.  Usually,  the  basis  functions 
are  of  the  subdomain  type,  which  means  they  are  non-zero  over  only  a  small  portion  of  the  domain 
of  the  unknown  function.  This  allows  great  flexibility  in  approximating  the  unknown  function, 
but  at  the  expense  of  needing  to  find  a  large  number  of  coefficients.  One  useful  choice  of  basis 
functions  [54]  uses  triangular  subdomains  to  model  the  unknown  function  over  arbitrary  three- 
dimensional  surfaces.  These  triangles  are  typically  0.1  A  on  a  side  (where  A  is  the  wavelength),  so 
that  one  square  wavelength  of  surface  requires  on  the  order  of  100  basis  functions.  It  is  easy  to  see 
that  the  problem  rapidly  becomes  unmanageable  as  the  domain  of  the  unknown  function  increases 
in  area  expressed  in  square  wavelengths. 

The  remainder  of  this  chapter  is  organized  as  follows.  The  geometry  of  a  special  class  of 
OCRED  geometry  is  defined  in  Section  4.1.  The  basis  functions  to  be  employed  in  the  method  of 
moments  are  introduced  in  Section  4.2,  while  the  testing  functions  and  inner  product  are  discussed 
in  Section  4.3.  The  details  of  calculating  the  elements  of  the  moment  method  matrix  are  given 
in  Section  4.4.  The  system  excitation,  or  forcing  function,  is  given  in  Section  4.5,  and  some 
computational  savings  reaped  by  exploiting  symmetry  properties  are  shown  in  Section  4.6. 

4.1  Body  of  Revolution  OCRED 

Since  the  backscatter  pattern  of  an  effective  OCRED  is  relatively  insensitive  to  small  changes 
in  aspect  angle,  we  will  restrict  our  attention  to  those  geometries  having  a  backscatter  pattern  that 
is  independent  of  azimuth  angle.  Such  OCREDs  are  symmetric  with  respect  to  the  z  axis,  and 
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are  termed  “body  of  revolution  ”  or  BOR,  OCREDs.  We  shall  see  also  that  this  restriction  greatly 
reduces  the  labor  of  modeling  the  OCRED  backscattering  problem. 


Figure  6  A  cone  frustum  generated  by  revolving  a  line  segment  about  the  z  axis. 

Consider  a  simple,  open  curve  that  lies  in  the  xz  plane  and  does  not  intersect  the  z  axis.  If  we 
revolve  this  curve  about  the  z  axis,  we  generate  an  open  surface  that  is  rotationally  symmetric  about 
the  z  axis.  Such  a  surface  is  called  a  surface  of  revolution,  and  the  curve  is  called  the  generating 
arc.  If  the  curve  is  a  straight  line  segment,  then  the  corresponding  surface  of  revolution  is  a  cone 
frustum,  as  shown  in  Figure  6.  If  we  let  the  endpoints,  and  only  the  endpoints,  of  the  curve  lie  on 
the  z  axis,  then  the  corresponding  surface  of  revolution  becomes  a  closed  surface,  and  the  region 
enclosed  by  the  surface  is  called  a  body  of  revolution.  If  the  OCRED  cavity  interior  D  is  a  body  of 
revolution,  we  call  the  OCRED  a  BOR  OCRED. 

The  generating  arc  for  the  BOR  OCRED  is  a  bounded,  piecewise-smooth  plane  curve  which 
begins  at  the  origin  and  ends  on  the  z-axis  somewhere  below  the  origin.  The  BOR  OCRED  gener¬ 
ating  arc  is  shown  in  Figure  7.  It  does  not  intersect  the  2:-axis  except  at  the  endpoints,  and,  without 
loss  of  generality,  we  restrict  the  generating  arc  to  lie  in  the  xar-plane.  We  parameterize  the  gen¬ 
erating  arc  by  its  arc  length  t  such  that  ^  =  0  corresponds  to  its  endpoint  at  the  origin  and  i  =  L, 
where  L  is  the  length  of  the  generating  arc,  corresponds  to  its  endpoint  on  the  negative  z  axis. 

If  r(£)  is  a  point  on  the  generating  arc,  then  we  see  that  r(^)  is  uniquely  determined  by  its 
coordinates  in  the  xz  plane;  that  is,  r(^)  =  xp(^)  -I-  z  z{tj,  as  shown  in  Figure  7.  Since,  for  the 
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Figure  7  The  generating  arc  of  a  body  of  revolution  (BOR)  OCRED.  The  dashed  portion  gener¬ 
ates  the  aperture  a  while  the  solid  portion  generates  the  cavity  surface  S.  The  arc  length  variable  i 
increases  from  0  to  L.  The  position  vector  r(£)  defines  a  point  on  the  generating  arc,  with  normal 
vector  n(£)  and  tangent  vector  £(£).  The  angle  a(£)  defines  the  orientation  of  £(£). 

BOR  OCRED  geometry,  no  part  of  the  generating  arc  lies  above  the  z  axis  or  to  the  left  of  the  x 
axis,  we  see  that  p(£)  >  0  and  z(£)  <  0  for  0  <  £  <  L. 

If  we  now  allow  r(^)  to  rotate  about  the  z  axis  by  an  amount  (j)  in  the  positive  sense,  then  we 
see  that  every  point  on  the  surface  of  revolution  is  specified  by  the  chosen  values  of  £  and  <j)  in  the 
ranges  0  <  £  <  L  and  0  <  ^  <  27r.  For  arbitrary  £  and  (j),  the  corresponding  point  on  the  surface 
of  revolution  is  then  given  by 

r{£,4))  =  xp(^)  cos4>  +  y  p{£)  sin^  -f  zz{£)  (40) 

It  is  obvious  that  p,  (j),  and  z  are  the  coordinates  of  the  circular  cylindrical  coordinate  system. 

We  have  established  that  the  coordinates  £  and  (p  uniquely  determine  a  point  r{£,  (f))  on 
the  surface  of  revolution.  In  addition,  they  also  determine  several  vector  quantities  of  interest. 
Specifically,  we  define  £{£)  to  be  the  unit  vector  in  the  direction  of  increasing  £,  ^{4>)  to  be  the 
unit  vector  in  the  direction  of  increasing  (f>,  and  h{£,  <j>)  to  be  the  unit  vector  normal  to  the  surface, 
pointing  away  from  the  interior  volume.  It  is  easy  to  see  that  £,  and  n  are  mutually  orthogonal. 
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Thus,  I  and  ^  are  tangent  to  the  surface  at  r(^,  <f)).  The  angle  a(£)  describes  the  orientation  of  £ 
and  n  in  the  xz-plane,  such  that  £  —  x  cos  a  —  y  sin  a  and  n  =  x  sin  a  +  y  cos  a.  The  radius  p(£) 
is  the  magnitude  of  the  projection  of  r(^)  onto  the  rry-plane  and  is  given  by  p(£)  —  |r(£)  •  It|. 


The  generating  arc  for  a  BOR  OCRED  can  have  any  amount  of  curvature  so  long  as  it 
generates  a  valid  BOR  OCRED.  To  make  the  analysis  more  tractable,  however,  we  shall  require 
the  generating  arc  to  be  piecewise  linear.  This  is  not  too  great  a  restriction  since  we  expect  an 
arbitrary  curved  generating  arc  may  be  well  approximated  by  a  piecewise  linear  one,  especially  if 
the  line  segments  are  small  compared  to  a  wavelength. 


Consider  a  piecewise  linear  BOR  OCRED  generating  arc  consisting  of  N  segments  joining 
N  +  1  nodes.  We  number  the  nodes  in  ascending  order  such  that  each  node  corresponds  to  a 
progressively  larger  value  of  the  arc  length  variable  £.  Then  we  may  write  the  relation 


0  =  ^0  <  <  ^2  <  •  ■  ■  <  ^N+1  —  L 


(41) 


The  number  of  segments  that  generate  the  aperture  a  is  Na  and  the  number  of  segments  that 
generate  the  cavity  surface  5  is  iST,,  such  that  N  =  Na  +  Ng.  Now  suppose  r(£)  is  a  point  on  the 
n**’  segment;  that  is,  £n-i  <  £  <  £n-  Let  pn  =  p{£n)  and  =  z{£n).  Then  any  point  r(^)  on  the 
n*  segment  is  fully  characterized  by  the  linear  relations 


p{£)^ 

z{£)  = 


Pn—\  {£n  Pn  (£  £n—l) 

£n  £n—\ 

Zn-l  {£n  -£)+Zn{£-  £„-l) 


t. 


n— 1 


for £n-\  <£<£n 

(42a) 

for  £n-l  <£  <£n 

(42b) 

The  orientation  angle  a{£)  is  constant  on  the  n***  segment  and  is  given  by 


a{£)  =  an 

=  arctan  — —  for  <  £  <  £n 

Pn  Pn—l 


(43) 
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Figure  8  A  piecewise-linear  generating  arc.  Na  segments  generate  the  aperture  a  and  Ng  seg¬ 
ments  generate  the  cavity  surface  S.  The  n*  segment  falls  between  r(^„_i)  and  r(^„). 

4.2  Basis  Functions 

The  first  step  in  the  method  of  moments  is  to  define  a  set  of  basis  functions  capable  of 
accurately  approximating  the  unknown  functions.  Referring  to  equations  (35)-(37),  we  see  that  the 
unknown  functions  we  wish  to  approximate  are  and  J<7(r)  for  r  €  a  and  3s{^)  for  r  e  S'. 

The  basis  functions  we  will  use  here  are  the  same  as  those  used  by  Mautz  &  Harrington  [41]  and 
Medgyesi-Mitschang  &  Putnam  [44],  These  basis  functions  have  a  piecewise  smooth  variation 
with  £  and  complex  exponential  (Fourier)  variation  in  <f>.  We  define  the  basis  functions  formally  as 

('Ma) 

(44b) 

where  I  =  x  cosa(^)  coscp  +  y  cosa(^)  sin(^  —  z  sina(£),  ^  =  — x  sin^  -|-  y  coscj),  and 
Tn{£)  is  a  “triangle”  function  as  shown  in  Figure  9  and  given  by  the  equation 
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Tn{i) 


Figure  9  Triangle  function,  T„(^).  T„(^)  is  the  piecewise  linear  function  having  unit  amplitude 
at  in  and  vanishing  at  all  other  nodes  on  the  generating  arc. 


Tn{i)  =  { 


/  for4-i<^<4 

for4<^<4+i 


0 


otherwise 


(45) 


When  n  =  Na,  v/e  will  treat  Tn{i)  to  be  either  the  left  half-triangle  or  the  right  half-triangle, 
depending  on  whether  Tn{i)  is  being  used  in  a  basis  function  with  support  on  the  aperture  or  the 
cavity  surface,  respectively.  This  is  necessary  to  allow  for  the  possible  discontinuity  of  the  electric 
and  magnetic  current  densities  in  the  neighborhood  of  the  aperture  rim.  Using  the  basis  functions 
of  equation  (44),  we  approximate  the  unknown  functions  by 

M  Na 

M,(£,,^)«Zo  E 

m=—M  n=l 
M  Na 

E  E 

m——M  n=l 

M  Na+N,-l 

Js{i,<P)^  E  [cL^Lni^A)  +  CL^in{i,<l>)\  (48) 

m=—M  n=Na 

where  M  is  the  highest-order  Fourier  mode  to  be  included  and  A^tn,  and  Cmt  are  complex 
coefficients  to  be  determined.  The  choice  of  M  will  be  discussed  in  Section  5.4.3,  but  it  should  be 
proportional  to  the  electrical  size  of  the  aperture,  that  is  M  oc  iNal^-  Likewise,  Na  and  Ng  will 


bL  <!>)  +  Bin  ^Li^,  4>)\  (47) 


'<n  </>)  +  Ain  ^)]  (46) 
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increase  linearly  with  frequency.  In  equations  (46)  and  (47),  0  <  £  <  and  0  <  ^  <  27r,  while 
in  equation  (48),  £jVa  <  ^  <  ^Na+N,  and  0  <  <^  <  27r.  To  enhance  numerical  stability,  we  scale  the 
expansion  of  by  Zq  so  that  all  expansion  coefficients  will  be  of  similar  orders  of  magnitude. 

The  behavior  of  the  fields  near  the  rim  of  the  aperture  provides  valuable  information.  As¬ 
suming  that  the  radius  of  the  aperture  rim  is  not  too  small,  the  electric  and  magnetic  fields  near  the 
aperture  rim  obey  the  same  conditions  as  the  fields  near  the  tip  of  a  two-dimensional  PEC  wedge, 
as  shown  in  Figure  10.  In  the  two-dimensional  wedge  geometry,  g  measures  the  distance  from 
the  wedge  tip  and  (p  measures  the  angle  from  one  of  the  wedge  faces.  This  geometry  is  studied 


Figure  10  Two-dimensional  wedge  geometry  corresponding  to  the  aperture  rim.  g  measures 
distance  from  the  wedge  tip.  ip  measures  angle  counterclockwise  from  a'^.  The  region  0  <  p  <  tt 
is  the  upper  half  space;  the  surface  y?  =  tt  is  the  aperture  cr;  the  region  tt  <  (/?  <  27r  —  a  is  the 
cavity  volume  D;  and  the  surface  y?  =  27r  —  a  is  the  cavity  surface  S. 


by  Meixner  [45],  who  shows  that  the  total  magnetic  field  parallel  to  the  wedge  tip  is  finite  and 
invariant  with  as  ^  >  0.  Then  by  the  relations  J^r  =  z  x  H  and  J5  =  n  x  H,  we  see  that 


lim  Jff  •  £  =  lim  3s  •  £;  that  is,  the  £  component  of  the  electric  current  density  on  5  U  a  is 
continuous  and  finite  in  the  neighborhood  of  the  aperture  rim.  In  contrast,  the  magnetic  field  per¬ 


pendicular  to  the  wedge  tip  is  unbounded  and  varies  with  p  as  g  ^  0,  so  that  lim 


Jcr  •  ^ 


and 


lim 


3s-4> 


are  unbounded  and,  in  general,  not  equal.  Turning  our  attention  to  the  electric  field, 


we  find  that  the  electric  field  parallel  to  the  wedge  tip  vanishes  as  p  — )•  0,  while  the  electric  field 


perpendicular  to  the  wedge  tip  is  unbounded  and  varies  with  as  p  — 0.  Then  using  the  relation 


M(r  =  — z  X  E,  we  see  that  lim 


M<,-£ 


vanishes,  while  lim 

£-^£Na 


is  unbounded.  This  a 


33 


priori  knowledge  of  the  behavior  of  J(j,  Jg,  and  in  the  neighborhood  of  the  aperture  rim  will 
be  extremely  valuable  in  later  calculations  (specifically,  Sections  4.4.2  and  4.4.3). 


4.3  Testing  Functions  and  Inner  Product 

The  second  step  in  the  method  of  moments  is  to  define  a  set  of  testing  functions  and  an 
appropriate  inner  product.  In  accordance  with  Galerkin’s  procedure,  we  will  use  testing  functions 
7  that  are  of  the  same  form  as  the  basis  functions. 


Tn'in  P 

p[e) 

Tn'in  ,, 

im'*' 


(49a) 

(49b) 


where  £'  =  H  cos  oi{£')  cos  (j)'  +  y  cos  Oi{£')  sin  (j)'  —  z  sin  a{i')  and  =  — x  sin  +  y  cos  (f)'.  As 

an  inner  product,  we  will  choose  the  following: 


(A,B)=  f  A(f')  •B*(r')  ds'  (50) 

Jsua 

where  B*  denotes  the  complex  conjugate  of  B  and  5  U  cr  =  dD  is  the  boundary  of  the  cavity 
interior.  This  choice  of  testing  function  and  inner  product  allows  us  to  realize  the  major  advantage 
of  the  BOR  symmetry  [41].  Since  =  0  for  m  7^  m',  we  see  that  the  basis  and 

testing  functions  are  orthogonal  with  respect  to  the  inner  product  when  their  Fourier  indices  are 
different.  In  other  words,  the  different  Fourier  modes  do  not  influence  each  other  and  effectively 
decouple.  This  will  allow  us  to  consider  each  Fourier  mode  separately. 

Applying  our  choices  of  basis  function,  testing  function,  and  inner  product  to  the  coupled  in¬ 
tegral  equations  (35)-(37)  results  in  a  system  of  linear  equations  which  can  be  expressed  compactly 
in  matrix  form  as 


A 

0 

Zo{£s^f,j)  ifci^i±«(^,7)  -Z^{Cd,1) 

B 

= 

^(J0,7) 

0  2  (£5^,  7)  (-^6^,7) 

C 

0 
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where 


X 

II 

f  Vxr(A:i;r,r')  , 

1^-  d. 

for  r '  G  cr 

(52) 

C2tf  =  Z  X 

f  V  X  T{ki-,f,f')  J 

//M-  -A  * 

for  r '  G  (7 

(53) 

Cz'if  =  z  X 

f  *-)  ^  ]  dc 

la  [  -3^1  -3ko  j 

for  r'  G  a 

(54) 

-3 

X 

II 

f  ip{r)  -Tiki^r,?')  ds 

Is 

for r'  E  a 

(55) 

Cs'ip  =  n'  X 

j  'if{Y)  ■T2{ki-,r,r')  da 

J  a 

for  t'eS 

(56) 

Ce-tj}  =  n'  X 

J^tp{r)  ■f2{ki]f,f')  ds  -  ^i?(r') 

forr'  G  S 

(57) 

Equation  (51)  is  often  written  in  compact  form  as  Z I  =  V,  where  Z  is  called  the  generalized 
impedance  matrix,  7  is  a  column  vector  containing  the  desired  expansion  coefficients,  and  V  is  a 
column  vector  containing  the  system  excitation.  The  dimensionality  of  Z  is  2(2M  +  l)(2iVa  + 
Ns),  which  is  equal  to  the  total  number  of  unknown  coefficients.  However,  by  exploiting  the 
Fourier  orthogonality,  equation  (51)  is  equivalent  to  2M  +  1  smaller  systems  {Zmlm  = 
m  =  0,  ±1 ...  ±  M),  each  corresponding  to  a  different  Fourier  index  m  and  having  dimension 
2{2Na  +  Ns).  This  decomposition  reduces  the  numerical  complexity  of  the  problem  in  two  ways. 
First,  the  number  of  matrix  elements  that  must  be  computed  is  reduced  by  a  factor  of  2M  + 
1.  Second,  under  the  assumption  that  the  computer  time  required  to  solve  a  system  of  N  linear 
equations  using  Gaussian  elimination  is  0{N^)  [6],  solving  the  (2M  +  1)  smaller  systems  can  be 
accomplished  (2M + 1)^  times  faster  than  the  equivalent  larger  system.  Thus,  exploiting  symmetry 
allows  us  to  solve  the  BOR  OCRED  scattering  problem  (2M  + 1)^  times  faster  than  solving  a  non- 
BOR  OCRED  scattering  problem  of  equivalent  size. 

4.4  Calculating  Elements  of  Zm 

Each  element  of  Zm  corresponds  to  the  coupling  between  a  basis  function  and  a  testing 
function.  Typically,  this  involves  calculating  an  iterated  surface  integral  where  the  dyadic  kernel 
is  V  X  T{k;  f,  r '),  T{k-,  r,  r '),  or  r2(A;;  r,  r').  It  is  convenient  to  partition  Zm  into  9  block  sub- 
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matrices; 


Zm.  — 


yll  7I2  7I3 

•^m 

721  722  723 

‘^m 

731  732  733 


(58) 


where  Z]l  denote  the  elements  of  Z^  that  come  from  testing  equation  (35),  (36),  or  (37)  (i  =  1, 2, 3 
respectively)  and  expanding  M<^,  3  a,  or  (j  =  1, 2, 3  respectively).  Examination  of  equation  (51) 
reveals  that  the  elements  comprising  vanish  identically.  The  other  eight  submatrices  are  treated 
individually  in  the  following  sections. 


Figure  11  Domains  of  •if{r)  and  7(r')  for  Z^,  Z^,  and  Z^.  Both  r  and  r'  are  defined 
on  annuli  in  the  xy-plane  formed  by  revolving  a  portion  of  the  aperture  generating  arc  about  the 
z-axis. 

4.4.1  Z^.  Z}^  is  the  2Na  x  2Na  submatrix  associated  with  the  testing  of  equation  (35) 

and  the  expansion  of  M(t  .  As  shown  in  Figure  1 1 ,  the  domain  of  both  the  basis  and  testing  functions 
is  confined  to  the  aperture  o.  The  elements  of  Z}^  are  found  from  the  inner  product 


(59) 
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where  m  =  0,  ±1,  ±2 . . .  ,  ±M,  1  <  n  <  iVa,  and  1  <  n'  <  Na,  Since  I'  J.  the  “cross  terms” 
Tin')  ^nd  Vanish.  We  are  left  with 


(V’mn )  Tmn' )  (V^mn’Tfnn') 

r-27r  piNa 


Tnji') 

Pin 


nl 


P(^0  J 


/9(^')  di'  dcj)' 


(60) 


We  note  that  pii')  =  (!  for  0  <  ^'  <  ^jv„  so  that  the  integrand  of  equation  (60)  is  a  rational 
function  of  polynomials  in  and  the  integral  can  be  solved  in  closed  form.  We  look  at  three  cases 
depending  on  the  values  of  n  and  n'.  If  n  and  n'  differ  by  two  or  more,  the  last  integral  vanishes 
because  the  product  of  the  two  triangle  functions  is  zero  when  their  respective  supports  do  not 
overlap.  If  |n'  —  nj  =  1,  then  we  have 


i'^Ln^fmn')  =  (^L-oin') 

-  27r  n 

(in'  -  ^nf'  Jin 


jin'  -iW- in) 


(4' -4)^  2  +'”'”''"4' 


for  In'  —  n|  =  1 


(61) 


If  v!  =  n,  then  we  have 


^  (^'-4-i)" 


=  27r 


lin-i  ^'(4-4 


=  27r 


f  4-34-1 
1 2  (4  -  4-1 ) 


+ 


4- 


n—l 


+ 


I 


in  4—1 
2 


In 


^'(4+1-4)^ 

in 


(62) 


f 

Ctj  _ 


n—l 


n+1 


in 

'2(4+1 -4)  '  L4+1-4 


■"¥} 


for  n'  =  n 


In  equation  (62),  the  following  special  cases  must  be  examined.  First,  if  n  =  1,  then  we  must 
evaluate  equation  (62)  in  the  limit  as  4  ^  0  to  avoid  taking  the  logarithm  of  infinity.  Second,  if 
n  =  1,  then  we  ignore  the  terms  involving  4+i  due  to  the  half-triangle  function  used  in  the  basis 
and  testing  functions  at  the  aperture  rim. 
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4.4.2  Z^.  is  the  2Na  x  2Na  submatrix  associated  with  the  testing  of  equa¬ 

tion  (35)  and  the  expansion  of  J^.  As  shown  in  Figure  11,  the  domain  of  both  the  basis  and 
testing  functions  is  confined  to  the  aperture  a.  From  equation  (13)  we  have  V  x  r(A:i;  r,  r ')  = 
—jki  (kfi  +  Vvl  G(fci;  r,  r '),  the  elements  of  Z^  can  be  found  from  the  inner  product 


■  [*  X  *t«]  }  ?')  dc  d^' 

+  (  [®(f)  VVG()!i;f,r')]}  daA/  (63) 


where  m  ==  0,  ±1,  ±2 . . .  ±M,1  <  n  <  Na,  and  1  <  n'  <  Na.  The  first  integral  in  equation  (63) 
is  a  convergent  integral,  though  improper  due  to  the  R~^  singularity  of  the  scalar  Green’s  function 
at  r  =  r'.  The  second  integral,  however,  does  not  converge  in  the  classical  sense  due  to  the 
singularity  of  VVG  at  f  =  r ',  and  must  be  interpreted  as  a  two  dimensional  Cauchy  principal 
value  (CPV)  integral  [46].  However,  this  interpretation  is  troublesome  when  r  lies  on  the  boundary 
of  a  because  the  CPV  concept  requires  the  division  of  a  into  two  parts,  one  of  which  is  a  disk 
centered  at  r  with  radius  ^  >  0.  To  give  meaning  to  the  second  integral,  we  regularize  it  by  using 
the  following  theorem  (proved  in  Appendix  C). 


Theorem  5  Let  A(r)  and  B(r)  be  continuous,  piecewise  differentiable  vector  functions  defined 
on  the  regular  surfaces  S\  and  S2,  respectively,  such  that  n'  •  A  =  0  on  and  h-'B  =  Q  on  S2 
where  n'  is  the  unit  normal  on  Si  and  n  is  the  unit  normal  on  82.  Then 


J  j  A(r')  •  n'  X  [B(r)  ■  VVG(/i:;  r,  r ')]  ds  ds'  = 

[  [  G{k-  f,  f ')  k  •  B(r)]  { V;  ■  [n'  x  A(r ')]  }  ds  ds' 
J  Si  J  S2 

+  [  G{k;r,r')  [Vj  •  B(r)l  A(r')  •  dc '  ds 

JS2  JdSi  ^  ^ 


(64) 


+  f  /  G{k-,r,r')  |  V' •  |^n' x  A(r')j  |  |^B(r)  x  n  ■  dc  ds' 

J  Si  J  dS2 

+  (f  (f  G{k;f,r')  [B(r)  x  nl  •  dc  A(r')  ■  dc ' 

Jasi  JdS2 

where  the  line  integrals  are  traversed  in  the  positive  sense  with  respect  to  the  normal  vectors. 
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Theorem  5  is  valid  for  any  regular  surfaces  Si  and  S2,  even  those  that  overlap.  In  particular, 
Theorem  5  is  valid  when  Si  =  S2-  If  we  associate  both  and  S2  with  cr,  and  let  A(r')  = 
^^*(r ')  and  B(r)  =  then  equation  (63)  becomes 

II  {ttnn  ■  [*  X  }  G{ki-,r,  r  0  da  da' 

J  j  (vi  •  [z  X  ott*(r')]  }  G{k-r,r')  da  da' 

/  fa  G{k\T,r')  ■  dc'  da 


+ 


(65) 


+ 


[  i  { Vi  •  [z  X  (r ')]  }  G{k-,  r,  r ')  x  z]  •  dc  da' 

J  o  J  do^ 

+  <f  (f  G{k] r, f ')  tet(r)  x  zl  •  dc  •  dc' 

Jda  Jda  ^ 

where  both  dc  and  dc'  are  differential  vectors  along  the  aperture  rim,  pointing  in  the  direction  of 
increasing  (f).  We  note  that,  with  Si=  a  and  S2  =  a. 


dc  =  ^  d(l) 
da  =  £  di  d(t) 


dc'  =  iNa  d(l>' 
da'  =  £'  dl'  d(t>' 


(66) 

(67) 


we  can  write  equation  (65)  as 


27r  ri^a  /■27r  f^Na 


r2n  rlNa  /•27r  ptNa  r  1 

+  j  J  j  J  [v,  •  V'(^,  <^)J  {Vi  -[ixf  {£',  cl>')]  ]Gtd£d(l>t'  d£'  d<t>' 

/>27r  /-fATa  /-Stt  P 

+  ^iv„  I  {Vi-  [ZX7*(^',,^')]}  G  [V’(^iv„,<^)xzJ  ■<j>d<j>t'd£'d(j)' 

+  £Naj^  G  \^{l,  0)  X  z]  •  <^  d(/>  TilNaA')  ■  d<i)' 


(68) 


where  'ip  =  'ip^,  7*  =  and  G  =  G{ki;T,  f ')  have  been  introduced  to  improve  the  readabil¬ 
ity  of  the  formula. 
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It  is  easy  to  see  that  all  the  integrands  in  equation  (68)  are  functions  of  I,  H! ,  and  (j)  -  <j)'. 
Employing  the  change  of  variables  9  =  (t)  —  (j)'  and  after  some  tedious  algebraic  manipulations,  we 
can  write  equation  (68)  as 


f-Na  2TriNa 


k 


III  {ni\0)-[ixi;ie,9)]}  ididoe'di' 

0  0  0 
t-Nai'KiNa 

III  {V,-[zxrie',0)]}  eduee'di' 


0  0  0 

iNa^TT 

+  iNa  1 1  [Vs-'ip{i,e)]  r{iN.,o)-yedede 


0  0 


C  C  r  T 

+  j  J  {v;-[zxr(^',o)]}  [p{iN^,e)xz\4edede' 


0  0 
27r 


/p-]k\R  p  _  n  ^ 


•  y  d9 


(69) 


where  R  =  i?(^,  i',  9)  =  y(7-r£^)^  +  2^€'(l  — cos0),  i  =  l^a  in  the  fourth  and  fifth  integrals, 
and  i'  =  ^Na  in  the  third  and  fifth  integrcils.  It  will  be  shown  in  Section  4.4.3  that  the  last  two 
integrals  in  equation  (69)  exactly  cancel  with  two  other  integrals  associated  with  and  thus  do 
not  contribute  to  Z^. 

The  weak  singularity  of  the  first  two  integrals  in  equation  (69)  can  be  extracted  and  integrated 
analytically  using  the  relation  [19] 


1 

2 


TT  b 


^ded9  = 
R 


0  a 


6£(()-«£(f) 
e'E{^)-e'E{p)  -  ^K{^)  + 


£'<  a  <  b 
a  <  i'  <b 
a  <b  <  i' 


(70) 


where  K  and  E  are  the  complete  elliptic  integrals  of  the  first  and  second  kinds,  respectively. 
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Figure  12  Domains  of  ^(r)  and  y(f')  for  and  Z^.  r  is  defined  on  an  annulus  in  the  xy- 
plane  formed  by  revolving  a  portion  of  the  aperture  generating  arc  about  the  2:-axis.  r '  is  defined 
on  a  cone  frustum  formed  by  revolving  a  portion  of  the  cavity  surface  generating  arc  about  the 
z-axis. 


4.4.3  Z^.  Z^  is  the  2Na  x  2Ns  submatrix  associated  with  the  testing  of  equation  (35) 
and  the  expansion  of  J5.  As  shown  in  Figure  12,  the  domain  of  the  basis  functions  is  confined  to 
the  cavity  surface  S  while  the  domain  of  the  testing  functions  is  confined  to  the  aperture  a.  From 
equation  (13)  we  have  V  x  r(A:i; r, r ')  =  —jki  (kfi  +  VV^  G{ki;r, f '),  the  elements  of  Z}^ 
can  be  found  from  the  inner  product 


+  dB<h'  (71) 

where  m  =  0,  ±1,  ±2 ...  ±  M,  Na  <  n  <  Na  +  Ng  —  1,  and  1  <  n'  <  Na-  Employing  the  same 
regularization  procedure  as  used  to  calculate  the  elements  of  Zj^,  we  obtain 
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(72) 


klJJ^  ■  [2  X  }  G{ku  r,  r ')  da  da' 

+  y*  [Vs  •  {v;  •  [2  X  7^t*(r')]  }  G(A:i;F,F')  da  da' 

+  jsfd  [^s"^mn(r)]  G{ki-T,f')if;fj{r')-dc' ds 

/  fas  ^  ^mt  (?)  X  z]  •  dc  c?a' 

<f  (f  G{ki]f, r')  tet(?)  X  zl  •  dc  ■  dc' 

Jda  JaS  ^ 


+ 

+ 


where  dc  is  a  differential  vector  along  the  aperture  rim,  pointing  in  the  direction  of  decreasing  4> 
and  dc '  is  a  differential  vector  along  the  aperture  rim,  pointing  in  the  direction  of  increasing  cj). 
We  note  that,  with  Si  =  a  and  S2  =  a. 


dc  =  —0  d(f) 
da  =  £  d£  d<j) 


dc'  =  £Na  d4>'  (73) 

da'  =  £'  d£'  dcj)'  (74) 


we  can  write  equation  (72)  as 


2n^Na  2ir  L 

kiJJfjG  [T{£',4>')-[zxxf{£,<f>)^^  £d£d(t)£' d£' d<f)' 

0  0 

2Tr  ^Na  27r  L 

+  /  j  j  J  G  [Vs-  i^{£,  0)]  { v;  •  [z  X  7*  {£',  (f)')]  }£d£d(l>  £'  d£'  d(f>' 

0  0  0£jva 
2tt  L  2w 

+  ^Na  1 1  Jo  [V,  ■rP{£,cf>)\  i*{£N^,(f>')4'dcl>'£d£dcf> 


0£iv„  0 
2n^Na  27r 


-iNa  I  j  I G  {v;  •  [z  X  r(^',  <f>')] }  [rfi^N^A)  X  z]  ■  0  d0  d£'  dcj^' 

000 

25r  27r 

-  I  Jg  [^{£N.,ct>)  X  z]  ■  ^  d<^  r  (^iV„,  0')  •  0'  dcl>' 


(75) 
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where  xf  =  7*  =  ^  ^0  have  been  introduced  to  improve  the  readabil¬ 

ity  of  the  formula. 

Proceeding  as  in  section  4.4.2,  we  make  the  change  of  variables  0  —  (f)  —  (j)'  and  write 
equation  (75)  as 


Jrarv 


kl 


-I- 


III 

0 

iNa2n  L 

in 

0  Of^Va 
L  27r 

ni 

^Na  ^ 
iNa^TT 


,6) 

2R{£J',e) 

^-jkiR(e,£',e) 

2R{£,1',6) 

]kiR{l,tNafi) 


2R{i,iN^,e) 

g— jfel  ^(!  ,0) 


|7*(£',  0)  •  [z  X  ^{t,  0)]  £dlde  £'  d£' 

[v,  •  •0(£,  0)]  { v;  -[2X7*  {£\  0)]  ]  £dtde  (!  di’ 


f{iN^,o)-y£dede 


(76) 


2R(e,  ,e,o)  W-[sxf«'.o)]}  [«£«.,«)  x  *] 

0  0  “ 

27r 

r  )^iVa  r  _*  1 

2Ji(£;,.,£,v.,9)  ['K'lv-S)  x2]  •';‘r«A-..0)-yd9 


where  i?(£,^',0)  =  ^[p{£)  -  4-  2p{£)£'  (1  -  cos0)  +  z{£)^. 

Because  of  the  continuity  of  J  -  ^  in  the  neighborhood  of  we  see  that  the  last  two  integrals 
of  equation  (69)  and  equation  (76)  exactly  cancel,  and  so  never  need  be  calculated.  In  fact,  the 
presence  of  these  integrals  in  the  first  place  is  an  artifact  of  the  separation  of  J  into  and  Js- 
This  separation  is  justified  because  of  the  discontinuity  of  J  •  (f)  across  the  aperture  rim,  but  J  •  £ 
could  easily  be  modeled  as  a  single  continuous  function  from  0  <  £  <  L. 

4.4.4  Z^.  is  the  2Na  x  2Na  submatrix  associated  with  the  testing  of  equation  (36) 

and  the  expansion  of  M^.  This  case  dictates  that  the  domain  of  both  the  basis  and  testing  functions 
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is  a.  The  elements  of  can  be  found  from  the  inner  product 


da  da'  (77) 

where  m  =  0,  ±1,  ±2 . . .  ±  M ,  Na  <  n  <  Na  +  Ng,  and  1  <  n'  <  Na- 
From  equation  (13),  V  x  r{k-,  f,  T')/{—jk)  can  be  written 


1.1. 


h 


Z  X 


Vxr(A:i;r,rO  V  x  r(A:o; ?, ?') 


-jh 


-jh 


Vxr(fc;f,rO 


RR 


-jk  ^{[^  +  ^jkR  +  ijkRf 


(78) 


where  R  =  (r  —  f ')  /R.  Then  the  dyadic  kernel  of  equation  (77)  can  be  written  as 


V  X  r(fci;r,fO  _  V  X  r(fco;r,rO  _ 

-jh  -jko 

f  [3  +  ZjkiR  +  (jkiR)^]  -  [3  +  ZjkoR  +  {jkoR)^]  \ 

1  47rR3 


I 


RR 


[1  +  jkiR  +  [jkiRf]  -  [l  +  jk^R  +  hk^R?]  = 

47rR3  f 


Expanding  the  complex  exponential  in  a  Taylor  series  about  i?  =  0, 


e-0i^R  =  l-jkR  +  O{R?), 


(80) 


we  find  that 


V  X  r(A:i;r,f')  _  V  x  r(fco;f,r') 

-jki  -jko 

for  arbitrary  vectors  A  and  B.  As  a  consequence  we  see  that,  although  each  of  the  dyadic  functions 
V  X  T{ki;r,r')/{—jki)  and  V  x  r(fco;r,r')/(— yfco)  has  a  R~^  singularity  at  r  =  r',  their 
difference  has  only  a  R~^  singularity  there.  Thus  the  integral  in  equation  (78)  is  merely  weakly 


B 


=  0{R-^), 


R^O  (81) 
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singular  and  can  be  handled  numerically  using  the  extraction  procedure  outlined  at  the  end  of 
Section  4.4.2.  Of  course,  in  the  degenerate  case  of  an  unfilled  OCRED  cavity,  =  ko  and  the 
dyadic  kernel  identically  vanishes. 

4.4.5  Z^.  is  the  2Na  x  2Na  submatrix  associated  with  the  testing  of  equation  (36) 
and  the  expansion  of  3cr.  This  case  dictates  that  the  domain  of  both  the  basis  and  testing  functions 
is  confined  to  the  aperture  a.  The  elements  of  Z^  are  found  from  the  inner  product  *) 

where  m  =  0,  ±1,  ±2 . . .  ±  M,  1  <  n  <  Na,  and  1  <  n'  <  JVq.  This  inner  product  is  given  in 
Section  4.4.1,  and  we  simply  restate  it  here  for  completeness. 


=  (^L>7mn')  =  0 

{'>PLn,lLn>)  = 


(82a) 


=  < 


2t:  I 


in  1 

2  {in  ~  in—\) 


+ 


“n—\ 


L^n  in—X 


In 


+ 


27r 


0 


in  3^n+l 
2  {in+1  ~  in) 

,p  p  in. 

2  ’^n^n'  f  , 


in+X 


in-X 

2 


L^n+1  in 


In 


in+X 


>  n  —  n 

\n'  —  n|  =  1 
In'  —  n|  >  1 


(82b) 


If  n  =  n'  =  1  in  equation  (82b)  then  the  term  involving  In  is  ignored,  while  if  n  =  n'  = 

Na  then  the  terms  involving  in+x  are  ignored,  for  the  same  reasons  as  given  in  Section  4.4.1. 


4.4.6  Z^.  Z^  is  the  2Na  x  2Ns  submatrix  associated  with  the  testing  of  equation  (36) 
and  the  expansion  of  J^.  This  case  dictates  that  the  domain  of  the  basis  functions  is  confined  to  the 
cavity  surface  S  while  the  domain  of  the  testing  functions  is  confined  to  the  aperture  a,  as  shown 
in  Figure  12.  The  elements  of  Z^  are  found  from  the  inner  product 


mn^  ii 


i)  =  /  (r')  •  ^  [^mn(?)  ■  r(A:i;  r,  r ')]  }  ds  da'  (83) 
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where  m  =  0,  ±1,  ±2 . . .  ±  M,  Na  <  n  <  Na  +  Ng  —  1,  and  1  <  n'  <  Na-  Writing  r(fc;  r,  r ')  = 
—jk  VG{k;  r,  f ')  x  I,  we  obtain 


=  -Jfc  ■  [vG(A:;r,r')  x  I 

=  -3k^f^Jni^xVG{k-,r,r') 


where  we  have  used  the  identity  A  •  x  x  •  C.  Explicitly  writing 

(85) 

(86) 

where  R  =  r  —  f equation  (83)  becomes 


p-jkR  ^ 

VG{k-v,r)  =  {-jkR-l)^R 


=  i-jkR-l) 


G{k-r,r'). 


{-jkR  -  1)  ds  da'  (87) 

Although  the  integrand  in  equation  (87)  has  a  0{R~‘^)  singularity  at  r  =  r since  the  cavity 
surface  S  and  aperture  a  do  not  overlap,  the  integral  converges.  We  now  calculate  each  of  the  four 
permutations  of  the  vector  portion  of  the  integrand  of  equation  (87).  We  have 


£  —  X  cos  O' cos  (/)  +  y  cos  a  sin  ^  —  zsina  (88a) 

^  =  —X  sin  </>  +  y  cos  (p  (88b) 

?  =  X  cos  (p'  +  y  sin  (p'  (88c) 

p'  =  -X  sin  (p'  +  y  cos  (p'  (88d) 

R  =  X  (/9 cos (p  —  p' cos 4>)  +  y  (p sin 0  —  p' sin p')  —iz  (88e) 
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from  which  we  find 


i  -  ix 

X  =  [pcos{^  —  4>)  —  p'^  sin  a  +  2rcosacos(</i  —  (j)') 

(89a) 

■  i  X 

X  =  [psina  +  ^cosa]  sin(^  —  (/>') 

(89b) 

a  -  zx  1 

X  R^  =  —z sin(<^  —  4>) 

(89c) 

-  z  X  1 

X  R^  =  z cos{<f)  —  4>) 

(89d) 

Figure  13  Domains  of  ip{r)  and  7(f ')  for  Z^.  r  is  defined  on  a  cone  frustum  formed  by 
revolving  a  portion  of  the  cavity  surface  generating  arc  about  the  z-axis.  r'  is  defined  on  an 
annulus  in  the  a;y-plane  formed  by  revolving  a  portion  of  the  aperture  generating  arc  about  the 
2r-axis. 

4.4.7  Z^.  Z^  is  the  2Ns  x  2Na  submatrix  associated  with  the  testing  of  equation  (37) 
and  the  expansion  of  J^.  As  shown  in  Figure  13,  the  domain  of  the  basis  functions  is  the  cavity 
surface  S,  while  the  domain  of  the  testing  functions  is  the  aperture  a.  The  elements  of  Z^  are 
found  from  the  iimer  product 


mn^  fmn‘ 


f^j  X  -riAri;?,?')  }  dads' 


(90) 


where  m  =  0,  ±1,  ±2, ...  ,  ±M,  1  <  n  <  Na,  and  Na  <  n'  <  Na  +  Ng  —  1.  Noting  that 
r(A:;  r,  r ')  can  be  written  as 


r(A;;r,r')  =  -yfcVG(A:;  r,  r ')  x  I 


(91) 
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and  using  the  identity  A  •  x  C j  =  [^A  x  Bj  ■  C,  we  can  write  equation  (90)  as 

=  dads'  (92) 


Making  use  of  equation  (85)  in  equation  (90)  we  obtain 


mn»  U 


n  = -11=1  ■  {n' 


**(r)xR]} 

Q-jkR 


(93) 


where  R  =  r  —  r Although  the  integrand  in  equation  (93)  has  a  0{R~^)  singularity  at  r  =  r 
since  the  aperture  a  and  cavity  surface  S  do  not  overlap,  the  integral  converges.  We  now  calculate 
each  of  the  four  permutations  of  the  vector  portion  of  the  integrand  of  equation  (93).  We  define  the 
relevant  vectors  as 


I  =  X  cos  0  +  y  sin  ^  (94a) 

^  =  —X  sm(f)  +  y  cos  0  (94b) 

/  =  X  cos  a'  cos  (j)'  +  y  cos  a'  sin  <j)'  —  z  sin  a'  (94c) 

=  —X  sin  +  y  cos  (f)'  (94d) 

R  =  X  (p cos (f)  —  p' cos 4>)  +  y  (p sin (f>  —  p' sin (j)')  +  zz'  (94e) 


from  which  we  find 

i'  -r^  X  {ix  R^  =  —z'  cos(^  —  (f)')  (95a) 

-n'  X  X  R^  =  —  (p'  sin  a'  +  z'  cos  a')  sin{(f>  —  ({>')  (95b) 

t'  -ii!  X  (^x  R^  =  z'  sin(0  —  </>')  (95c) 

-n'  X  X  R^  =  [p  —  p^  cos(^  —  ^')]  sin  a'  —  z'  cos  a'  cos(<^  —  (j)')  (95d) 

4.4.8  Z^.  is  the  2Ns  x  2Ns  submatrix  associated  with  the  testing  of  equation  (37) 

and  the  expansion  of  J5:.  As  shown  in  figure  14,  the  domain  of  both  the  basis  and  testing  functions 
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Figure  14  Domains  of  ip{f)  and  7(f ')  for  Z^.  Both  r  and  r '  are  defined  on  cone  frusta  formed 
by  revolving  a  portions  of  the  cavity  surface  generating  arc  about  the  z-axis. 


is  the  cavity  surface  S.  The  elements  of  are  found  from  the  inner  product 
=  •r2(fci;f,f')]}  dads' 

(?')*'  (96) 

where  m  =  0,  ±1,  ±2, . . .  ,  ±M,  Na  <  n  <  Na  +  Ng  —  1,  and  Na  <  n'  <  Na  +  Ng  —  1.  Noting 
that  r2(A:;  r,  r ')  can  be  written 


r2(A:;f,r')  =  -jk  [vC?(fc;f,r')  x  I  -  VG(fc;  r,  f/)  x  1^] 
and  using  the  identity  A  •  j^B  x  cj  =  j^A  x  fij  •  C,  we  can  write  equation  (96)  as 

^  X  VG'(A:i;  f,  ? ')]  }  ds  ds' 

J  S  J  s 

~  f  [  X  X  VG(/i;i;r,r/)l  }  ds  ds' 

^  S  *'  s  ^ 


(97) 


(98) 
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We  may  express  the  gradients  of  the  scalar  Green’s  functions  as 


VG{k-,  r,  r ')  =  {-jk  |r  -  r '|  -  l) 
VG{k;r,Ti)  =  (-jA;|r-r/[  -  l) 


exp(-jfc|r-r'|) 

Atx  I?  —  r'l^ 
exp(-jA:|r-r/|) 
47r|r-ri'p 


(r-r') 

(?-?/) 


(99a) 

(99b) 


We  have  defined  R  =  |r  -  ?'[  and  R  =  f  -  f';  we  will  now  define  Ri  =  |r  -  r/l  and  ^  = 
r  —  ?/.  Note  that  this  definition  of^  is  a  deviation  from  the  convention  that  the  subscript  i  on 
a  vector  indicates  the  image  of  that  vector!  These  quantities  are  shown  in  Figure  15.  We  rewrite 


Figure  15  Definition  of  Rj.  Ri  =  r  —  Fj'  where  r  denotes  a  point  on  S  and  ri  denotes  a  point 
on  Si. 


equation  (98)  as 


-3h  {iif  (r')  ■  [ii'  X  {*t(f)  X  r}]  }  ^-Jk^R  -  1)  ^  *  *' 

+]h  ^  {ttn’if)  ■  [*'  X  X  R,}  J  }  {-jhRi  -  1)  ds  *'  (100) 
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Using  elementary  vector  analysis,  we  may  show  that  A  •  x  Cj j  =  —  Aj  •  x  C j  so  equa¬ 
tion  (100)  becomes 

Imn'i  ~ 

-jki  X  r}]  }  i-jhR  -  1)  ds  ds' 

-jki  /  [“i  ^  X  Rj}]  }  {-jkiR  -  1)  ds  ds'  (101) 

The  second  integral  in  equation  (101)  converges  because  Rj  >  0  for  f  G  5  and  r'  G  S.  The  third 
integral  also  converges  since  the  integrand  is  bounded  on  the  (finite)  region  of  integration.  The  first 
integral  in  equation  (101)  must  be  handled  carefully  because  of  the  (apparent)  0{R~^)  singularity 
of  the  integrand  at  r  =  r However,  we  can  show  that 

7(r ')  •  n'  X  |t^(r)  x  r}]  |  =  0{R^),  r  ->■  r (102) 

where  7(r')  is  either  £'  or  and  i^(r)  is  either  £  or  Thus  the  integrand  of  the  first  integral 
in  equation  (101)  has  only  a  weak  (D{R~^)  singularity  at  r  =  f '  and,  hence,  the  first  integral 
converges. 

4.5  Forcing  Function 

In  the  matrix  equation  Zmlm  =  the  right  side,  V^,  contains  the  effect  of  the  m*  Fourier 
mode  of  the  incident  field.  1^  is  a  column  vector  given  by  [ui  V2  . . .  where  N  =  ANa+2Ns 
is  the  system  size.  The  elements  of  Vm  correspond  to  the  right  sides  of  equations  (35)-(37).  The 
right  side  of  equation  (35)  is  zero,  which  causes  the  first  2Na  elements  of  Vm  to  identically  vanish. 
Likewise,  the  right  side  of  equation  (37)  is  zero,  causing  the  last  2Ns  elements  of  Vm  to  vanish. 
Only  the  elements  of  Vm  which  correspond  to  the  testing  of  equation  (36)  do  not  identically  vanish; 
these  are  elements  i  =  2Na  +  1  through  i  =  4Na.  The  effect  of  the  incident  field  is  manifested  via 
the  “physical  optics”  current  density  J°(r')  =  2z  x  H‘"®(r'),  r'  G  cr. 
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We  assume  that  the  system  excitation  is  caused  by  a  uniform  plane  wave  having  a  propaga¬ 
tion  unit  vector  k,  as  shown  in  Figure  16.  Without  loss  of  generality,  we  restrict  k  to  lie  in  the 


xz-plane,  which  we  will  call  the  plane  of  incidence.  Then  we  define  the  unit  vector  ex  to  be  per¬ 
pendicular  to  the  plane  of  incidence  and  the  unit  vector  ey  to  be  parallel  to  the  plane  of  incidence. 
Both  ex  and  ey  are  perpendicular  to  k,  and  satisfy  the  relation  k  =  ey  x  ex-  Defining  9  as  the 
angle  between  z  and  —  k,  we  find 


=  —X  sin  9  —  i  cos  9 

(103a) 

=  y 

(103b) 

=  —X  cos  9 +  2.  sui9 

(103c) 

We  assume  that  the  electric  field  vector  E*"®  of  the  incident  plane  wave  has  unit  amplitude 
and  zero  phase  at  the  origin,  and  is  either  completely  parallel  to  the  plane  of  incidence  (||  polariza¬ 
tion)  or  completely  perpendicular  to  the  plane  of  incidence  (-L  polarization): 


where  e  =  e||  for  parallel  polarization  and  e  =  for  perpendicular  polarization.  Then  the 
incident  magnetic  field  Tb  k  x  can  be  written  as 


(105) 


where  e  =  ej.  for  parallel  polarization  and  e  =  — ey  for  perpendicular  polarization.  The  “physical 
optics”  current  density  J°(r')  =  2z  x  H‘''‘^(r '),  z'  =  0,  is  given  by 


_  2  yjj  gjfcop' sine  cos  1^'  J 


— X 


y  COS0 


for  II  Pol 
for  ±  Pol 


(106) 


where  p'  =  \Jx'^  +  y'^  and  (j)'  =  axctsji{y' Jx'). 

To  find  Vi,  2Na  +  1  <  *  <  4Na,  we  calculate  the  inner  product  of  J°  and  corresponding 
testing  function  (r '); 


Vi  = 


jkpZo 

2 

jkoZo 


/jO  \ 
(•'  5  Tmn' ' 


J<x 


(107) 


where  1  <  n'  <  Na-  The  testing  function  was  introduced  in  Section  4.3;  we  repeat  it 

here  for  convenience. 


(108a) 

(I08b) 

Since  the  testing  functions  that  contribute  to  V  are  defined  only  for  r '  G  a,  the  arc-length  variable 
£'  is  equivalent  to  the  polar  radial  variable  p{£'),  and  we  can  write  £'  =  xcos^'  -1-  y  sini^'  and 
0  =  —X  sin  <f>'  +  y  cos  (j)'.  Expressing  the  inner  product  as  an  integral  in  polar  coordinates,  we  have 

^27r  f^Na 

/  f{e',cl>')9{e,4>')di'd<t>'  (109a) 

Jo  Jo 
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where 


=  2  Fo  eJ^o^'sinScos^' 

—  COS  </>'  for  7  =  7^  and  ||  Pol 
sin  (j)'  for  7  =  7'^  and  1 1  Pol 

9{0,^')  =  \ 

cos  0  sin  (j)'  for  7  =  7^  and  J.  Pol 
cos  6  cos  (/>'  for  7  =  7"^  and  J.  Pol 

7-"  , 

Using  basic  trigonometric  identities  and  the  relation  Jn{z)  = -  /  cos(n0)  dcj),  where 

TT  70 

Jn  is  the  Bessel  function  of  the  first  kind  of  order  n,  we  can  write  the  inner  product  as 

_  .  ,  r^Na 

<J°>oit)  =  2  7ry'o  /  Tn>{e')  fmin  di'  (110a) 

Jo 

where 

7  5m  ^  ')  for  7  =  7^  and  II  Pol 

for  7  =  7'^  ^0  II  Pol 

un  =  { 

—  cos  6  gm  {£')  for  7  =  7^  and  ±  Pol 
—7  cos  0  5m  ^  0  for  7  =  7*^  and  ±  Pol 
5m  ^(^')  =  7”*  {Jm-i(M'sin5)  ±  Jm+i(M'sin5)} 

We  note  that  5^^  =  5^^  and  5^^  =  -5^^- 
4.6  Symmetry 

Up  to  this  point,  we  have  been  working  with  the  matrix  equation  Zm  Im  =  Pm  where 
m  =  0,  ±1, . . .  ±  M.  It  so  happens  that  we  need  only  explicitly  solve  the  matrix  equation  for  m  = 
0, 1, . . .  M.  The  solution  vector  Im  for  negative  m  can  be  found  with  no  additional  computation 
by  exploiting  the  symmetry  between  Zm /m  =  PmandZ_mf-m  =  VLm- We  may  symbolically 


(110b) 


(110c) 


(109b) 


(109c) 
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write  the  matrix  equation  Zm  Im  =  Vm  as 


In  equation  (111),  the  block  submatrix  contains  the  elements  of  corresponding  to  £- 

directed  basis  functions  and  ^-directed  testing  functions;  the  other  block  submatrices  are  defined 
analogously.  The  block  subvector  contains  the  expansion  coefficients  of  equation  (46)  for  i- 
directed  basis  functions;  the  other  subvectors  of  Im  are  defined  analogously.  The  block  subvector 
contains  the  result  of  testing  equation  (36)  with  ^-directed  testing  functions;  Vm  is  defined 
analogously. 

Each  submatrix  of  Zm  is  either  odd  or  even  with  respect  to  the  Fourier  index  m.  By  exam¬ 
ining  the  formulations  in  Sections  4.4. 1^.4. 8,  we  find  that 

'fs  0\  fo  fo 
\0  e)  \e  o)  o)  \At) 

fo  fE  o\  /e 
\E  o)  e)  e)  \Bi) 

k  0\  fE  o\  /e  o\  /ci\ 

\0  0J  eJ  e)\  [[cij 

Equation  (112)  repeats  equation  (111)  except  that  each  submatrix  of  Zm  in  equation  (1 1 1)  has  been 
replaced  with  an  “E”  if  it  is  even  with  respect  to  m,  an  “O”  if  it  is  odd  with  respect  to  m,  and  an 
“0”  if  it  is  identically  zero.  From  section  4.5  we  see  that  and  =  —V^m  for  parallel 

polarization,  while  =  -V-m  =  Vlm  for  perpendicular  polarization.  From  this  we 
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may  draw  the  following  conclusions  regarding  the  subvectors  of  Im- 

B*  =  (113) 

Ot,  = 

where  the  top  sign  is  taken  for  parallel  polarization  and  the  bottom  sign  is  taken  for  perpendicular 
polarization.  Exploiting  the  relations  of  equation  (113),  we  may  then  re-write  equations  (46)-(48) 


to  sum  over  only  non-negative  Fourier  indices: 

M  Na  rp  fp\ 

M,(^,  «  -^0  E  E  /(</>)  ^  +  Atn  9{<f>)  i]  (114a) 

m=0  n=l  ’ 

M  Na  rp  /p\  . 

</>)  «  E  E  ^  +  BL  f{<p)  ^  (114b) 

m=0  n=l  ^1  1 

Na'^Ns  1  rp  /^\  -  . 

</')  «  E  S  9{<I>)^  +  Cir,f  {<!>)(]  (114c) 

m=0  n=Na  ' 

where  0  <  <^  <  27r,  0  <  £  <  ^iVa  for  and  3^,  iNa  <  ^  <  ^Na+N^  for  J^,  and 

/(^)  =y  sin(m^)  (115a) 

g{<j))  =  cos{m(f>)  (115b) 

<• 

1  for  m  =  0 

em=<  (115c) 

2  for  m  7^  0. 


Cl  =  ±C^-m 
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V.  Numerical  and  Experimental  Results 

A  computer  program  was  written  to  implement  the  method  of  moments  solution  outlined  in  Chap¬ 
ter  rv  for  the  BOR  OCRED  scattering  problem.  For  simplicity,  the  program  was  written  only  for 
OCRED  geometries  in  which  the  cavity  volume  is  unfilled.  As  noted  previously,  when  the  cavity 
volume  is  unfilled,  only  the  electric  current  density  on  the  aperture  (Jo-)  and  the  electric  current 
density  on  the  cavity  surface  (J5)  need  to  be  found.  Once  Ja  and  Js  are  determined,  the  program 
calculates  the  monostatic  radar  cross  section  for  the  user-supplied  angle  of  incidence. 

To  validate  the  results  generated  by  the  program,  several  physical  BOR  OCRED  models 
were  fabricated  and  their  radar  cross  sections  measured.  The  measurements  are  discussed  in  the 
next  section,  and  the  results  for  each  model  are  discussed  in  subsequent  sections. 

5.1  BOR  OCRED  Geometries 


(a)  Spherical  OCRED.  £0  =  0,  (b)  Cylindrical  OCRED.  £0  =  0,  (c)  Conical  OCRED.  £0  =  0, 

£l=Ra,£2  =  Ra+Rs{TT-a).  £1  =  Ra,  £2  =  Ra  +d.  £3  =  £i  =  R^,  £2  =  R^  + . 

2Ra  -I-  d.  There  are  two  cylindri¬ 
cal  OCREDs,  deep  and  shallow. 

Figure  17  Generating  arcs  for  the  four  BOR  OCRED  geometries  investigated  in  this  dissertation. 
The  cavity  volume  D  is  carved  out  by  rotating  the  generating  arc  about  the  2:-axis.  The  dashed 
portion  forms  the  aperture  a  while  the  solid  portion  forms  the  cavity  surface  S. 

Figure  17  shows  the  generating  arcs  for  the  four  BOR  OCRED  geometries.  Each  BOR 
OCRED  geometry  is  formed  by  revolving  the  generating  arc  about  the  ^r-axis.  The  portion  of  the 
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generating  arc  on  the  p-axis  generates  the  OCRED  aperture  a,  while  the  portion  below  the  p-axis 
generates  the  cavity  surface  S.  The  ground  plane  (the  z  =  0  plane  less  the  aperture  a)  is  implied. 

Each  generating  arc  shown  in  Figure  17  can  be  mathematically  defined  in  terms  of  the  arc 
length  parameter  I,  which  measures  the  distance  along  the  generating  arc,  starting  from  the  origin. 
We  formally  define  each  generating  arc  mathematically  as  follows. 


Spherical  OCRED.  The  generating  arc  for  the  spherical  BOR  OCRED  is  shown  in  Fig¬ 
ure  17(a).  The  equations  defining  the  generating  arc  are 


0<i<Ra 

p{tj  =  i 

\ 

(116) 

Rg  sin  ( 

Ra<£<Ra  +  Rs{t^  -  a) 

* 

0, 

0<£<Ra 

z(£)  -  < 

Rg cos  ^ 

—  Rg  cosa,  Ra  <£  <  Ra  +  -Rs(7r  —  a) 

(117) 

where  Rg  =  7/8  inches  =  0.022225  meters  is  the  radius  of  the  spherical  cavity,  Ra  =  5/8  inches  = 
0.015875  meters  is  the  radius  of  the  aperture,  and  a  —  arcsin(i?o/i?s)  «  45.6°.  The  domain  of 
p{t)  and  z{i)  is  0  <  ^  <  i?a  +  Rs{t^  —  «)• 


Cylindrical  OCREDs.  The  generating  arc  for  the  generic  cylindrical  BOR  OCRED  is  shown 
in  Figure  17(b).  Two  cylindrical  BOR  OCREDs  are  considered  in  this  research,  a  shallow  one  and 
a  deep  one.  The  equations  defining  the  generating  arc  are 


i,  0<e<Ra 

Pi^)  ^  Rai  Ra  ^  ^  <  Ra  +  d 

2R(i  -|-  ol  —  £,  Ra  -f-  d  ^  ^  2/{q  -|-  d 

f 

0,  o<e<Ra 

~  ^  Ra  ~  Ra  ^  ^  <  Ra  +  d 

— d,  Ra  -f"  d  ^  ^  ^  -|-  d 


(118) 


(119) 
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where  Ra  =  5/8  inches  =  0.015875  meters  is  the  radius  of  the  aperture,  and  d  is  the  depth  of  the 
cylindrical  cavity.  For  the  deep  cylindrical  OCRED,  d  =  3  inches  =  0.0762  meters,  while  for  the 
shallow  cylindrical  OCRED,  d  =  1.5  inches  =  0.0381  meters.  The  domain  of  p{i)  and  z{i)  is 
0  <e<  2Ra  +  d. 


Conical  OCRED.  The  generating  arc  for  the  conical  BOR  OCRED  is  shown  in  Figure  17(c). 
The  equations  defining  the  generating  arc  are 


P{f)  =  { 


0<e<Ra 


i?„|l-^=£i=j,  Ra<i<Ra  +  yM+^ 


0, 


0<e<Ra 


zii)  =  <  p_Ji  ^ _ 

-d-=^,  Ra<e<Ra  +  ^^WT^ 


(120) 


(121) 


where  Ra  =  I  inch  =  0.0254  meters  is  the  radius  of  the  aperture,  and  d  =  1.75  inches  =  0.04445 
meters  is  the  depth  of  the  conical  cavity.  The  domain  of  p{£)  and  z{£)  isO  <  £  <  Ra  +  y/Rl  +  (P. 


5.2  Measurements 

A  physical  model  of  each  BOR  OCRED  geometry  described  in  Section  5.1  was  created  for 
the  purpose  of  measuring  its  RCS  at  various  frequencies,  polarizations,  and  aspects.  Each  BOR 
OCRED  model  was  machined  out  of  a  solid  block  of  aluminum  and  mounted  in  a  square,  1/16  inch 
thick  aluminum  plate.  The  models  are  shown  in  Figure  18.  While  the  conductivity  of  aluminum 
is  not  infinite,  it  is  sufficiently  high  to  model  the  ideal  perfectly  conducting  nature  of  the  OCRED 
cavity  surface  and  ground  plane. 

Measurement  Process.  Each  OCRED  model  was  then  mounted  in  a  test  fixture,  as  shown  in 
Figure  19.  A  metallic  adapter  was  used  to  transition  from  the  7-inch  square  OCRED  model  to  the 
30-inch  diameter  hole  in  the  test  fixture.  Metallized  tape  was  used  to  ensure  electrical  continuity 
between  the  adapter  and  the  OCRED  model,  and  between  the  adapter  and  the  test  fixture.  The 
test  fixture  had  a  large,  flat  surface  which  was  metallized  to  simulate  an  infinite  conducting  ground 
plane,  and  was  specially  shaped  to  minimize  its  monostatic  radar  cross  section. 
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(c)  Conical  OCRED  model. 


(d)  Shallow  Cylindrical  OCRED  model. 


Figure  18  OCRED  models. 


The  RCS  of  each  OCRED  model  was  measured  in  Wright  Laboratory’s  Advanced  Compact 
Range.  Figure  20  shows  the  test  fixture,  with  OCRED  model  installed,  mounted  on  the  target  sup¬ 
port  pylon.  The  RCS  of  the  OCRED  model  is  measured  using  a  coherent  background  subtraction 
technique  [33].  This  technique  consists  of  measuring  the  far-zone  scattered  electric  field  under 
four  separate  configurations  of  the  compact  range.  In  the  first  configuration,  the  scattered  electric 
field  of  the  test  fixture  with  the  BOR  OCRED  model  installed  is  measured  (£;OCRED+tf^. 
second  configuration,  the  scattered  electric  field  of  the  test  fixture  with  the  OCRED  aperture  re¬ 
placed  with  a  conductor  is  measured  in  the  third  configuration,  the  scattered  electric  field 

of  a  metal  sphere  is  measured  and  in  the  fourth  configuration,  the  scattered  electric  field 
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Figure  19  The  test  fixture  with  an  OCRED  model  installed.  The  OCRED  model  is  mounted  in 
an  adapter  having  a  30-inch  diameter,  which  is  in  turn  mounted  in  the  test  fixture.  Metallic  tape  is 
used  at  the  junctions  to  ensure  electrical  continuity. 


Figure  20  The  test  fixture  mounted  on  the  support  pylon  in  the  Advanced  Compact  Range. 
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Figure  21  The  adapter  with  one  of  the  OCRED  models  installed.  The  other  three  models  are 
shown  in  the  foreground. 


of  the  empty  chamber  is  measured  chamber^  jjig  jj^g  qcRED  is  computed  as 


croCRED  = 


^OCRED+TF  _  ^TF 
^Sphere  _ j^Empty  Chamber 


2 

•^Sphere 


(122) 


where  csphere  is  the  calculated  RCS  of  the  conducting  sphere.  The  quantities  in  equation  (122) 
are  implicit  functions  of  aspect  angle,  polarization,  and  frequency.  The  Advanced  Compact  Range 
measures  these  quantities  to  within  0.25  dB  of  their  true  values  over  the  frequency  range  of  2-18 
GHz.  Neglecting  the  presence  of  noise  in  the  receiver,  this  measurement  technique  is  exact  pro¬ 
vided  the  OCRED  and  the  test  fixture  do  not  interact,  and  its  accuracy  degrades  as  such  interactions 
become  more  significant.  The  nature  of  this  degradation  is  as  follows. 

Figure  22  illustrates  how  the  background  subtraction  process  works.  and  £;OCREd+tf 
are  the  complex  quantities  (scattered  electric  fields)  directly  measured  in  the  compact  range,  and 
their  difference  is  ^ocred  ^g  ^Quid  really  like  to  measure  £;0CRED  ^  ^tf^  ^g 

find  the  true  value  £;OCREd  _  ^^ocred  ^tf^  _  giv  ^ocred+tf  (jjffgrs  from  + 

E^^  by  an  amount  £!interaction^  which  is  caused  by  the  interactions  between  the  test  fixture  and 


Figure  22  Complex  (phasor)  representation  of  the  scattering  from  the  OCRED  model  and  the 
test  fixture.  is  the  scattered  electric  field  from  the  test  fixture;  is  the  scattered  electric 

field  from  the  OCRED  model;  £:OCREd+tf  jg  scattered  electric  field  from  the  test  fixture  with 
the  OCRED  model  installed;  is  the  scattered  electric  field  from  the  interactions  between 

the  test  fixture  and  the  OCRED  model;  and  E^^^^  is  the  difference  between  E^  and  £;OCRED+tf 

the  OCRED  model.  It  is  generally  not  possible  to  directly  measure  the  amplitude  or  phase  of 
^interaction  However,  if  is  small  relative  to  £;0CRED  +  ^tf^  ^jOCRed  +  ^^tf  ^ 

^OCRED+TF  ^OCRED  ^  ^OCRED  Generally,  £j'"'eraction  jg  largest  near  grazing  incidence  for 
parallel  polarization. 

A  problem  that  can  cause  errors  in  the  measurement  of  ^^cred  occurs  when  E^^  is  very 
much  larger  than  £:Ocred  (.j^jg  g^gg^  ^jjg  calculation  of  ^ocred  _  ^ocred+tf  _  ^tf  becomes 
numerically  sensitive,  such  that  small  uncertainties  in  measurement  can  produce  large  uncertainties 
in  J;OCRED  figure  23  shows  the  bounds  on  the  uncertainty  in  calculating  =  j5;OCRED+tf  _ 

E^^  as  a  function  of  the  relative  magnitude  of  £;0CREE)  ^tf  ^g  ggg  ^bat  even  though  £:0CRED 
and  E^^  are  measured  to  be  within  0.25  dB  of  their  true  values,  an  uncertainty  of  1 .0  dB  in  ^ocred 
requires  E^  to  be  at  least  3  dB  below  £;0CRED  jfig  situation  is  encountered  at  near-normal 
incidence  when  scattering  from  the  flat  portion  of  the  test  fixture  is  very  large. 

5.2.1  The  Spherical  OCRED.  The  first  OCRED  model  to  be  considered  is  spherical  in 
shape.  Its  generating  arc  is  shown  in  Figure  17(a)  and  described  in  Section  5.1.  Because  the  spher¬ 
ical  OCRED  bears  a  strong  resemblance  to  a  spherical  cavity,  we  expect  the  frequency  response  of 
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Figure  23  Uncertainty  in  estimating  ^ocred  ^ocred+tf  _  ^tf  Uncertainty  of  £;0Cred+tf 
and  is  0.25  dB. 

the  RCS  of  the  spherical  OCRED  to  exhibit  singular  behavior  near  the  resonant  frequencies  of  a 
spherical  cavity  of  the  same  radius.  Indeed,  we  can  think  of  the  spherical  OCRED  as  a  spherical 
cavity  with  an  aperture  that  serves  to  reduce  the  Q  (quality  factor)  of  the  cavity.  Using  a  classical 
separation  of  variables  approach,  the  lowest  resonant  frequencies  of  a  spherical  cavity  having  a 
0.875  inch  radius  can  be  determined  [21,  page  271].  These  frequencies,  and  their  corresponding 
modes,  are  shown  in  Table  1.  (The  triply  subscripted  modes,  e.g.  TMpqr,  are  defined  by  the  r’*’  zero 
of  the  spherical  Bessel  function  (TE)  or  its  derivative  (TM)  of  order  q,  along  with  the  associated 
Legendre  function  Pq.) 

Table  1  Lowest  resonemt  frequencies  of  a  perfectly  conducting  spherical  cavity  with  radius  of  7/8 
inch.  Each  resonant  frequency  is  associated  with  multiple  modes,  each  associated  with  a  particular 
value  of  p,  which  is  allowed  to  vary  from  0  to  the  value  of  the  second  index. 


Modes 

iim 

TEpii 

TMp3i 

Resonant  Frequency  [GHz] 

5.89 

8.31 

9.65 

10.7 

The  measured  RCS  of  the  spherical  OCRED  model  is  plotted  in  Figure  24.  The  RCS  is 
plotted  versus  frequencies  ranging  from  2  to  18  GHz,  for  9  =30°,  45°,  and  60°  off  normal  and  for 
parallel  and  perpendicular  polarizations.  Several  features  are  evident.  First,  there  is  a  pronounced 
null  in  the  RCS  near  6  GHz  at  all  angles  and  polarizations.  This  null  is  likely  associated  with  the 
resonant  frequency  of  the  TMpu  mode  of  the  spherical  cavity  (see  Table  1).  The  null  is  shifted 
about  5%  due  to  the  effects  of  the  aperture.  (If  the  aperture  were  smaller,  that  is  if  a  in  Figure  17(a) 
were  less  than  45°,  then  we  would  expect  the  null  to  be  shifted  back  toward  the  TMpn  resonant 
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Figure  24  The  measured  RCS  of  the  spherical  OCRED  model  as  a  function  of  frequency. 


frequency  since  the  Q  of  the  cavity  would  be  higher).  Second,  at  frequencies  above  7  GHz,  the  RCS 
at  different  angles/polarizations  appears  to  be  uncorrelated.  This  behavior  is  indicative  of  complex 
interactions  among  the  higher-order  cavity  modes  which  are  excited  differently  with  respect  to 
angle  and  polarization.  Third,  at  frequencies  below  5  GHz,  the  RCS  increases  nearly  monotonically 
for  all  angles/polarizations.  In  fact,  it  increases  proportional  to  the  square  of  the  electrical  area  of 
the  aperture.  For  example,  if  the  frequency  increases  by  a  factor  of  2,  then  the  electrical  area  of 
the  aperture  quadruples  and  the  the  OCRED  RCS  increases  by  20  log  4  =  12  dB.  Indeed,  this 
behavior  agrees  with  Rayleigh  scattering  [33]  in  which  RCS  is  proportional  to  the  fourth  power  of 
the  frequency. 

From  these  observations,  we  conclude  that  the  spherical  OCRED  RCS  falls  into  one  of  three 
categories,  depending  on  the  frequency  relative  to  the  lowest  resonant  frequency  fro  of  the  cor¬ 
responding  spherical  cavity.  When  the  frequency  is  significantly  below  fro,  then  the  OCRED 
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Table  2  Lowest  cutoff  frequencies  of  a  perfectly  conducting  cylindrical  waveguide  with  radius 
of  5/8  inch. 


Modes 

TEii 

TMoi 

TE21 

TMii 

Cutoff  Frequency  [GHz] 

5.54 

7.23 

9.19 

11.5 

behaves  as  a  Rayleigh  scattered  When  the  frequency  is  near  fro,  the  RCS  exhibits  a  deep  null  that 
is  nearly  independent  of  angle  and  polarization.  When  the  frequency  is  significantly  above  fro, 
multiple  cavity  modes  are  excited  in  the  OCRED  and  they  coherently  combine  to  create  an  RCS 
that  does  not  follow  an  easily  predictable  pattern. 

5.2.2  The  Deep  Cylindrical  OCRED.  The  second  OCRED  model  to  be  considered  is 
cylindrical  in  shape.  Its  generating  arc  is  shown  in  Figure  17(b)  and  described  in  Section  5.1.  For 
the  deep  cylindrical  cavity,  the  depth  is  c?  =  3.0  inches  =  0.0762  meters.  Because  the  cylindrical 
OCRED  bears  a  strong  resemblance  to  a  cylindrical  waveguide,  we  expect  the  frequency  response 
of  the  RCS  of  the  cylindrical  OCRED  to  exhibit  singular  behavior  near  the  cutoff  frequencies  of 
a  cylindrical  waveguide  of  the  same  radius.  Using  a  classical  separation  of  variables  approach, 
the  lowest  cutoff  frequencies  of  a  cylindrical  waveguide  having  a  0.625  inch  radius  can  be  de¬ 
termined  [21].  These  frequencies,  and  their  corresponding  modes,  are  shown  in  Table  2.  The 
measured  RCS  of  the  deep  cylindrical  OCRED  model  is  plotted  in  Figure  25.  The  RCS  is  plot¬ 
ted  versus  frequencies  ranging  from  2  to  18  GHz,  for  9  =30°,  45°,  and  60°  off  normal  and  for 
parallel  and  perpendicular  polarizations.  Several  features  are  evident.  First,  there  are  pronounced 
nulls  in  the  RCS  near  5.9  and  6.9  GHz  at  all  angles  and  polarizations.  These  nulls  are  associated 
with  the  group  velocities  of  the  lowest-order  mode  of  the  equivalent  cylindrical  waveguide  (see 
Table  2),  and  will  be  discussed  in  more  detail  in  Section  5.3.1.  Second,  at  frequencies  well  above  8 
GHz,  the  RCS  at  different  angles/polarizations  appears  to  be  uncorrelated.  This  behavior  is  indica¬ 
tive  of  complex  coherent  interference  among  the  higher-order  waveguide  modes  which  are  excited 
differently  with  respect  to  angle  and  polarization.  Third,  at  frequencies  below  5  GHz,  the  RCS  in¬ 
creases  nearly  monotonically  for  all  angles/polarizations,  again  displaying  the  Rayleigh  scattering 
behavior  seen  with  the  spherical  OCRED. 

From  these  observations,  we  conclude  that  the  deep  cylindrical  OCRED  RCS  falls  into  one 
of  three  categories,  depending  on  the  frequency  relative  to  the  lowest  cutoff  frequency  fco  of  the 


66 


T - 1 - 1 - 1 - 1 - 1 - r 


Frequency  [GHz] 

Figure  25  The  measured  RCS  of  the  deep  cylindrical  OCRED  model  as  a  function  of  frequency. 

corresponding  cylindrical  waveguide.  When  the  frequency  is  significantly  below  fco,  the  OCRED 
behaves  as  a  Rayleigh  scatterer.  When  the  frequency  is  near  fco,  the  RCS  exhibits  deep  nulls  near 
the  cutoff  frequencies  of  the  lowest-order  cylindrical  waveguide  modes.  These  nulls  are  nearly 
independent  of  angle  and  polarization.  When  the  frequency  is  significantly  above  fco,  multiple 
waveguide  modes  propagate  in  the  OCRED  and  coherently  combine  to  create  an  RCS  that  does 
not  follow  an  easily  predictable  pattern. 

5.2.3  The  Conical  OCRED.  The  third  OCRED  model  to  be  considered  is  conical  in 
shape.  Its  generating  arc  is  shown  in  Figure  17(c)  and  described  in  Section  5.1.  The  measured 
RCS  of  the  shallow  cylindrical  OCRED  model  is  plotted  in  Figure  26.  The  RCS  is  plotted  ver¬ 
sus  frequencies  ranging  from  2  to  18  GHz,  for  9  =30°,  45°,  and  60°  off  normal  and  for  parallel 
and  perpendicular  polarizations.  Unlike  the  previous  OCRED  geometries,  the  RCS  of  the  conical 
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OCRED  does  not  have  any  pronounced  nulls.  In  the  parallel  polarization,  the  RCS  of  the  con¬ 
ical  OCRED  is  remarkably  flat,  staying  within  5  dB  of  -22  dBsm  over  nearly  all  the  2-18  GHz 
bandwidth.  In  the  perpendicular  polarization,  the  RCS  for  6  =45°  and  9  =60°  off  normal  appear 
smooth  and  well  correlated,  while  the  RCS  for  0  =30°  exhibits  a  null  near  10  GHz.  This  null  is 
likely  caused  by  the  destructive  interference  between  the  scattering  from  the  near  and  far  edges  of 
the  aperture. 

5.2.4  The  Shallow  Cylindrical  OCRED.  The  fourth  OCRED  model  to  be  considered  is 
cylindrical  in  shape.  This  OCRED  model  is  the  same  as  the  second  one,  discussed  in  Section  5.2.2, 
except  the  depth  is  changed  to  d  =  1.5  inches  =  0.0381  meters.  Like  the  deep  cylindrical  OCRED, 
we  expect  the  frequency  response  of  the  RCS  of  the  shallow  cylindrical  OCRED  to  exhibit  singular 
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behavior  near  the  cutoff  frequencies  of  a  cylindrical  waveguide  of  the  same  radius.  These  cutoff 
frequencies  are  shown  in  Table  2. 


T - 1 - 1 - 1 - 1 - 1 - r 


Frequency  [GHz] 


Figure  27  The  measured  RCS  of  the  shallow  cylindrical  OCRED  model  as  a  function  of  fre¬ 
quency. 

The  measured  RCS  of  the  shallow  cylindrical  OCRED  model  is  plotted  in  Figure  27.  The 
RCS  is  plotted  versus  frequencies  ranging  from  2  to  18  GHz,  for  6  =30°,  45°,  and  60°  off  normal 
and  for  parallel  and  perpendicular  polarizations.  We  see  the  similar  features  as  those  for  the  deep 
cylindrical  OCRED,  except  that  the  null  at  approximately  5.9  GHz  is  absent.  It  will  become  clear 
in  Section  5.3.3  why  this  null  is  not  present,  and  we  shall  further  learn  that  the  location  of  the  null 
is  a  function  of  cylindrical  cavity  depth  d. 
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5.5  Calculations  -  RCS 


The  computer  program  implementing  the  numerical  algorithm  described  in  Chapter  IV  was 
run  on  the  four  BOR  OCRED  geometries  described  in  Section  5.1.  Since  the  algorithm  is  a 
frequency-domain  one,  each  computer  run  generated  data  for  a  single  frequency,  so  generating 
RCS  versus  frequency  plots  like  those  shown  in  Section  5.2  is  impractical.  Instead,  calculations 
for  each  geometry  were  made  at  several  judiciously  chosen  frequencies,  and  the  results  compared 
to  measurements.  In  all  cases,  the  moment  method  solution  is  calculated  using  20  segments  per 
free-space  wavelength. 

5.3.1  The  Spherical  OCRED.  The  spherical  OCRED  poses  a  special  challenge  in  that 
its  generating  arc  is  curved,  but  the  computer  program  was  written  to  handle  only  piecewise  linear 
generating  arcs.  Thus,  the  curved  part  of  the  generating  arc  has  been  replaced  by  a  piecewise 
linear  approximation,  as  shown  in  Figure  28.  Through  numerical  experiments,  it  was  found  that  an 
eleven-segment  piecewise  linear  approximation  was  adequate  for  the  highest  frequency  tested  (10 
GHz),  though  more  segments  might  be  necessary  at  higher  frequencies. 


Figure  28  A  piecewise  linear  approximation  to  the  generating  arc  for  the  spherical  OCRED. 
As  shown,  the  curved  part  of  the  generating  arc  is  approximated  by  three  line  segments.  In  the 
calculations,  eleven  line  segments  are  used. 

Figure  29  shows  the  resultant  calculated  RCS  for  the  spherical  OCRED  as  a  function  of 
angle  from  normal  at  three  discrete  frequencies.  The  corresponding  measured  RCS  values  are 
plotted  as  well,  and  good  agreement  is  observed  between  calculations  and  measurements  at  each 
frequency.  The  three  frequencies  chosen  are  4,  6.25,  and  10  GHz.  4  GHz  was  chosen  because  the 
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spherical  OCRED  behaves  as  a  Rayleigh  scatterer  at  that  frequency.  6.25  GHz  was  chosen  because 
measurements  show  that  the  RCS  of  the  spherical  OCRED  has  a  null  near  there  (see  Figure  24), 
which  implies  that  the  spherical  OCRED  is  near  resonance  at  6.25  GHz.  This  is  a  challenging 
frequency  for  the  calculations.  10  GHz  was  chosen  to  test  the  accuracy  of  the  moment  method 
program  with  an  electrically  large  problem. 

As  can  be  seen  in  Figure  29,  the  measured  and  calculated  values  agree  quite  well  at  each 
frequency  and  polarization.  Measured  data  only  exists  for  angles  from  9  =30°  to  9  =90°  from  nor¬ 
mal.  Measured  data  was  not  taken  at  normal  or  near-normal  aspects.  (The  background-subtraction 
measurement  technique  significantly  degrades  there  due  to  the  excessively  large  physical  optics 
return  from  the  test  fixture.)  The  RCS  at  4  GHz  is  nearly  constant  with  9  in  parallel  polarization, 
but  falls  off  rapidly  with  9  in  perpendicular  polarization,  a  direct  result  of  the  cos  9  factors  appear¬ 
ing  in  the  incident  field,  equation  (110),  and  the  radiation,  equation  (39).  At  6.25  GHz,  the  deep 
null  in  perpendicular  polarization  is  evident,  although  it  is  interesting  how  narrow  in  9  it  is.  The 
calculations  show  that  the  RCS  at  the  two  polarizations  converge  as  0  0,  as  expected.  The  mea¬ 

sured  data  is  significantly  corrupted  near  grazing  aspects  due  to  significant  interactions  between  the 
OCRED  and  the  edges  of  the  test  fixture.  As  expected,  this  corruption  is  most  evident  in  parallel 
polarization  and  decreases  as  frequency  increases. 

Figures  30  and  31  show  the  measured  and  calculated  RCS  of  the  spherical  OCRED  as  a 
function  of  frequency  and  angle  for  parallel  and  perpendicular  polarizations,  respectively.  The  fre¬ 
quency  varies  from  5.9  to  6.5  GHz  while  the  angle  varies  from  30°  to  85°  from  normal.  The  deep 
null  in  parallel  polarization  is  clearly  evident  (see  Figure  24),  and  we  note  that  the  location  of  the 
null  in  angle  varies  with  frequency.  The  measured  data  exhibits  contamination  near  grazing  inci¬ 
dence,  especially  in  parallel  polarization,  as  discussed  in  Section  5.2.  This  contamination  “ripples” 
with  frequency,  which  is  to  be  expected  since  the  electrical  distance  between  the  OCRED  model 
and  the  edge  of  the  test  fixture  varies  with  frequency. 

5.3.2  The  Deep  Cylindrical  OCRED.  The  generating  arc  (see  Figure  17(b))  for  the 
deep  cylindrical  OCRED  is  piecewise  linear,  and  thus  can  be  modeled  exactly  by  the  computer 
program.  Figure  32  shows  the  resultant  calculated  RCS  for  the  deep  cylindrical  OCRED  as  a 
function  of  angle  from  normal  at  four  discrete  frequencies.  The  corresponding  measured  RCS 
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Figure  29  The  RCS  of  the  spherical  OCRED  as  a  function  of  aspect  angle  9  at  three  discrete 
frequencies. 


veilues  are  plotted  as  well,  and  good  agreement  is  observed  between  calculations  and  measurements 
at  each  frequency.  The  four  frequencies  chosen  are  4,  5.9,  6.9  and  10  GHz.  4  GHz  was  chosen 
because  the  deep  cylindrical  OCRED  behaves  as  a  Rayleigh  scatterer  at  that  frequency.  5.9  and 
6.9  GHz  were  chosen  because  measurements  show  that  the  RCS  of  the  deep  cylindrical  OCRED 
has  nulls  near  those  frequencies  (see  Figure  25),  which  implies  that  the  deep  cylindrical  OCRED 
is  near  resonance  there.  These  are  challenging  frequencies  for  the  calculations  since  the  RCS  is 
very  sensitive  to  the  frequency  of  the  incident  field.  10  GHz  was  chosen  to  test  the  accuracy  of  the 
moment  method  program  with  an  electrically  large  problem  in  which  multiple  waveguide  modes 
are  expected  to  propagate  in  the  cavity. 

As  can  be  seen  in  Figure  32,  the  measured  and  calculated  values  agree  quite  well  at  each 
frequency  and  polarization.  Measured  data  only  exists  for  angles  from  6  =30°  to  6  =90°  from 
normal.  Measured  data  was  not  taken  at  normal  or  near-normal  aspects.  At  4  GHz,  the  RCS  is 
nearly  flat  for  parallel  polarization  and  decays  as  cos^  9  for  perpendicular  polarization.  At  5.9 
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Figure  30  RCS  versus  angle  from  normal  and  frequency  for  the  spherical  OCRED  for  parallel 
polarization.  Measured  RCS  plotted  on  left,  calculated  RCS  plotted  on  right. 


Figure  31  RCS  versus  angle  from  normal  and  frequency  for  the  spherical  OCRED  for  perpen¬ 
dicular  polarization.  Measured  RCS  plotted  on  left,  calculated  RCS  plotted  on  right. 

GHz,  both  polarizations  decay  more  rapidly  with  9  than  at  4  GHz.  At  6.9  GHz,  a  null  appears  at 
6  =  55°  in  parallel  polarization.  At  10  GHz,  the  roll-off  of  RCS  with  9  is  nearly  monotonic.  As 
was  the  case  for  the  sphericed  OCRED,  the  measured  data  is  significantly  corrupted  near  grazing 
aspects  due  to  significant  interactions  between  the  cavity  and  the  edges  of  the  test  fixture,  and  this 
corruption  is  most  evident  in  parallel  polarization  and  decreases  as  frequency  increases. 

5.3.3  The  Conical  OCRED.  The  generating  arc  for  the  conical  OCRED  (see  Fig¬ 
ure  17(c))  is  piecewise  linear,  and  thus  can  be  modeled  exactly  by  the  computer  program.  Fig¬ 
ure  33  shows  the  resultant  calculated  RCS  for  the  conical  OCRED  as  a  function  of  angle  from 
normal  at  two  discrete  frequencies.  The  corresponding  measured  RCS  values  are  plotted  as  well, 
and  good  agreement  is  observed  between  calculations  and  measurements  at  each  frequency.  The 
two  frequencies  chosen  are  4  and  10  GHz.  4  GHz  was  chosen  to  maintain  consistency  with  the 
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Figure  32  The  RCS  of  the  deep  cylindrical  OCRED  as  a  function  of  aspect  angle  6  at  four 
discrete  frequencies. 


other  OCRED  geometries.  10  GHz  was  chosen  to  test  the  accuracy  of  the  moment  method  program 
with  an  electrically  large  problem  in  which  significant  field  penetration  into  the  cavity  is  expected. 

As  can  be  seen  in  Figure  33,  the  measured  and  calculated  values  agree  quite  well  at  each 
frequency  and  polarization.  Measured  data  only  exists  for  angles  from  6  =30°  to  9  =90°  from 
normal.  Measured  data  was  not  taken  at  normal  or  near-normal  aspects.  At  4  GHz,  the  cos^  6 
dependence  is  evident  at  perpendicular  polarization.  The  RCS  at  normal  incidence  is  higher  than 
for  the  spherical  or  deep  cylindrical  OCREDs  due  to  the  larger  aperture  of  the  conical  OCRED.  At 
10  GHz,  there  is  a  fairly  pronounced  null  at  32.5°  in  perpendicular  polarization,  indicating  possible 
destructive  coherent  scattering  from  the  rim  and  wall  of  the  cavity.  As  with  previous  geometries, 
the  measured  data  is  significantly  corrupted  near  grazing  aspects  due  to  significant  interactions 
between  the  OCRED  and  the  edges  of  the  test  fixture,  and  this  corruption  is  most  evident  in  parallel 
polarization  and  decreases  as  frequency  increases. 
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Figure  33  The  RCS  of  the  conical  OCRED  as  a  function  of  aspect  angle  9  at  two  discrete  fre¬ 
quencies. 


5.3.4  The  Shallow  Cylindrical  OCRED.  The  generating  arc  for  the  shallow  cylindrical 
OCRED  (see  Figure  17(b))  is  piecewise  linear,  and  thus  can  be  modeled  exactly  by  the  computer 
program.  Figure  34  shows  the  resultant  calculated  RCS  for  the  shallow  cylindrical  OCRED  as 
a  function  of  angle  from  normal  at  four  discrete  frequencies.  The  corresponding  measured  RCS 
values  are  plotted  as  well,  and  good  agreement  is  observed  between  calculations  and  measurements 
at  each  frequency.  The  four  frequencies  chosen  are  4,  5.9,  6.9  and  10  GHz.  4  GHz  was  chosen 
because  the  shallow  cylindrical  OCRED  behaves  as  a  Rayleigh  scatterer  at  that  frequency.  6.9 
GHz  was  chosen  because  measurements  show  that  the  RCS  of  the  shallow  cylindrical  OCRED  has 
a  null  near  that  frequency  (see  Figure  27),  which  implies  that  the  shallow  cylindrical  OCRED  is 
near  resonance  there.  This  is  a  challenging  frequency  for  the  calculations.  5.9  GHz  was  chosen  as 
a  comparison  with  the  deep  cylindrical  OCRED.  As  noted  previously,  the  deep  cylindrical  OCRED 
RCS  has  a  null  at  5.9  GHz  that  is  curiously  absent  in  the  shallow  cylindrical  OCRED’s  RCS.  10 
GHz  was  chosen  to  test  the  accuracy  of  the  moment  method  program  with  an  electrically  large 
problem. 

As  can  be  seen  in  Figure  34,  the  measured  and  calculated  values  agree  quite  well  at  each 
frequency  and  polarization.  Measured  data  only  exists  for  angles  from  6  =30°  to  6  =90°  from 
normal.  Measured  data  was  not  taken  at  normal  or  near-normal  aspects.  At  4  GHz,  the  RCS  of  the 
shallow  cylindrical  OCRED  is  almost  identical  to  that  of  the  deep  cylindrical  OCRED,  suggesting 
the  RCS  is  independent  of  cavity  depth  at  that  frequency.  At  5.9  GHz,  the  RCS  has  nearly  the  same 
shape  as  for  the  deep  cylindrical  OCRED  but  is  nearly  12  dB  higher.  At  6.9  GHz,  the  RCS  has 
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Figure  34  The  RCS  of  the  shallow  cylindrical  OCRED  as  a  function  of  aspect  angle  9  at  four 
discrete  frequencies. 


nulls  in  both  polarizations,  but  at  slightly  different  angles.  At  10  GHz,  the  RCS  rolls  off  gently  in 
parallel  polarization,  and  has  a  null  in  perpendicular  polarization. 


5.4  Calculations  -  Other  Attributes 

In  Section  5.3,  the  ability  of  the  numerical  algorithm  to  accurately  model  the  radar  cross  sec¬ 
tion  of  several  BOR  OCRED  geometries  was  validated  against  measurements.  However,  the  RCS 
is  a  derived  quantity  inasmuch  as  it  depends  on  the  electric  current  density  on  the  cavity  surface. 
Therefore  it  is  prudent  to  evaluate  the  quality  of  the  (approximate)  electric  current  density,  which  is 
the  direct  product  of  the  numerical  algorithm.  Since  we  have  no  way  to  measure  the  electric  current 
density,  we  will  instead  check  that  the  numerical  algorithm  produces  electric  current  densities  3 a 
and  Js  that  satisfy  our  a  priori  knowledge  of  those  quantities.  We  will  also  evaluate  the  numerical 
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stability  of  the  moment  method  matrices  used  to  find  and  Js,  and  investigate  which  Fourier 
modes  must  be  included  in  the  solution  to  produce  accurate  RCS  data. 


5.4.1  The  Electric  Current  Density.  The  electric  current  density  J(r),  r  G  cr  U  S', 
will  exhibit  certain  properties.  As  noted  in  Chapter  III,  to  satisfy  Meixner’s  edge  condition  at  the 
aperture  rim,  the  £  component  of  J(r)  must  be  continuous  and  bounded  as  r  crosses  from  a  to  S, 
and  the  0  component  of  J(r)  must  be  discontinuous  and  unbounded  as  r  crosses  from  a  to  S.  At 
points  away  from  the  rim,  we  expect  J(r)  to  be  a  smooth  function  of  position. 

The  electric  current  density  near  the  center  of  the  aperture  is  a  continuous  function  of  position 
there,  and  so  can  be  expressed  in  a  Taylor  series  about  the  center  as 


J(r)=J(0)  +  O(|r|) 

=  xJa;(0) +yjj,{0) +  C>(|r|)  for  |r|  ^  0 


(123) 


However,  on  the  portion  of  the  aperture  associated  with  the  portion  of  the  generating  arc  between 
£o  and  ^i,  the  current  density  is  expressed  as  (see  equation  (47)) 

M 

J(f)=  YL  Kl«‘+«mi#] 

m=—M 
M 

=  E 
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y 

Equating  equations  (123)  and  (124),  we  find  that  and  vanish  for  |m|  ^  1  and  that 


Bml  siTKf) 


^mcj) 
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Bmi  sin  ^-5^1  cos(^ 


I 


for  |r|  <  li 


BLi  =  jBti  for  m  =  -1  (125a) 

BLi  =  -3Bti  form  =  l  (125b) 


These  relations  are  consistent  with  the  symmetry  relations  of  Section  4.6,  equation  (113).  Another 
way  of  stating  this  is  that  only  the  first-order  Fourier  modes  do  not  identically  vanish  at  ^  =  0. 
Using  an  identical  analysis  as  above  near  the  bottom  of  the  OCRED  cavity,  we  find  that  only  the 
first-order  Fourier  modes  of  J  do  not  identically  vanish  at  £  =  L. 
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In  the  rest  of  this  section,  we  examine  the  electric  current  density  as  a  function  of  the  arc 


length  variable  1.  We  define 


Mi)  =  1^'' 

to  be  the  m*  Fourier  components  of  the  ^-directed  and 
spectively. 


(126) 

(127) 

(^-directed  electric  current  density  J,  re- 


Figure  35  The  electric  current  density  plotted  as  a  function  of  arc  length  i  for  the  spherical 
OCRED;  /  =  4  GHz,  6  =  30°.  The  aperture  rim  corresponds  to  ^  =  1.59  cm. 


The  Spherical  OCRED.  At  a  frequency  of  4  GHz,  the  spherical  OCRED  behaves  as  a 
Rayleigh  scatterer,  as  can  be  seen  from  Figure  24.  A  characteristic  of  a  Rayleigh  scatterer  is  that  its 
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Figure  36  The  electric  current  density  plotted  as  a  function  of  arc  length  £  for  the  spherical 
OCRED:  /  =  4  GHz,  6  =  60°,  and  m  =  1.  The  aperture  rim  corresponds  to  ^  =  1.59  cm. 


scattering  cross  section  varies  sharply  with  respect  to  its  electrical  size,  but  not  with  respect  to  its 
shape.  This  implies  that  electromagnetic  fields  inside  the  cavity  are  small.  This  is  further  suggested 
by  the  fact  that  the  frequency  is  below  the  lowest  resonant  frequency  of  the  corresponding  closed 
spherical  cavity  (see  Table  1).  Figure  35  shows  the  magnitude  of  the  electric  current  density  on 
dD  for  6  =  30°  and  Fourier  modes  m=0,  1,  2,  and  3.  This  shows  that  all  Fourier  components  of 
the  electric  current  density  are  largest  near  the  rim  and  rapidly  decay  for  ^  As  expected, 

only  the  m  =  1  Fourier  component  has  a  non-zero  value  at  ^  =  0.  Figure  36  shows  the  magnitude 
and  phase  of  J(  and  for  6  —  60°  and  m  =  1.  Here  we  see  that  is  bounded  and  continuous 
near  ^  In  contrast,  the  magnitude  of  peaks  strongly  near  £  —  while  the  phase  of 

exhibits  a  180°  phase  shift  there.  This  behavior  is  entirely  consistent  with  the  edge  condition  at  the 
aperture  rim. 

At  6.25  GHz,  the  frequency  is  above  the  resonant  frequency  of  the  lowest-order  mode  of 
the  corresponding  closed  spherical  cavity.  Furthermore,  as  shown  in  Figure  24,  the  spherical  OC¬ 
RED  no  longer  behaves  as  a  Rayleigh  scatterer  at  6.25  GHz.  This  suggests  that  the  fields  should 
significantly  penetrate  the  cavity  interior.  Figure  37  shows  the  magnitude  of  the  electric  current 
density  components  for  6  =  30°  and  m=0,  1,  2,  and  3.  The  m=0  and  m=l  components  show  that 
there  is  indeed  significant  field  penetration  into  the  OCRED  interior.  However,  the  m=2  and  m=3 
(and  higher)  components  do  not  significantly  penetrate  the  cavity  interior,  indicating  that  only  the 
low-order  Fourier  modes  are  excited  at  this  frequency.  We  also  see  that  the  necessary  conditions 
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Figure  37  The  electric  current  density  plotted  as  a  function  of  arc  length  i  for  the  spherical 
OCRED:  /  =  6.25  GHz,  Q  =  30°.  The  aperture  rim  corresponds  to  £  =  1.59  cm. 


are  satisfied  at  ^  =  0,  £  =  ,  and  £  =  L.  That  is,  only  the  m=l  Fourier  component  is  non-zero  at 

the  aperture  center  and  cavity  bottom,  and  the  and  components  satisfy  the  edge  condition  at 
the  aperture  rim.  Figure  38  shows  the  magnitude  and  phase  of  and  for  9  =  60°  and  m  =  1. 
The  magnitude  and  phase  of  Je  are  continuous  and  bounded  at  £  =  while  the  magnitude  of 
J(f,  is  unbounded  and  the  phase  of  is  discontinuous  there.  Also,  note  the  null  and  180°  phase 
shift  of  at  ^  Ri  3.6  cm;  this  is  suggestive  of  a  standing  wave  pattern  caused  by  the  interference 
between  two  spherical  cavity  modes.  This  phenomenon  is  especially  likely  because  there  are  mul¬ 
tiple  modes  that  resonate  at  fro,  as  shown  in  Table  1.  Another  interesting  feature  is  the  “scalloping” 
evident  in  the  plot  of  J^.  This  is  caused  by  the  eleven-segment  piecewise  linear  approximation  to 
the  curved  portion  of  the  spherical  OCRED  generating  arc. 
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Figure  38  The  electric  current  density  plotted  as  a  function  of  arc  length  i  for  the  spherical 
OCRED:  /  =  6.25  GHz,  6  —  60°,  and  m  =  1.  The  aperture  rim  corresponds  to  ^  =  1.59  cm. 

At  10  GHz,  the  frequency  is  above  the  resonant  frequency  of  the  three  lowest-order  modes 
of  the  corresponding  closed  spherical  cavity  (see  Table  1).  We  expect  to  see  a  more  complex  field 
distribution  in  the  spherical  OCRED  than  was  evident  at  6.25  GHz  when  only  one  cavity  mode 
was  active.  Figure  39  shows  the  magnitude  of  the  electric  current  density  components  for  9  —  30° 
and  m=0,  1,  2,  and  3.  The  m=0,  1,  and  2  Fourier  components  penetrate  the  cavity  interior,  but  the 
m=3  component  does  not.  As  was  the  case  at  6.25  GHz,  all  the  necessary  conditions  are  satisfied 
at  ^  =  0,  ^  and  i  =  L.  Figure  40  shows  the  magnitude  and  phase  of  Je  and  for  9  =  60° 

and  m  =  1.  The  magnitude  and  phase  of  Ji  are  continuous  and  bounded  at  £  =  .  We  expect 

to  be  unbounded  at  ^  and,  indeed,  for  perpendicular  polarization  it  is,  but  not  for  parallel 

polarization.  However,  as  shown  in  Figure  41,  is  unbounded  in  magnitude  and  its  phase  shifts 
180°  at  ^  for  m=2.  Thus  we  see  that  the  total  obeys  the  edge  condition  at  £  = 

although  each  of  its  individual  Fourier  components  may  not.  As  was  the  case  at  6.25  GHz,  we  see 
the  distinctive  “scalloping”  effect  of  J(f,  due  to  the  eleven  segment  piecewise  linear  approximation 
to  the  curved  portion  of  the  spherical  OCRED  generating  arc. 

The  Deep  Cylindrical  OCRED.  At  a  frequency  of  4  GHz,  the  deep  cylindrical  OC¬ 
RED  behaves  as  a  Rayleigh  scatterer,  as  can  be  seen  from  Figure  25.  Table  2  shows  that  the  deep 
cylindrical  OCRED  will  not  support  a  propagating  waveguide  mode  at  4  GHz,  and  we  expect  to  see 
little  field  penetration  into  the  cavity  interior.  Figure  42  shows  the  magnitude  of  the  electric  of  the 
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Figure  39  The  electric  current  density  plotted  as  a  function  of  arc  length  t  for  the  spherical 
OCRED:  /  =  10  GHz,  0  =  30°.  The  aperture  rim  corresponds  to  £  —  1.59  cm. 


electric  current  density  components  for  9  =  30°  and  m=0, 1,  2,  and  3.  This  figure  indicates  that  all 
of  the  Fourier  components  of  Je  and  decay  for  ^  Additionally,  all  Fourier  components 

obey  the  edge  condition  at  the  aperture  rim,  emd  only  the  m=l  Fourier  component  does  not  vanish 
at  £  =  0  and  £  =  L. 

At  a  frequency  of  5.9  GHz,  the  deep  cylindrical  OCRED  will  support  the  TEn  cylindrical 
waveguide  mode  and  we  expect  to  see  significant  field  penetration  into  the  cavity  interior.  Figure  44 
shows  the  magnitude  of  the  electric  current  density  for  9  =  30°  and  m=0,  1,  2,  and  3.  Each 
Fourier  component  decays  for  ^  except  for  the  m=l  component  which  shows  considerable 

penetration  into  the  cavity  interior.  The  electric  current  density  for  m=l  is  evidently  caused  by 
an  interference  pattern  set  up  by  a  downward  propagating  TEn  wave  and  its  upward  propagating 
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Figure  40  The  electric  current  density  plotted  as  a  function  of  arc  length  I  for  the  spherical 
OCRED:  /  =  10  GHz,  Q  =  60°,  and  m  =  1.  The  aperture  rim  corresponds  to  ^  =  1.59  cm. 

Table  3  Longitudinal  wavelengths  the  four  lowest-order  modes  of  a  perfectly  conducting  cylin¬ 
drical  waveguide  with  radius  of  5/8  inch  at  5.9,  6.9,  and  10  GHz.  A  real  value  indicates  a  propa¬ 
gating  mode;  an  imaginary  value  indicates  an  evanescent  mode. 


Mode 

Az  [cm] 

f=5.9  GHz 

/=6.9  GHz 

/=10  GHz 

TEii 

14.7 

7.3 

3.6 

TMoi 

Jl.2 

7I3.8 

4.3 

TE21 

j4.3 

74.9 

7.6 

TMii 

73.2 

75.2 

reflection.  At  5.9  GHz,  the  “longitudinal  wavelength”  of  the  TEn  mode  is  given  by  Az  = 
27r/\/A:^  —  (1.841/i?a)^  =  14.7  cm  where  1.841  is  the  first  zero  of  J[  [21,  p.  205].  Table  3  shows 
the  “longitudinal  wavelength”  for  the  three  lowest-order  waveguide  modes  at  5.9,  6.9  and  10  GHz. 
This  value  is  very  nearly  twice  the  depth  of  the  deep  cylindrical  OCRED,  or,  in  other  words, 
the  OCRED  is  nearly  one-half  Az  at  5.9  GHz.  This  explains  the  deep  null  in  Figure  25  there. 
Figure  45  shows  the  magnitude  and  phase  of  the  electric  current  density  for  9  =  60°  and  m=l. 
A  similar  standing- wave  pattern  is  evident,  except  that  the  TEn  mode  is  excited  more  strongly  at 
9  —  30°.  Additionally,  the  and  components  are  nearly  ±90°  out  of  phase  for  ^  with 

Ji  experiencing  a  180°  phase  shift  at  its  null  near  £  =  5.5  cm. 

At  6.9  GHz,  the  deep  cylindrical  OCRED  still  propagates  only  the  TEn  cylindrical  waveg¬ 
uide  mode,  so  we  expect  to  see  a  field  distribution  in  the  cavity  interior  similar  to  that  at  5.9  GHz. 
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Figure  41  The  electric  current  density  plotted  as  a  function  of  arc  length  t  for  the  spherical 
OCRED:  /  =  10  GHz,  6  =  60°,  and  m  =  2.  The  aperture  rim  corresponds  to  ^  =  1.59  cm. 


Figure  46  shows  the  magnitude  of  the  electric  current  density  for  6  =  30°  and  m=0,  1,  2,  and 
3.  These  components  look  very  much  Idee  those  shown  in  Figure  44  for  5.9  GHz.  One  obvious 
difference,  however,  is  the  “double  hump”  pattern  of  the  m=l  Fourier  component  instead  of  the 
“single  hump”  pattern  present  at  5.9  GHz.  The  longitudinal  wavelength  of  the  TEn  mode  at  6.9 
GHz  is  \z  =  7.3  cm,  netirly  half  the  longitudinal  wavelength  at  5.9  GHz.  This  indicates  that  the 
deep  cylindrical  OCRED  is  nearly  one  longitudinal  wavelength  deep  at  6.9  GHz,  causing  the  deep 
null  in  Figure  25  there.  Figure  47  shows  the  magnitude  and  phase  of  the  electric  current  density 
for  9  =  60°  and  m=l.  The  same  “double  hump”  standing- wave  pattern  is  evident.  The  Jg,  and 
components  are  nearly  ±90°  out  of  phase  for  £  >  and  both  experience  a  180°  phase  shift  at 
their  nulls  {I  «3.8  and  7.4  cm  for  ,  i  «5.6  cm  for  J^). 

At  10  GHz,  the  deep  cylindrical  cavity  will  support  three  propagating  modes  (see  Table  2). 
Figure  48  shows  the  magnitude  of  the  electric  current  density  for  6  =  30°  and  m=0,  1,  2,  and  3.  At 
this  frequency,  the  three  lowest-order  waveguide  modes  (TEn,  TMqi,  and  TE21)  propagate,  each 
with  a  different  longitudinal  wavelength  as  given  in  Table  3.  Using  this  information,  we  see  that 
the  m  =  0  Fourier  component  contains  the  TMqi  mode,  the  m  =  1  Fourier  component  contains 
the  TEii  mode,  and  the  m  =  2  Fourier  component  contains  the  TE21  mode.  This  is  evident  from 
the  fact  that  the  null-to-null  distance  in  the  standing  wave  patterns  correlates  to  one-half  A^  of 
these  modes.  Figure  49  shows  the  magnitude  and  phase  of  the  electric  current  density  for  6  =  60°, 
m  =  1.  This  figure  shows  that  Ji  under  parallel  polarization  does  not  have  a  nicely  formed  lobing 
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Figure  42  The  electric  current  density  plotted  as  a  function  of  arc  length  t  for  the  deep  cylindri¬ 
cal  OCRED:  /  =  4  GHz,  6  =  30°.  The  aperture  rim  corresponds  Xo^—  1.59  cm. 

structure  near  the  top  of  the  OCRED  sidewall  i  +  1),  as  it  did  at  0  =  30°.  This  can 

be  attributed  to  the  TMn  waveguide  mode,  which  is  weakly  evanescent  at  10  GHz.  (Table  2  shows 
that  the  TMn  mode  cuts  on  at  11.5  GHz.)  Comparing  Figures  48  and  49  shows  that  the  TMn 
mode  is  excited  more  strongly  at  d  =  60°  than  at  6  =  30°,  probably  due  to  the  different  incident 
field. 


The  Conical  OCRED.  At  a  frequency  of  4  GHz,  the  conical  OCRED  is  in  transition 
between  the  Rayleigh  and  resonance  regimes.  This  implies  that  the  electromagnetic  field  does  not 
penetrate  as  strongly  into  the  cavity  as  at  higher  frequencies.  Figure  50  shows  the  magnitude  of 
the  electric  current  density  for  0  =  30°,  Tn=0-3.  As  expected,  none  of  the  Fourier  components 
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Figure  43  The  electric  current  density  plotted  as  a  function  of  arc  length  ^  for  the  deep  cylin¬ 
drical  OCRED:  /  =  4  GHz,  0  =  60°,  and  m  =  1.  The  aperture  rim  corresponds  to  ^  =  1.59 
cm. 


penetrates  strongly  into  the  cavity  interior,  and  all  of  them  obey  the  edge  condition  at  ^ 
Furthermore,  only  the  m  =  1  Fourier  component  does  not  vanish  at£  —  0. 

At  10  GHz,  the  conical  OCRED  does  not  behave  as  a  Rayleigh  scatterer,  and  we  expect  to 
see  significant  field  penetration  into  the  interior  of  the  cavity.  Figure  52  shows  the  magnitude  of  the 
electric  current  density  for  0  =  30°,  m=0-3.  In  contrast  to  4  GHz,  there  is  considerable  penetration 
into  the  cavity  interior,  especially  for  the  lowest  three  Fourier  components.  Figure  53  shows  the 
magnitude  and  phase  of  the  electric  current  density  for  9  =  60°,  m  =  1.  Notice  the  fairly  clean 
null  and  corresponding  180°  phase  shifts  at  £  w  4.1  cm  and  £  «  5.2  cm.  This  is  evidence  of 
a  well-formed  interference  pattern,  though  only  for  under  perpendicular  polarization  and 
under  parallel  polarization.  The  satisfaction  of  equation  (125)  at  ^  =  0  is  quite  visible  in  Figure  53. 

The  Shallow  Cylindrical  OCRED.  At  a  frequency  of  4  GHz,  the  shallow  cylindrical 
OCRED  behaves  as  a  Rayleigh  scatterer,  as  can  be  seen  from  Figure  27.  Table  2  shows  that  the 
shallow  cylindrical  OCRED  will  not  support  a  propagating  waveguide  mode  at  4  GHz,  and  we 
expect  to  see  little  field  penetration  into  the  cavity  interior.  Figure  54  shows  the  magnitude  of  the 
electric  current  density  components  for  9  =  30°  and  m=0, 1,  2,  and  3.  This  figure  indicates  that  all 
of  the  Fourier  components  of  Ji  and  decay  for  £  >  Additionally,  all  Fourier  components 
obey  the  edge  condition  at  the  aperture  rim,  and  only  the  m=l  Fourier  component  does  not  vanish 
at  ^  =  0  and  £  =  L. 
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Figure  44  The  electric  current  density  plotted  as  a  function  of  arc  length  t  for  the  deep  cylindri¬ 
cal  OCRED:  /  =  5.9  GHz,  Q  =  30°.  The  aperture  rim  corresponds  to  ^  =  1.59  cm. 


At  a  frequency  of  5.9  GHz,  the  shallow  cylindrical  OCRED  will  support  the  TEn  cylindrical 
waveguide  mode  and  we  expect  to  see  significant  field  penetration  into  the  cavity  interior.  Figure  55 
shows  the  magnitude  of  the  electric  current  density  for  0  =  30°  and  m=0,  1,  2,  and  3.  Each 
Fourier  component  decays  for  ^  except  for  the  m=l  component  which  shows  considerable 

penetration  into  the  cavity  interior.  The  electric  current  density  for  m=l  is  evidently  caused  by 
an  interference  pattern  set  up  by  a  downward  propagating  TEi  i  wave  and  its  upward  propagating 
reflection.  Table  3  shows  the  “longitudinal  wavelength”  for  the  three  lowest-order  waveguide 
modes  at  5.9  GHz,  from  which  we  see  that  =  14.7  cm  for  the  TEn  mode.  This  value  is  very 
nearly  four  times  the  depth  of  the  shallow  cylindrical  OCRED,  or,  in  other  words,  the  OCRED  is 
nearly  one-quarter  A^  at  5.9  GHz,  rather  than  one-half  Az  as  was  the  case  for  the  deep  cylindrical 
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Figure  45  The  electric  current  density  plotted  as  a  function  of  arc  length  t  for  the  deep  cylin¬ 
drical  OCRED:  /  =  5.9  GHz,  6  =  60°,  and  m  —  The  aperture  rim  corresponds  to  ^  =  1.59 
cm. 


OCRED.  This  explains  the  absence  of  a  null  in  Figure  27  there.  Figure  56  shows  the  magnitude 
and  phase  of  the  electric  current  density  for  6  =  60°  and  m=l.  The  same  standing- wave  pattern  is 
evident.  Additionally,  the  and  components  are  nearly  ±90°  out  of  phase  for  I  >  ,  with 

Ji  experiencing  a  180°  phase  shift  at  its  null  near  i  =  5.5  cm. 

At  6.9  GHz,  the  shallow  cylindrical  OCRED  still  propagates  only  the  TEn  cylindrical 
waveguide  mode,  so  we  expect  to  see  a  field  distribution  in  the  cavity  interior  similar  to  that  at 
5.9  GHz.  Figure  57  shows  the  magnitude  of  the  electric  current  density  for  6  =  30°  and  m=0, 
1,  2,  and  3.  These  components  look  very  much  like  those  shown  in  Figure  55  for  5.9  GHz.  One 
obvious  difference,  however,  is  the  “single  hump”  pattern  of  the  m=\  Fourier  component  instead 
of  the  “half  hump”  pattern  present  at  5.9  GHz.  The  longitudinal  wavelength  of  the  TEn  mode 
at  6.9  GHz  is  =  7.3  cm,  nearly  half  the  longitudinal  wavelength  at  5.9  GHz.  This  indicates 
that  the  shallow  cylindrical  OCRED  is  nearly  one-half  longitudinal  wavelength  deep  at  6.9  GHz, 
causing  the  deep  null  in  Figure  27  there.  Figure  58  shows  the  magnitude  and  phase  of  the  electric 
current  density  for  6  =  60°  and  m=l.  The  same  “single  hump”  standing-wave  pattern  is  evident. 
The  j£  and  components  are  necu-ly  ±90°  out  of  phase  for  £  >  and  experiences  a  180° 
phase  shift  at  its  null  at  £  «3.8  cm. 

At  10  GHz,  the  shallow  cylindrical  cavity  will  support  three  propagating  modes  (see  Table  2). 
Figure  59  shows  the  magnitude  of  the  electric  current  density  for  9  =  30°  and  m=0, 1, 2,  and  3.  At 
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Figure  46  The  electric  current  density  plotted  as  a  function  of  arc  length  I  for  the  deep  cylindri¬ 
cal  OCRED:  /  =  6.9  GHz,  6  =  30°.  The  aperture  rim  corresponds  to  ^  =  1.59  cm. 


this  frequency,  the  three  lowest-order  waveguide  modes  (TEn,  TMoi,  and  TE21)  propagate,  each 
with  a  different  longitudinal  wavelength  Az  as  given  in  Table  3.  Using  this  information,  we  see  that 
the  m  =  0  Fourier  component  contains  the  TMoi  mode,  the  m  =  1  Fourier  component  contains 
the  TEii  mode,  and  the  m  =  2  Fourier  component  contains  the  TE21  mode.  Figure  60  shows  the 
magnitude  and  phase  of  the  electric  current  density  for  9  =  60°,  m  =  1.  This  figure  shows  that 
j£  under  parallel  polarization  does  not  have  a  nicely  formed  lobing  structure  near  the  top  of  the 
OCRED  sidewall  <  ^  <  +  1),  as  it  did  at  9  =  30°.  This  can  be  attributed  to  the  TMn 

waveguide  mode,  which  is  weakly  evanescent  at  10  GHz.  Comparing  Figures  59  and  60  shows  that 
the  TMii  mode  is  excited  more  strongly  at9  =  60°  than  at  9  =  30°,  probably  due  to  the  different 
incident  field. 
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Figure  47  The  electric  current  density  plotted  as  a  function  of  arc  length  £  for  the  deep  cylin¬ 
drical  OCRED:  /  =  6.9  GHz,  6  =  60°,  and  m  =  1.  The  aperture  rim  corresponds  to  £  =  1.59 


Figure  48  The  electric  current  density  plotted  as  a  function  of  arc  length  £  for  the  deep  cylindri¬ 
cal  OCRED:  /  =  10  GHz,  9  =  30°.  The  aperture  rim  corresponds  to  £  =  1.59  cm. 
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Figure  49  The  electric  current  density  plotted  as  a  function  of  arc  length  £  for  the  deep  cylin¬ 
drical  OCRED:  /  =  10  GHz,  6  =  60°,  and  m  =  1.  The  aperture  rim  corresponds  to  £  =  1.59 
cm. 


Figure  50  The  electric  current  density  plotted  as  a  function  of  arc  length  £  for  the  conical  OC¬ 
RED:  /  =  4  GHz,  9  =  30°.  The  aperture  rim  corresponds  to  £  =  2.54  cm. 
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Figure  51  The  electric  current  density  plotted  as  a  function  of  arc  length  £  for  the  conical  OC¬ 
RED:  /  =  4  GHz,  6  =  60°,  and  m  =  1.  The  aperture  rim  corresponds  to  £  =  2.54  cm. 


Figure  52  The  electric  current  density  plotted  as  a  function  of  arc  length  £  for  the  conical  OC¬ 
RED:  /  =  10  GHz,  6  =  30°.  The  aperture  rim  corresponds  to  £  =  2.54  cm. 
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Figure  53  The  electric  current  density  plotted  as  a  function  of  arc  length  I  for  the  conical  OC¬ 
RED:  /  =  10  GHz,  0  =  60°,  and  m  =  1.  The  aperture  rim  corresponds  to  ^  =  2.54  cm. 


Figure  54  The  electric  current  density  plotted  as  a  function  of  arc  length  I  for  the  shallow  cylin¬ 
drical  OCRED:  /  =  4  GHz,  0  =  30°.  The  aperture  rim  corresponds  to  ^  =  1.59  cm. 
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Figure  55  The  electric  current  density  plotted  as  a  function  of  arc  length  i  for  the  shallow  cylin¬ 
drical  OCRED:  /  =  5.9  GHz,  Q  =  30°.  The  aperture  rim  corresponds  Xot  —  1.59  cm. 


Figure  56  The  electric  current  density  plotted  as  a  function  of  arc  length  t  for  the  shallow  cylin¬ 
drical  OCRED:  /  =  5.9  GHz,  Q  =  60°,  and  m  =  1.  The  aperture  rim  corresponds  to  ^  =  1.59 
cm. 
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Figure  57  The  electric  current  density  plotted  as  a  function  of  arc  length  £  for  the  shallow  cylin¬ 
drical  OCRED:  /  =  6.9  GHz,  6  =  30°.  The  aperture  rim  corresponds  to  £  =  1.59  cm. 


Figure  58  The  electric  current  density  plotted  as  a  function  of  arc  length  £  for  the  shallow  cylin¬ 
drical  OCRED:  /  =  6.9  GHz,  9  =  60°,  and  m  =  1.  The  aperture  rim  corresponds  to  £  =  1.59 
cm. 
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Figure  59  The  electric  current  density  plotted  as  a  function  of  arc  length  £  for  the  shallow  cylin¬ 
drical  OCRED:  /  =  10  GHz,  0  —  30°.  The  aperture  rim  corresponds  to  £  =  1.59  cm. 


Figure  60  The  electric  current  density  plotted  as  a  function  of  arc  length  I  for  the  shallow  cylin¬ 
drical  OCRED:  /  =  10  GHz,  0  =  60°,  and  m  =  1.  The  aperture  rim  corresponds  to  ^  =  1.59 
cm. 
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5.4.2  The  Condition  of  the  Moment  Method  Matrix.  The  effectiveness  of  the  moment 
method  solution  technique  requires  the  generalized  impedance  matrix  to  be  well-conditioned.  If 
the  matrix  is  well-conditioned,  then  small  changes  in  the  forcing  vector  cause  small  changes  in 
the  solution  vector,  and  the  linear  transformation  defined  by  the  matrix  is  stable.  If  the  matrix 
is  poorly  conditioned,  then  small  changes  in  the  forcing  vector  may  produce  large  changes  in  the 
solution  vector,  and  the  linear  transformation  defined  by  the  matrix  is  unstable.  The  condition  of  the 
moment  method  matrix  is  inherently  limited  by  the  integral  equations  to  which  the  moment  method 
is  applied,  and  a  well-conditioned  moment  method  matrix  is  an  indication  of  a  well-conditioned 
set  of  integral  equations. 

The  condition  of  an  iV  x  N  matrix  is  measured  by  the  so-called  “condition  number”  k.  It  is 
a  measure  of  the  sensitivity  of  the  linear  system  =  V  in  the  following  sense.  If  the  elements  of 
Z  and  V  are  in  error  by  some  small  amount,  then  the  resulting  solution  I  will  also  be  in  error,  such 
that 


{Z  +  SZ)  (/-h(52:)  =  (F-f5V) 


(128) 


The  error  of  the  solution  I  is  bounded  by  [18,  page  79] 


M<  JM  ,  Ml 
ii/||  -  Lmi  ii^iiJ 


+  0{5) 


(129) 


The  condition  number  k  is  thus  dependent  on  the  matrix  norm  employed.  If  the  2-norm  is  chosen, 
then  K  is  given  by  [18] 


On 

K  =  — 


0-1 


(130) 


where  ai  is  the  smallest  singular  value  of  the  matrix  and  ajv  is  the  largest.  (Recall  that  the  singular 
values  of  a  complex  matrix  are  always  real,  non-negative,  and  Op  <  a^+i.)  A  perfectly  conditioned 
matrix  {e.g.,  the  identity  matrix)  has  a  condition  number  k  =  1,  while  a  condition  number  much 
greater  than  unity  reflects  a  poorly  conditioned  matrix. 

Tables  4-7  show  the  condition  numbers  of  the  moment  method  matrix  used  to  calculate  the 
electric  current  density  J(r)  for  each  combination  of  OCRED  geometry,  frequency,  and  Fourier 
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Freq  [GHz] 

N 

Fourier  Index,  m 

0 

1 

2 

3 

4 

5 

6 

7 

4 

62 

111 

327 

rai 

131 

6.25 

86 

172 

368 

201 

199 

197 

10 

114 

219 

522 

368 

408 

294 

278 

273 

271 

Table  4  Condition  number  Kof  N  x  N  moment  method  matrix  Zm  for  the  spherical  OCRED. 


Freq  [GHz] 

N 

Fourier  Index,  m 

0 

1 

2 

3 

4 

5 

6 

7 

4 

74 

66 

191 

72 

72 

5.9 

100 

92 

241 

106 

104 

103 

6.9 

116 

110 

262 

129 

125 

124 

124 

10 

186 

159 

354 

198 

187 

178 

176 

175 

175 

Table  5  Condition  number  k  of  N  x  N  moment  method  matrix  Zm  for  the  deep  cylindrical 
OCRED. 

index.  It  should  be  stressed  that  these  condition  numbers  relate  to  using  the  moment  method  as 
applied  to  equations  (36)  and  (37)  to  find  and  3s  since  the  cavities  considered  in  this  chapter  are 
unfilled.  Several  trends  are  immediately  clear.  First,  the  condition  number  k  generally  increases 
with  frequency,  the  dependence  being  nearly  linear.  Second,  the  condition  number  corresponding 
to  m  =  1  exceeds  the  condition  number  for  any  other  Fourier  mode,  except  under  special  circum¬ 
stances  that  will  be  described  later  in  this  section.  Third,  as  the  Fourier  index  m  increases  past 
unity,  the  condition  number  asymptotically  approaches  a  finite  limit.  Most  importantly,  however, 
in  all  cases  the  condition  number  is  quite  benign  in  the  sense  that  it  does  not  vary  greatly  with 
geometry,  frequency,  or  Fourier  mode.  In  fact  the  largest  observed  k  was  632  while  the  smallest 
was  66,  a  variation  of  less  than  a  factor  of  ten. 

Figure  61  shows  the  condition  number  k  for  the  spherical  OCRED  as  a  function  of  frequency. 
Data  for  the  first  four  Fourier  modes  (m  =  0  —  3)  are  shown.  The  m  —  0  Fourier  mode  shows 
a  peak  in  its  condition  number  at  approximately  6.01  GHz  while  the  condition  numbers  for  the 
other  Fourier  modes  are  nearly  constant  across  the  band.  The  peak  in  the  m  =  0  Fourier  mode  can 
be  explained  as  follows.  The  lowest-order  mode  of  the  corresponding  spherical  cavity  resonates 
at  5.89  GHz  (See  Table  1).  The  presence  of  the  aperture  in  the  spherical  OCRED  causes  two 
things  to  happen:  the  resonant  frequency  becomes  complex  due  to  the  radiation  losses  through  the 
aperture,  and  the  volume  of  the  cavity  is  reduced  by  a  function  of  the  aperture  angle  a  and  the  cavity 
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Freq 

[GHz] 

N 

Fourier  Index,  m 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

4 

54 

150 

244 

109 

108 

108 

107 

10 

116 

377 

598 

263 

261 

267 

263 

257 

254 

253 

253 

253 

253 

Table  6  Condition  number  nof  N  x  N  moment  method  matrix  Zm  for  the  conical  OCRED. 


Freq  [GHz] 

N 

Fourier  Index,  m 

0 

1 

2 

3 

4 

5 

6 

7 

4 

54 

69 

166 

75 

75 

5.9 

70 

92 

273 

106 

104 

103 

6.9 

80 

110 

249 

129 

125 

124 

124 

10 

108 

161 

353 

180 

190 

181 

178 

178 

177 

Table  7  Condition  number  k  of  N  x  N  moment  method  matrix  Zm  for  the  shallow  cylindrical 
OCRED. 

radius  Rs  (see  Figure  17(a)).  In  effect,  we  consider  the  spherical  OCRED  to  be  a  perturbation  of 
a  spherical  cavity.  The  volume  of  the  corresponding  spherical  cavity  is  Kavity  =  while 

the  volume  of  the  interior  of  the  spherical  OCRED  is  Vocred  =  (|  +  cos  a  —  ^  cos^  a)  ttR^. 
With  a  «  45.6°,  we  find  *=5^  0.939.  A  spherical  cavity  having  a  volume  equal  to  Vocred 
would  have  a  radius  of  approximately  equal  to  0.979i?s,  and  thus  would  resonate  at  approximately 
5.89/0.979=6.015  GHz.  From  this  we  see  that  the  condition  number  of  the  m  =  0  Fourier  mode 
peaks  at  a  frequency  corresponding  to  the  resonant  frequency  of  a  spherical  cavity  having  the  same 
volume  as  the  spherical  OCRED  interior.  We  expect  this  relation  to  be  valid  when  the  aperture 
angle  a  is  small  enough  such  that  the  fields  inside  the  OCRED  are  similar  to  the  fields  inside  the 
corresponding  spherical  cavity.  As  the  aperture  angle  a  is  reduced  toward  zero,  we  would  expect 
the  peak  in  Figure  61  to  become  higher  and  narrower,  and  to  move  toward  5.89  GHz. 

5.4.3  Relative  Importance  of  Fourier  Modes.  The  accuracy  of  the  RCS  produced  by  the 
numerical  algorithm  is  a  function  of  the  number  of  Fourier  modes  included  in  the  solution.  It  is  not 
known  a  priori  how  many  Fourier  modes  must  be  included  to  achieve  a  desired  level  of  accuracy, 
but  a  rule-of-thumb  [41,44]  is  to  include  all  Fourier  modes  up  to  the  M*,  where  M  =  \2kpm!oi\, 
k  is  the  wavenumber,  p„ax  is  the  maximum  distance  between  a  point  on  the  generating  arc  and  the 
axis  of  rotation  and  [a;]  is  the  least  integer  greater  than  x. 
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Figure  61  Condition  number  for  the  spherical  OCRED  as  a  function  of  frequency,  k  for  the 

m  —  0  Fourier  mode  is  plotted  as - ;  /c  for  the  m  =  1  Fourier  mode  is  plotted  as - ;  k  for  the 

m  =  2  Fourier  mode  is  plotted  as - ;  k  for  the  m  =  3  Fourier  mode  is  plotted  as . . 


We  may  assess  the  relative  importance  of  each  Fourier  mode  in  the  following  manner.  In¬ 
stead  of  looking  at  the  contribution  of  the  Fourier  mode  to  the  total  RCS,  we  will  look  at  its 
contribution  to  the  total  electric  current  density  J(r).  The  importance  metric  we  will  employ  will 
be  the  integral  of  the  m***  Fourier  component  of  |J(r)|,  integrated  along  the  generating  arc.  This 
can  be  expressed  directly  in  terms  of  the  expansion  coefficients  as 


0  0  7Va+i^— 1 


\^n+l  f'Ti—l 
mnl  o 


n=l 

Na 


n=l 


ft — iVo 
Na-\-Ns-\ 

+  E 

n=Na 


•  ^n+l—^n-1  ,  l^(A  <^n+l  —  ^n-1 


(131) 

(132) 


We  will  normalize  this  metric  by  the  contribution  from  the  dominant  Fourier  mode,  such  that 


fi  _ 

iL 

■*m 

P  = 

it 

max  <  1^  > 

t  J  m=0 

(133) 

(134) 


Figures  62-65  show  the  normalized  metric  Im  plotted  versus  Fourier  index  m  for  each  of  the  four 
BOR  OCRED  geometries.  The  aspect  angle  is  held  constant  atO  —  45°,  and  data  is  shown  for  both 
II  and  X  polarizations.  Using  this  metric,  we  see  that  the  contribution  of  the  M*  Fourier  mode 
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Figure  62  Normalized  metric  Im  versus  Fourier  index  for  the  spherical  OCRED,  6  =  45°.  Plots 
for  4  GHz  (solid),  6.25  GHz  (dashed),  and  10  GHz  (dash-dot).  M  =  \2kpmax  \  is  4  at  4  GHz,  6  at 
6.25  GHz,  and  10  at  10  GHz. 


Figure  63  Normalized  metric  Im  versus  Fourier  index  for  the  deep  cylindrical  OCRED,  0  =  45° . 
Plots  for  4  GHz  (solid),  5.9  GHz  (dashed),  6.9  GHz  (dash-dot),  10  GHz  (dotted).  M  =  r2A:pmaxl 
is  3  at  4  GHz,  4  at  5.9  GHz,  5  at  6.9  GHz,  and  7  at  10  GHz. 
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Figure  64  Normalized  metric  Im  versus  Fourier  index  for  the  conical  OCRED,  9  =  45°.  Plots 
for  4  GHz  (solid)  and  10  GHz  (dashed).  M  =  \2kpmax.]  is  5  at  4  GHz  and  11  at  10  GHz. 


Figure  65  Normalized  metric  Im  versus  Fourier  index  for  the  shallow  cylindrical  OCRED,  6  = 
45°.  Plots  for  4  GHz  (solid),  5.9  GHz  (dashed),  6.9  GHz  (dash-dot),  10  GHz  (dotted).  M  = 
[2A:pniaxl  is  3  at  4  GHz,  4  at  5.9  GHz,  5  at  6.9  GHz,  and  7  at  10  GHz. 
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is  several  orders  of  magnitude  below  the  contribution  of  the  dominant  mode.  Furthermore,  for  m 
greater  than  approximately  M/2,  the  integrated  current  magnitude  decays  exponentially  with  m, 
at  a  rate  of  7-8  dB  per  unit  increase  in  m.  Surprisingly,  this  rate  of  decay  is  nearly  independent  of 
geometry  and  of  frequency,  suggesting  that  this  may  be  an  intrinsic  property  of  all  BOR  OCREDs. 
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VI.  Conclusions 


This  dissertation  documents  an  investigation  into  the  development  of  an  integral  equation-based 
solution  for  the  determination  of  the  electromagnetic  scattering  from  a  cavity-backed  aperture  em¬ 
bedded  in  an  infinite,  perfectly  conducting  ground  plane.  The  cavity  may  be  filled  with  a  homo¬ 
geneous,  isotropic,  linear  material.  Such  a  geometry  is  of  interest  in  designing  an  effective  RCS 
enhancement  device  to  be  installed  on  airborne  vehicles  while  not  compromising  aerodynamic  per¬ 
formance  properties.  A  substantial  body  of  work  exists  in  the  open  literature  documenting  research 
into  the  scattering  properties  of  open  cavities,  but  none  incorporates  the  integral  equation  technique 
developed  here.  The  coupled  set  of  integral  equations  developed  here  provides  an  accurate  predic¬ 
tion  for  the  RCS  of  the  open  cavity,  and  the  results  support  the  assertion  that  the  integral  equations 
possess  a  unique  solution  at  all  real,  positive  frequencies. 

In  the  previous  chapters,  a  comprehensive  theoretical,  numerical,  and  experimental  treat¬ 
ment  is  presented  for  the  aforementioned  scattering  problem.  Fundamental  theorems  governing 
the  behavior  of  the  electromagnetic  fields  in  the  cavity  interior  and  upper  half-space  are  introduced 
and  proved.  These  theorems  are  used  to  generate  a  coupled  set  of  three  integral  equations  in  which 
the  unknown  quantities  are  the  tangential  magnetic  field  on  the  cavity  surface  and  the  tangential 
electric  and  magnetic  fields  on  the  aperture.  It  is  shown  that,  for  the  case  of  unfilled  cavities,  the 
formulation  degenerates  into  a  coupled  set  of  two  integral  equations  involving  just  the  tangential 
magnetic  fields.  The  system  excitation  is  provided  by  the  incident  tangential  magnetic  field  on 
the  aperture,  which  may  be  radiated  by  near-field  or  far-field  sources.  Integral  representations  are 
given  for  the  total  fields  in  the  cavity  interior,  and  for  the  scattered  fields  in  the  upper  half-space. 

A  numerical  algorithm,  based  on  the  method  of  moments,  is  presented  to  find  an  approximate 
solution  to  the  set  of  integral  equations  for  the  important  case  of  axisymmetric  (body  of  revolution) 
cavities.  The  unknown  tangential  fields  are  expanded  using  piecewise-linear  functions  and  complex 
exponentials  (Fourier  modes).  Methods  to  treat  the  numerical  computation  of  the  resultant  integrals 
are  given,  as  well  as  a  novel  transformation  to  reduce  the  order  of  the  singularity  of  a  dyadic  kernel. 
Various  aspects  of  the  resultant  approximate  solution  are  discussed.  The  far-field  scattered  fields 
are  compared  to  measurements  of  several  full-scale  physical  models,  and  excellent  agreement  is 
observed.  The  approximate  tangential  fields  on  the  cavity  surface  and  aperture  are  shown  to  be 
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well-behaved  and,  in  particular,  to  satisfy  the  edge  condition  at  the  aperture  rim.  The  stability  of 
the  numerical  solution  is  discussed. 

This  research  represents  the  first  use  of  this  integral  equation-based  technique  to  the  open 
cavity  problem.  The  success  of  the  numerical  algorithm  suggests  several  areas  in  which  future 
research  might  concentrate.  First,  the  numerical  algorithm  could  easily  be  extended  to  handle 
non-axisymmetric  cavities.  This  would  necessitate  using  a  different  set  of  basis  functions;  the 
so-called  “Rao-Wilton-Glisson”  [54]  basis  functions  might  prove  useful  here.  Such  an  extension 
would  enable  the  analysis  of  several  interesting  cavity  shapes,  such  as  rectangular  ones.  Second, 
the  numerical  algorithm  could  be  verified  with  dielectric-filled  cavities.  Such  verification  should 
incorporate  a  measurement  program  with  several  types  of  material  fillings  and  cavity  shapes.  A 
two-dimensional  version  should  be  verified  before  a  three-dimensional  version  is  attempted.  Third, 
an  optimization  scheme  could  be  developed  to  find  the  “best”  open  cavity  shape  for  a  given  set  of 
scattering  requirements.  This  would  provide  an  alternative  to  the  costly  “design-fabricate-measure” 
iterative  loop  that  would  currently  have  to  be  employed  today.  Fourth,  the  integral  equations  could 
be  generalized  to  handle  radiation  from  an  open  cavity  using  an  excitation  source  inside,  or  on  the 
surface  of,  the  cavity.  This  would  allow  the  technique  to  be  used  for  antenna  analysis  as  well  as 
scattering  analysis.  Finally,  the  integral  equation  technique  could  be  used  as  part  of  a  larger,  hybrid 
system  to  handle  installation  issues  in  an  airborne  application.  For  instance,  if  an  open  cavity  were 
embedded  in  the  upper  surface  of  an  aircraft  wing,  the  integral  equation  technique  might  be  used 
for  the  cavity  and  physical  optics  might  be  used  for  the  rest  of  the  airframe. 
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Appendix  A.  Proof  of  the  Fundamental  Theorems 

In  this  appendix,  we  prove  the  fundamental  theorems  presented  in  Chapter  III.  The  proofs  employ 
the  classical  methods  of  infinitesimal  vector  calculus.  For  Theorem  1,  we  also  show  that  r(r,  r ') 
satisfies  a  dyadic  form  of  the  reduced  wave  equation  in  the  distributional  sense. 


Proof  of  Theorem  1 

In  this  section,  we  prove  Theorem  1  using  the  classical  methods  of  infinitesimal  vector  cal¬ 
culus.  We  recognize  that  r(r,  r ')  can  be  expressed  as  the  sum  of  three  components,  namely, 

_  ^  =(£) 

r(r,r')  =  ^r  (?,?')  (135a) 

£=1 

=(£) 

r  {r,f')  =  —jkVG{k;f,r')xeiei  for£=  1,2,3  (135b) 


where  ei  =  x,  §2  =  y,  and  63  =  z.  Now  we  introduce  Green’s  Second  Identity  for  mixed 
vector-dyadic  fields  [61], 


Xl( 


V.V.a)  B-A  (VXV.B 


dv  = 

j  fi  •  j^A  X  X  -I-  X  X  Bj  ds  (136) 


where  F  is  a  volumetric  region  with  regular  boundary  dV,  h  is  the  outward  unit  normal  vector 
on  dV,  and  A  =  A(r)  and  B  =  B(r)  are  arbitrary  vector  and  dyadic  functions,  respectively, 

twice-differentiable  in  V.  In  particular,  we  choose  A(r)  to  satisfy  the  homogeneous  wave  equa- 

—  =(f)  ={£) 

tion  (2)  in  V,  and  let  B(f)  =  F  (r,  r ')  with  £.  and  r'  fixed.  It  can  be  shown  that  F  (r,  r ')  also 

satisfies  equation  (2)  for  r  ^  f'.  With  these  choices,  it  is  easy  to  see  that  the  volume  integral  in 
equation  (136)  vanishes  when  r '  0  V,  where  V  is  the  closure  of  V.  Thus 


A  X  V 


=(^) 
X  F 


-b 


(Vx  a) 


=(«) 
X  F 


ds  =  0 


for  r '  0  F 


(137) 


If  r '  G  F,  we  exclude  from  F  a  small  volume.  Vs,  containing  r'  in  its  interior.  We  choose 
Vs  to  be  a  sphere  of  radius  5  centered  at  r as  shown  in  Figure  66.  We  apply  equation  (136)  to  the 
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Figure  66  A  volume  V  and  volume  V^.  Vg  is  spherical  with  radius  S,  centered  at  r '  6  V. 


region  V  —  Vj.  We  see  that  the  volume  integral  in  equation  (136)  vanishes  and  we  have 


/  “ 

Jav 


A  X  V  X  r  + 


(v  X  a)  X  r 


ds  = 
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Now  we  are  in  the  position  to  evaluate  the  right  side  of  equation  (138)  in  the  limit  as  (5  — 0.  To  do 
so,  we  consider  the  two  integrals  on  the  right  side  of  equation  (138)  separately. 


The  First  Integral.  Consider  the  first  integral  on  the  right  side  of  equation  (138), 


Javs 


A(r)  X 


ds 


(139) 


We  establish  a  local  spherical  coordinate  system  (g,  d,  ip)  centered  at  r Without  loss  of  generality, 
we  let  ^  =  3,  so  that  the  “north  pole”  of  the  local  spherical  coordinate  system  is  aligned  with  63  = 
z,  as  shown  in  Figure  67.  The  spherical  unit  vectors  are:  p  =  sin  cos  ipit  +  sin  d  sin  ipy  +  cos  dz. 


=  cos  cos  (/?x  +  cos  sin  (py  —  sin  i9z,  and  (p  =  —  sin  (px  +  cos  (py.  Using  the  identity  [47] 
a  •  (b  X  E)  =  — b  •  (a  X  c)  =  (a  X  b)  •  c  and  noting  that  n  =  — ^  on  dVs,  we  can  write 


h{S)  ^  X  A(r)] 


=(£) 

VxP  (r,r') 


ds 


(140) 


Since  A(r)  satisfies  the  homogeneous  wave  equation  in  the  neighborhood  of?',  we  have  that  A(r) 
is  infinitely  differentiable  at  r  =  r '  and  thus  can  be  expanded  in  a  Taylor  series  about  f '  in  the 
neighborhood  of  r Indeed,  expressing  A(r)  in  terms  of  its  components  in  Cartesian  coordinates 
as  A{r)  =  Ax{x,y,z)x  +  Ay{x,y,z)y  +  Az (x,  y,  z)z,  we  can  write 
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+  0{6‘^) 


where  the  field  components  and  their  derivatives  are  evaluated  at  r  =  ?'.  Expressed  in  the  same 
=(3) 

coordinates,  the  curl  of  F  (?,?')  is 


_ g — ^  r  1  "I 

VxF  {r,r')  =  —jk-^ — ^  |2  (jA:p+ l)cosi?y+  (jfcp)^  +  1  sini9i9| 


47ry^ 


(142) 


We  apply  these  last  two  equations  to  the  integral  in  equation  (140).  Noting  that  the  differential 
surface  ds  =  5“^  sini9  dd  dp,  we  find 


h{S)=jk 


e-jks 
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(jkd)^  +  jk6  +  1 
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(143) 
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Figure  67  Local  spherical  coordinate  system  centered  at  r ' 


Letting  (5  -4  0,  we  have 


lim/i((5)  =  X  A(f')  -zz 
<5->0  3 


The  Second  Integral.  Consider  the  second  integral  on  the  right  side  of  equation  (138), 


l2{S)  =  -l^v^'  (vx  A(r))  xr^®\f,r') 


ds 


Using  the  same  identities  as  used  in  the  previous  section,  l2{6)  can  be  written  as 


l2iS)  =  -l^v  (vxA(?))-  exr^'^(r,r') 


ds 


where  p  x  \r,  ?')  is 


=(3) 

gxT  (f,r')=jk(jkg+l)- — ^-sinWz 


as  expressed  in  the  local  spherical  coordinates.  Expanding  V  x  A(r)  in  a  Taylor  series  about 
r  =  r',  l2{6)  becomes 

/2(^)  =  $  r  sin^  ■d  d-d  [v'  X  A(r ')]  ■  zz  +  0{6)  (148) 

2  yo  ^ 


Letting  S  ^  0,'we  have 


lim/2(<5)  =  ^V'x  A(r')-zz 
(5^0  3 


By  taking  the  limit  as  Vs  becomes  an  infinitesimally  small  sphere  (that  is,  S  0),  and  incorporat¬ 
ing  equations  (144)  and  (149)  into  equation  (138),  we  can  write 

j  n  •  A  X  X  +  (v  X  a)  X  F^^^  ds  =  jkV  x  A(r ')  •  zz  for  r'  eV  (150) 

=(1)  =(2) 

Similarly,  equation  (150)  holds  for  F  and  F  ,  with  zz  replaced  by  xx  and  yy,  respectively. 
Hence  we  may  write 

J  n  •  A  X  X  F^^^^  +  (v  X  a)  X  F^^^  ds  =  jkV  x  A(r ')  •  for  r'  eV 


where  £=1,2, 3.  By  equation  (135a),  it  follows  from  equation  (151)  that 

J  n  ■  A  X  (v  X  F)  +  (v  X  a)  X  F]  ds  =  jkV  x  A(r')  for  r '  G  F  (152) 
Combining  equations  (137)  and  (152),  we  have 


J  fi  •  |a  X  I^V  X  Fj  -I-  I^V  X  a  X  f|  ds  = 


jkV  X  A(r')  for r'  G  F 
0  for  r '  0  F 


and  the  proof  of  Theorem  1  is  complete.  Notice  that  Theorem  1  yields  no  information  for  the  case 
in  which  r '  G  dV. 
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A  distributional  differential  equation  satisfied  by  F.  In  this  section,  we  show  that  F  satisfies 
the  distributional  differential  equation 


V  X  V  X  F  -  fc^F  =  -jkVS{r-r')  x  I 


First,  let  /a(r ')  be  given  by 

Ia{r')  =  -jkJ  A  -  [V(J(r  -  ?')  x  !]  dv  (155) 

Using  the  identity  a  •  x  =  a  x  b,  the  above  equation  becomes 

/a(r')  =  —jk  f  Ax  V^(r  —  r')  dv  (156) 

Jv 

Noting  V(5(r  —  r ')  =  — V'(5(r  —  r '),  using  the  vector  identity  V'  x  (ua)  =  uV'  x  a  —  a  x  V'u, 
and  noting  V'  x  A(f)  =  0,  we  have 

Ia{r')  =  -Jk  j^Vx  [j(f  -  r ')  a]  dv  (157) 

The  curl  operates  on  primed  variables,  and  can  be  brought  in  front  of  the  integral  sign,  giving  us 


fa(r')  =  -jkV  X  f  S{r-r')A 
Jv 


Exploiting  the  sifting  property  of  the  delta  function,  we  have 


/a(?')  = 


—jkV'  X  A(r ')  for r'  gV 
0  forr'^F 


Now,  let  Ib{r')  be  given  by 


/j,(r')  =  J  A(r)  •  [v  X  V  X  F  -  At^f] 
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Noting  that  A  satisfies  equation  (2), 


V  X  V  X  A(r)  -  k^A{r)  =  0 


(161) 


we  have 


/i(r ')  =  J  [a  •  (v  X  V  X  r)  -  (v  X  V  X  a)  •  r] 


dv 


Using  equation  (136),  we  have 


/^(r ')  =  -  y  n  ■  [a  X  (v  X  r)  +  (v  X  a)  X  r] 


ds 


and  in  light  of  equation  (153) 


Ib{r')  =  { 


-jkVxA{f')  forr'€U 
0  for  r '  ^  F 


(162) 


(163) 


(164) 


Comparing  equations  (159)  and  (164)  we  see  that  lair')  =  h{r')  for  r'  ^  dV  and,  due  to  the 
arbitrary  nature  of  A  leads  to 


V  X  V  X  r  -  =  -jkV6{r-r')  x  I 


(165) 


Proof  of  Theorem  2 

In  this  section,  we  prove  Theorem  2.  We  will  make  use  of  some  of  the  intermediate  results 
of  the  first  proof  of  Theorem  1  as  presented  in  Section  A.  In  particular,  we  will  require  the  vector 
function  A(r)  and  the  dyadic  function  r(r,  r ')  to  satisfy  the  same  requirements  as  in  Section  A. 

Examination  of  equations  (14a)  and  (14b)  along  with  equation  (135b)  reveals  that  we  can 
write  the  auxiliary  dyadic  functions  as 


Let  V  be  regular  volumetric  region  and  Vi  be  its  image  with  respect  to  the  xy-plane.  Let  F  U  Vi  be 
empty;  that  is,  V  lies  completely  on  one  side  of  the  rcy-plane.  Then  we  can  consider  three  cases. 


Case  1:  r '  G  V.  Applying  equation  (151)  in  a  straightforward  manner  to  equation  (166) 
allows  us  to  write 


V  X  Ti  + 


V  X  A 


X  Ti}  ds  =  jkW'  X  A(f') 


forr'eF  (167) 


Case  2:  Fj'  G  V.  Equation  (151)  does  not  apply  directly  to  the  present  case  in  which 
Ti  G  V.  Replacing  z  by  —z  in  the  development  of  equation  (151),  we  can  derive  the  related 
equation 


ds 


=  jkV'i  X  A(ri')  •  eeBi  for  ri  G  V 


(168) 


where 


Vi  X  A  = 


dA^  dA^ 
+ 


dy  dz 


X  + 


dAx  dA, 


dz  dx 


y  + 


dAy  _  dAx 

dx  dy 


(169) 


Applying  equation  (168)  to  equation  (166)  allows  us  to  write 


/ 


■  A  X 


VxEi 


+  VxA  xri|ds  =  jkV'i  X  A(fi')  •  (xx  +  yy  —  zz) 


jk 


V'  X  A(ri') 


(170) 


forfi'  G  F 


Case  3:  r'  ^  F,  Vi  ^  F.  When  neither  r  nor  Fi  lie  in  the  interior  of  F,  then  it  is  clear  from 
equations  (151),  (168),  and  (166)  that  the  surface  integral  on  the  boundary  of  F  vanishes.  That  is. 


J  n-|Ax  V  X  Fi  +  VxAjxri|ds  =  0  for  F' 0  F  and  Fi' 0  F 


(171) 
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Combining  these  three  cases,  we  can  write  the  single  equation, 


V  X  Ti 


+ 


V  X  A 


]xr,} 


ds 


jkV  X  A{r ') 


=  < 


jk  [V'  X  A(r/) 
0 


for r'  eV 
for  r '  G  Vi 

for  r'  and  r '  0  1^ 


and  the  proof  of  Theorem  2  is  complete. 


(172) 


Proof  of  Theorem  3 

The  proof  of  Theorem  3  is  exactly  the  same  as  the  proof  of  Theorem  2  except  that  Ti  is 
replaced  by  r2.  This  leads  to 


/ 

JdV  '■ 


A  X 


VxTa 


+ 


V  X  A 


%} 


X  To  S  ds 


=  < 


jkV  X  A(f  0 
-jk  [V'  X  A(r/)] . 


for  f'  eV 
forr'  G  Vi 

for  r '  0  F  and  r '  0  F 


(173) 
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Appendix  B.  Application  of  the  Fundamental  Theorems 

In  this  appendix,  we  present  the  results  of  applying  Theorems  1, 2,  and  3  to  the  OCRED  geometry. 
Each  application  represents  a  choice  of  one  from  each  of  the  following  three  sets: 

•  Theorem  1,  Theorem  2,  or  Theorem  3. 

•  E  is  the  cavity  volume,  D,  or  the  upper  half  space  (UHS),  2:  >  0. 

•  A  is  the  electric  or  magnetic  field. 

The  table  at  the  bottom  of  this  page  provides  an  index  and  gives  the  equation  number  correspond¬ 
ing  to  each  combination.  Significant  algebraic  manipulations  are  needed  to  arrive  at  the  equations 
that  follow,  and  will  not  be  given  here.  The  interested  reader  can  verify  any  equation  by  exploiting 
vector  and  dyadic  identities  (all  of  which  can  be  found  in  [62]),  along  with  the  boundary  condi¬ 
tions  satisfied  by  E  on  5  and  and  the  boundary  conditions  satisfied  by  Fi  and  r2  on  a  (see 
equations  (15)-(16)). 


Theorem  1,  V  =  D,  A  =  E.- 


/ 


z  X  E(r)]  •  V  X  r(A:i,r, r')  da  —  jkiZi 


/ 


X  H(r)  •  r(/:i,r, r')  ds 


{A:?ZiH(r')  forf'eD 

(174) 

0  for  >  0 


Table  8  Equation  index  for  application  of  the  fundamental  theorems  to  the  OCRED  problem. 


Theorem  1 

Theorem  2 

Theorem  3 

E  =  D,  A  =  E 

(174) 

(178) 

(182) 

F  =  D,  A  =  H 

(175) 

(183) 

(176) 

■ikbM 

(184) 

(177) 

(181) 

(185) 
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Theorem  1,  V  =  D,  A  =  H; 


f  z  X  E(r)  ■r{ki,r,r')  da  +  ^  f  fn  x  H(f)  -V  xr{ki,r,r')  ds 

J <7  ^  jkiri  J SUa  ^  ■' 


jki'E{r')  forr'eD 

0  for  z'  >0 


Theorem  1,V  =  UHS,  A  =  E; 


/  z  X  E(f)l  •  V  X  r{ko,  r,  r ')  da  -  jkoZo  [  z  x  H*(r)l  •  T{ko,r,  r')  da 
J<T  ^  -*  JaUcr<^  ^ 


0  for  r'  e  D 

-k‘^ZoW{r')  forz'>0 


Theorem  1,V  =  UHS,  A  =  H; 


[  z  X  E(r)  ■r{ko,r,r')  da  + —^  [  [z  x  H®{r)l  ■  V  x  r(fco,  r,  ?')  dcr 
Ja  ^  JI^O^O  JaU(T<^  ^ 


I  for r'  E  D 

-jkoW{r')  for2r'>0 
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Theorem  2,  V  =  D,  A  =  E; 


2  y  [z  X  E(r)]  •  V  X  r(A:i,  r,  r ')  da  -  jk^Zi  J  [n  x  H(r)]  •  ri(fci,  r,  r ')  ds 

klZiH{r')  forr'eD 

fc^ZiHi(rj')  forr'GA  (^^8) 

0  for  r'  E  Df 

Theorem  2,  V  =  D,  A  =  H; 


2  J  z  X  H(f)  •  V  X  r(A:i ,  r,  f')  c?(T  +  J  j^n  x  H(r)  •  V  x  Fi  (A:i ,  f,  f ')  ds 

-k‘fYiE{r')  forr'eD 

-k^YiEii^i)  forr'eA  (179) 

0  forr'eD/ 

Theorem  2,V  =  UHS,  A  =  E^; 


Theorem  2,V  =  UHS,  A  =  H®.- 


/ 


z  X  H*(?) 

j  •  V  X  r(fco,?,?')  da  =  < 

k 

for  r'  £  D 
for  2:'  >  0 


Theorem  3,  V  =  D,  A  =  E; 


2  J  j^z  X  H(r)  •  r{ki,r,f')  da  +  j  n  x  H(r)  •  r2(fei,  r,  ?')  ds 


=  < 


3hn{v') 

-jkiRiin') 

0 


for  r'  e  D 
forr'  e  Di 
for?'  e  Df 


Theorem  3,  V  =  D,  A  —  H; 


2J  zxE(r)  •r(A;i,f,r')d(7+-^^  [nxH(r)]  ■Vxr2(A:i,r,r')ds 


jfciE(r')  forr'eD 

=  \  -jkiEi{Ti')  for?'  6  Di 

0  for?'  G  Df 


(181) 


(182) 


(183) 
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Theorem  3,V  =  UHS,  A  =  W: 


p  I  J  'J  ^ 

/  [zxH^(r)1  ■F(A:o,r,r')da=  i  ^ 

J  (TUcr’^  ^  ■'  I  ~.lko 


^H?(ri')  forr'eD 


=^W{r')  forz'>0 


Theorem  3,  V  =  UHS,  A  =  H^- 


f  [z  X  E(f)]  ■r(A:o,r,rO  da  = 


^Ef(ri')  forr'sD 

^W{r')  forz'>0 
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Appendix  C.  Regularization  of  a  Singular  Integral 

When  formulating  the  expressions  for  the  elements  of  and  Z^,  we  must  deal  with  an  iterated 
surface  integral  of  the  form 


VVG(fc;r,r')]  ds  ds' 


(186) 


where  5i  and  S2  are  regular  surfaces,  A(r ')  is  defined  on  B(r)  is  defined  on  S2,  n'  is  the  unit 
normal  vector  on  Si  at  the  point  defined  by  r and  n  has  a  similar  relationship  with  S2  and  r.  The 
unit  normal  vectors  vary  continuously  almost  everywhere  on  their  domains.  We  further  suppose 
that  n'  ■  A(r')  =  0  on  and  that  n  •  B(f)  =  0  on  82',  that  is,  A  and  B  have  only  tangential 
components  on  their  respective  domains.  The  kernel  of  equation  (186)  is  VVG(A:;  r,  ?'),  which 
can  be  written  as 


-jkR 


VVG(A:;  r,  f ')  =  {  [3  +  3jA:i?  +  [jkRf  RR  -  (1  +  jkR)  l}  (187) 


where  R=  |r  —  r'|  and  R  =  (r  —  f ')  /R.  From  equation  (187)  we  see  that 


VVG(fc;  f,  r ')  ~  (3RR  -  l)  for 


r  r 


(188) 


Let  i?min  be  the  minimum  value  of  R  for  any  f '  6  Si  and  r  G  82-  If  Si  and  S2  do  not 
overlap,  then  Rmin  >  0,  VVG(A:;  f,r')  is  bounded,  and  the  integral  in  equation  (186)  exists  and 
is  uniformly  convergent.  However,  if  Rmin  is  small,  then  the  absolute  value  of  the  kernel  of  equa¬ 
tion  (186)  will  be  C?(Rj^^)  when  r'  is  near  r,  and  numerical  integration  techniques  may  have  poor 
convergence  properties.  If  Si  cuid  S2  overlap,  then  the  integral  in  equation  (186)  must  be  treated 
very  carefully.  We  may  consider  the  overlapping  case  to  be  the  limit  of  a  non-overlapping  case 
as  Rmin  — >■  0.  Wilson  [65]  establishes  the  existence  of  this  integral,  but  its  numerical  evaluation 
is  cumbersome  and  computationally  expensive.  We  now  show  a  method  by  which  equation  (186) 
may  be  transformed  into  a  form  that  is  much  more  amenable  to  numerical  computation. 

Let  Si  and  S2  be  two  regular  surfaces  that  do  not  overlap,  as  shown  in  Figure  68.  Also,  let 
A(r ')  and  B(r)  be  continuous,  differentiable  vector  functions  defined  on  Si  and  S2,  respectively. 
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Figure  68  Geometry  used  in  regularization  procedure.  Regular  surfaces  and  S2  do  not  over¬ 
lap.  r  G  iS'2  and  fi  is  the  unit  normal  on  S2  at  f.  r '  G  and  n'  is  the  unit  normal  on  5i  at 


Because  the  surface  integrals  in  equation  (186)  are  uniformly  convergent,  we  may  write 

/  =  /  A(f ')  •  n' X  C(r ')  (189a) 

Jsi 


where 


C(r')  =  [  B(r)  •  VVG(A:;r,r')  ds  (189b) 

JS2 

Since  B(f)  is  everywhere  timgent  to  S2,  then  B(r)  =  n  x  B(r)  x  n  so  that 

C(r ')  =  y*  {n  X  [B(r)  x  n]  }  •  VVG(A:;  r,  f ')  ds  (190) 

Using  the  identity  x  bj  ■  c  =  a  •  x  cj  with  a  =  n,  b  =  B(r)  x  n  and  c  =  VVG(fc;  r,  r '), 
equation  (190)  becomes 

C{r')  =  j  n-{[B(r)xn]  X  VVG(A:;r,f')}  ds  (191) 
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Using  another  identity,  namely  V  x  =  (V  x  a)  b  —  a  x  Vb  with  a  =  B(r)  x  n  and 

b  =  VG{k;  r,  r '),  equation  (191)  becomes 


C(r')  =  f  n 

JS2 


•  V  X 


B(r)  X  nj  VG(fc;r,r')  ds 

-  j  n  •  V  X  I  [B(f)  X  fi]  VG(fc;  r,  r ')}  ds  (192) 

So 


We  introduce  a  dyadic  form  of  Stokes’  Theorem: 


/nVxads=®  dca 

Js  JdS 


(193) 


where  the  line  integral  is  traversed  in  the  positive  sense  with  respect  to  the  unit  normal  vector  n 
(i.e.,  in  accordance  with  the  right-hand  rule).  Applying  equation  (193)  to  the  second  integral  in 
equation  (192),  we  have 

C(f')  =  J  n  ■  V  X  [B(f)  X  n]  ^G{k]r,r')  ds 

-  I  VG{t,f,f')  [B(r)  X  nl  •  dc  (194) 

JdSi  *■ 


The  first  integral  in  equation  (194)  may  be  simplified  by  noting  that 


n  ■  V  X  B(r)  x  n  =  — V*  •  B(r) 


(195) 


where  “Vj-”  is  the  surface  divergence  operator  [62,  page  501],  Then  equation  (194)  becomes 


C(r')  =  -  [  V,  •  B(r)  VG(A:;r,r')  ds 

JS2 


l  VG(A;;r,f')  [B(r)  X  nl  -dc  (196) 

JdS2  ^  ■' 
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Combining  equation  (196)  with  equation  (189a),  and  noting  VG{k-,r,r')  =  — VG(fe;r,  r'),  we 
obtain 


1=  [  [  A(r ')  •  [n'  X  VG{k- r,  r ')]  V,  •  B(r)  ds  ds' 

J  Si  J  S2 

+  [  (f  A(f')  •  [n'  X  V'G(A:;r,r')]  [B(f)  x  n]  ■  dc  ds'  (197) 

Js,  Jasi  ^ 


or 


1=  f  D{r)  Vs  •  B(r)  ds  +  (f  i7(r)  : 

J  S2  J  ^S2 


Bfr)  X  n 


dc 


(198a) 


with 


£)(?)=  f  A(r')-n'x  V'G(fc;r,r')ds' 

Jsi 

Using  the  identity  a  •  x  c j  =  b  •  (c  x  a),  equation  (198b)  becomes 

D{t)=  f  n'- [v'G(A;;f,fO  X  A(r')l  ds' 
Jsi  •- 


(mb) 


(199) 


Using  another  vector  identity,  namely  V'  x  {<pa)  =  V'^  x  a  +  (V'  x  a.)  cj)  with  0  =  G{k-,  r,  r ') 
and  a  =  A(f '),  equation  (199)  becomes 

D{r)  =  J  h' -V' X  j^G(A:;f,  f')A(r')j  ds' 

-  [  h'-V'  X  A(r ')  G{k- r,  f ')  ds'  (200) 

Js, 


Using  a  vector  form  of  Stokes’  Theorem, 


equation  (200)  may  be  written  as 


£)(?)=/  G{k;r,T')  A{f')  ■  dc' —  f  h'  -V'  x  A{r')  G{k-,r,r')  ds' 

Jasi  Jsi 


(202) 


The  second  integral  in  equation  (202)  may  be  simplified  by  noting  that  A(r ')  is  everywhere  tangent 
to  Si,  so  that  A(r ')  =  n'  x  A(r')  x  n'.  Then,  incorporating  equation  (195),  equation  (202) 
becomes 


D{r)  =  i  G{k-,  f,  r ')  A(f ')  ■  dc' +  f  G{k;  r,  r ')  V,'  •  [n'  x  A(r ')]  ds'  (203) 
JdSi  J  Si 


Finally,  we  combine  equation  (203)  with  equation  (198a)  to  get 


(204) 


This  proves  Theorem  5. 

If  we  compare  equation  (204)  with  equation  (186),  we  see  that  the  order  of  the  singularity 
of  the  integrand  has  been  reduced  from  to  R~^.  This  has  been  accomplished  at  the  cost  of 
requiring  that  the  vector  functions  A  and  B  have  bounded  derivatives  everywhere  on  their  domains 
and  evaluating  additional  integrals  on  the  boundaries  of  and  82-  These  additional  integrals  need 
not  always  be  numerically  evaluated.  They  must  be  considered  only  if  A(r')  has  a  non-vanishing 
tangential  component  on  dSi,  or  if  B(f)  x  n  has  a  non- vanishing  tangential  component  on  882. 

The  advantage  of  equation  (204)  over  equation  (186)  can  be  illustrated  by  the  following 
numerical  example.  Let  5i  be  a  lA-radius  disk  in  the  xy-plane,  centered  at  the  origin.  Let  82  be  a 
lA-radius  disk  parallel  to  Si,  centered  at  (0, 0,  ^),  ^  >  0.  This  geometry  is  shown  in  Figure  69. 
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Figure  69  Geometry  for  numerical  example,  is  a  lA-radius  disk  centered  at  (0, 0, 0).  S2  is 
a  lA-radius  disk  centered  at  (0, 0,  (5).  Both  are  parallel  to  the  3:y-plane.  The  vectors  are  such  that 
r '  G  Si  and  r  €  S^- 

Define  A(f ')  and  B(r)  as 


II 

1— * 

1 

1 

(205a) 

B(f')  =  ^ 

(205b) 

P 


with  p'  =  |r'|,  p  =  |r  —  Sz\,  =  — xsin</>'  +  ycos^',  p  =  xcos0  +  y sinc/i.  Since  both 
A(r ')  and  B(r)  vanish  on  the  boundaries  of  their  respective  domains,  only  the  first  integral  in 
equation  (204)  survives.  The  surface  divergences  are 


V;  •  [n'  X  A(r')]  =  (206a) 

Vs-B(r')  =  ±-e''^  (206b) 

P 

where  the  positive  sign  is  taken  when  p'  or  p  is  less  than  0.5A  and  the  negative  sign  is  taken  other¬ 
wise.  A  computer  program  was  written  to  calculate  I  using  the  formulation  of  equation  (186)  and 
of  equation  (204).  Figure  70  shows  the  results  for  various  values  of  the  separation  distance  6  and 
the  Gauss-Legendre  quadrature  order  Ng.  With  6  held  constant  the  formulation  of  equation  (186) 
converges  very  slowly  as  Ng  increases,  and  with  Ng  held  constant  it  diverges  as  S  decreases.  In 
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Ng  l/S 


Figure  70  Magnitude  of  I  for  as  a  function  of  quadrature  order  Ng  with  5  =  0.02,  and  as  a 
function  of  6  with  Ng  =  30.  Equation  (186)  formulation  is  plotted  as  a  dashed  line;  equation  (204) 
formulation  is  plotted  as  a  solid  line. 

contrast,  the  formulation  of  equation  (204)  is  quite  well-behaved.  Thus,  equation  (204)  requires 
far  fewer  computational  resources  to  evaluate  than  equation  (186). 
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